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Epigraph 



We are at the very beginning of time for the human race. It is 
not unreasonable that we grapple with problems. But there are 
tens of thousands of years in the future. Our responsibility is to 
do what we can, learn what we can, improve the solutions, and 
pass them on. 

- Richard Feynman 
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Figure 6.14: Simulated wakefield and energy spectra for the experimental data of the 
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Abstract 



Particle accelerators are the tools that physicists use today in order to probe the fundamental 
forces of Nature, by accelerating charged particles such as electrons and protons to high 
energies and then smashing them together. For the past 70 years the acceleration schemes have 
been based on the same technology, which is to place the particles onto radio-frequency electric 
fields inside metallic cavities. However, since the accelerating gradients cannot be increased 
arbitrarily due to limiting effects such as wall breakdown, in order to reach higher energies 
today's accelerators require km-long structures that have become very expensive to built, and 
therefore novel accelerating techniques are needed to push the energy frontier further. 

Plasmas do not suffer from those limitations since they are gases that are already 
broken down into electrons and ions. In addition, the collective behavior of the particles in 
plasmas allows for generated accelerating electric fields that are orders of magnitude larger 
than those available in conventional accelerators. Such wakefields have been demonstrated 
experimentally, typically by feeding either single electron bunches or laser beams into high 
density plasmas. As such plasma acceleration technologies mature, one of the main future 
challenges is to monoenergetically accelerate a second trailing bunch by multiplying its energy in 
an efficient manner, so that it can potentially be used in a future particle collider. 

The work presented in this dissertation is a fruitful combination of theory, simulations 
and experiments that analyzes the use of multiple electron bunches in order to enhance certain 
plasma acceleration schemes. Specifically, the acceleration of a trailing electron bunch in a high- 
gradient Wakefield driven by a preceding bunch is demonstrated experimentally for the first 
time by using bunches short enough to sample a small phase of the plasma wakes. Additionally, 

xxi 



it is found through theoretical analysis and through simulations that by using multiple bunches 
to drive the wakefields, the energy of a trailing bunch could be efficiently multiplied in a single 
stage, thus possibly reducing the total length of the accelerator to a more manageable scale. 
Relevant proof-of-principle experimental results are also presented, along with suggested 
designs that could be tested in the near future. Furthermore, electron beam and plasma 
diagnostics are analyzed and presented, which are necessary for properly completing and 
understanding any plasma wakefield experiment. Finally, certain types of plasma sources that 
can be used in related experiments are designed, diagnosed and tested in detail. 
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In science one tries to tell people, in such a way as to be 
understood by everyone, something that no one ever knew 
before. But in poetry, it's the exact opposite. 

- Paul Dirac 



Chapter 1 
Introduction 

1-1 Physics and Particle Accelerators 

"\Ne are trying to prove ourselves wrong as quickly as possible, because only in that way can we 
find progress". This is how the celebrated physicist Richard Feynman described the evolution of 
physics in his famous series of lectures he gave at Cornell University in 1965 regarding the 
character of physical laws [60]. Indeed, the experimenters are not very interested in confirming 
already-established laws under regular conditions. Rather, they prefer to look carefully in those 
places where the laws are most likely to be falsified. Consider for example the experiments that 
measured the blackbody spectrum [151] and the bending of light near the sun [46] in the early 
20^'' century which challenged our understanding of nature, thus allowing Max Planck and Albert 
Einstein, respectively, to suggest that energy is quantized [142] and that space is not flat [48], 
ultimately giving birth to the modern theories of Quantum Mechanics and General Relativity. 

Today and for the past 70 or so years, the most successful tools that we use to look into 
those new regimes are particle accelerators [138]. These machines use strong electric fields to 
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accelerate charged particles such as electrons and protons to high energies which are then 
smashed against each other at near the speed of light\ By observing the behavior of the 
products of those collisions, the elementary constituents of matter can be analyzed and new 
fundamental laws can be formed [11]. The ongoing success of these colliders has relied on the 
technological developments that allow them to accelerate particles at ever higher energies: 
From Einstein's famous mass-energy equivalence equation, E = mc^ , higher energies translate 
to more massive particles that can be created and detected in the accelerators. The progress of 
the accelerator energy over the past decades is depicted in Figure 1.1 . 

The energy of accelerators has increased almost 10 orders of magnitude since the first 
cyclotrons where constructed in the early 1930s, and this scaling has produced some amazing 
discoveries in the field of particle physics as the higher energies allowed for heavier particles to 



^ If a proton from the 14 TeV Large Hadron Collider was to race a photon around the earth, it would finish 
the race in 145 |is, lagging only 1 mm behind the photon. 



be detected. Most of the quark structure inside protons and neutrons was identified in the 
1970s and 1980s at the Stanford Linear Accelerator Center (SLAC) [8], with masses that range 
between 1 MeV/c^ (for the up quark) to 5 GeV/c^ (for the bottom quark). The carriers of the 
electroweak force, W and Z bosons, were discovered in 1983 at CERN's Super Proton 
Synchrotron (SPS) and were found to have masses 80 GeV/c^ and 90 GeV/c^ respectively. The 
heaviest elementary particle known to date, the top quark, was discovered only in 1995 at 
Fermilab's Tevatron accelerator [39] with a mass of 171 GeV/c^. Without particle accelerators, 
our understanding of the four fundamental forces of nature would have been limited to 
Einstein's theory of gravity [48] from 1915 and Maxwell's equations for electromagnetism [121] 
from 1865. 

This success though has not happened without cost, as accelerators have come a long 
way since Lawrence's first 1 MeV cyclotron in 1932 which was a mere 10 cm long [103]. The 
price, complexity and size of the latest particle accelerators have skyrocketed along with their 
energy. This dissertation is written in 2008, the year that the Large Hadron Collider (LHC) is 
commissioned. The LHC is the largest and most complicated machine ever constructed by 
humans, a circular 14 TeV proton - proton collider built underground in a 27 km tunnel at 
CERN [1], near the French-Swiss borders^ This worldwide effort of 2,000 scientists from 34 
countries will look for, among other things, the Higgs Boson with an upper expected mass of 
144GeV/c^ (which may provide the tools to unify the three of the four fundamental forces), 
evidence for supersymmetric particles (which may explain the excess of non-baryonic matter in 



^ 14 TeV is roughly the energy needed to pick up a piece of hair from the floor and place it on a table (this 
is mostly kinetic energy - the rest mass energy of a proton is E = mc^ -IGeV). The total energy 
circulating in the 2x2800 bunches of LHC with 10^^ protons each is 724 MJ or the energy of a TGV train 
running at full speed (320 km/h). 



the universe) and string theory's extra dimensions (which may explain why gravity is so weak 
compared to the other forces in nature). The price tag for LHC is between $5 - 10 billion. 

LHC - scale accelerators are very close to the limit of what we can practically afford to 
build using conventional technologies, even collaboratively. The International Linear Collider 
(ILC) [3], a proposed 1 TeV lepton smasher, is estimated to cost $7 billion and extend over 30 km 
in length, more than 10 times the length of SLAC. The Compact Linear Collider (CLIC) [175] is 
another scenario of a proposed lepton collider at 3 TeV with a length of 33 km and similar cost. 
Any such projects will be extremely hard to finance and may ultimately follow the path of the 
Superconducting Super Collider (SSC), a 87 km circumference accelerator that was cancelled 
after the projected cost exceeded $12 billion in 1993\ 

The reason why the cost has been scaling with the collider energy is that the accelerating 
gradients (i.e., the energy gained per unit length) have more or less remained constant over the 
past few decades, in the order of 10 - 100 MV/m. Therefore, the only way to scale into higher 
energies is to simply make the accelerating portion longer, thus increasing the construction and 
maintenance costs at the same time. Those conventional techniques rely on radio-frequency 
electric fields that propagate inside metallic cavities. One of the major limitations of this 
technology is that if the accelerating fields are increased beyond a certain point, the modules 
are damaged [23] because the energy from these fields is approaching the energy that binds 
electrons to the atoms in the cavity walls and tears them apart. This is why the ILC gradient can 
only be 40 MV/m [3] (limited by the loss of superconductivity in the niobium cavities due to 
heat), only two times larger than the 20 MV/m of the 40 year old SLAC. Clearly, a revolution is 



^ 24 km of the underground tunnel were bored in Texas after 2 years of construction before cancellation. 
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needed here; a new acceleration technology that will do for the accelerators what the transistor 
did for the electronics 60 years ago: miniaturize them. 

1,2 Plasma Accelerators 

1.2.1 Basic Principles 

The new particle acceleration technology that is the subject of this dissertation is plasma 
accelerators [50, 51]. The main principle behind it is that one can use a relativistic electron beam 
(Plasma Wakefield Acceleration - PWFA) or a laser pulse (Laser Wakefield Acceleration - LWFA) 
to excite relativistically moving waves in plasmas (ionized gases). The electric fields excited in 
those waves are at least 2-3 orders of magnitude stronger than the ones supported by the 
conventional metallic cavities, and thus, when utilized properly, they are used to accelerate 
particles to high energies in much smaller distances (thus reducing the cost). In addition, since 
plasma accelerating cavities (which are created dynamically every time) consist of already 
ionized particles, they do not suffer from the material breakdown damage like the conventional 
metallic cavities. Figure 1.2 illustrates the basic physical principles of a plasma accelerator. 

As the relativistic electron beam or the laser beam propagates inside the plasma, the 
electric field of the beam displaces the plasma electrons^ while the plasma ions are much 
heavier and stay static. Eventually the electrons are drawn back on the axis from the force of the 
positively charged ions, overshoot the axis, and they start to oscillate (perpendicular to the 
propagation direction of the beam in this figure). This sets up a charge density wave behind the 
beam, since the regions where the electrons cross the axis are more negatively charged than the 
regions where the electrons are far away from the axis. This charge difference creates a very 



^ The electric field of a relativistic beam is mostly transverse to its direction of motion. 
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Figure 1.2: The main physical picture of a plasma accelerator. A relativistically moving electron beam or 
laser pulse is moving to the right inside a uniform plasma. The electric field of the driver displaces the light 
plasma electrons, which are then drawn back on axis from the static ions, setting up an electron density 
oscillation. The corresponding electric field that is generated on axis is plotted in the bottom. 



strong electric field that propagates at the velocity of the beam. As a result, if a witness electron 
beam is placed at the proper phase, it will experience this accelerating field and be propelled to 
high energies. 

A simple estimate for the strength of the electric field oscillation amplitude can be made 
by using Gauss's law: 



W.E=P- 



(1) 



+ik,^z 



Assuming a one - dimensional plane wave perturbation of the charge density, p-etSne ' , 
where e is the electron charge, A/z is the plasma density perturbation, and k =co I c is the 
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plasma wavenumber, then the electric field perturbation is E^zE^e ' , with W = zdl dz (cis 



the speed of light in vacuum). For a density perturbation on the order of the neutral plasma 
density, A^z ~ fz^, the electric field amplitude E^ is roughly 
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Note that the plasma frequency co^ is given by 



%-i^ (3) 



For a plasma with a neutral density n^ =10 cm , equation (2) indicates that an accelerating 
electric field (wakefield) of 10 GV/m could be supported by the plasma. These extremely large 
wakefields are the single most important reason that has motivated the plasma accelerator 
community over the last 30 years to investigate the details of using plasmas as accelerating 
structures. 

The physical reason why plasmas can support those high gradients is the collective 
effect of the plasma electrons. As opposed to the other three states of matter where the 
particles and molecules are distributed in a fuzzy and incoherent way, inside a plasma its billions 
of freed electrons can be manipulated together and forced to act coherently. The success of the 
future plasma accelerator will rely on maintaining this nice collective behavior over long enough 
distances. Moreover, it should be pointed out that the plasma accelerator is merely an energy 
transformer. It does not provide energy; it may only transfer the energy of an existing beam to a 
trailing beam. As such, it still requires a pre-created driver to provide the initial energy to be 
extracted. Therefore, the plasma accelerator will not replace existing accelerators; it will rather 
extend them to higher energies [81]. 



1.2.2 Progress of Plasma Accelerators 

In order to put the work presented in this dissertation into perspective, a review of the history 
of plasma accelerators until today is necessary [42] and the purpose of this section. 

Plasmas where identified as early as 1879 [40] (dubbed radiant matter), by Langmuir's 
time in the 1920s [102] they were known to support propagation of electron waves [169] and 
became well understood by the 1950s [5, 41]. Fainberg in 1956 was the first to suggest [56] that 
if one would generate plasma waves that propagated at relativistic speeds ( v^^ - c ) then 

particles could be accelerated by sampling the relativistic electric fields inside the plasma. 
Although some first demonstrations of the interaction of an electron bunch train with a plasma 
happened in the mid-1970s by Berezin et al. [17, 57] showing 250 kV/m gradients, it wasn't 
until the rigorous work of Tajima and Dawson in their seminal 1979 paper [168] that a significant 
effort to conduct research on plasma accelerators was seriously initiated [43, 82]. In this paper 
the idea of using a laser to excite the relativistic waves in the plasma was first proposed. By 1985 
the excitation of the first GV/m unloaded relativistic waves in plasmas was reported by Clayton 
et al. [36] using laser beat - wave techniques^ 

Earlier in the same year Chen et al. [31, 93] analyzed the excitation of plasma waves by 
relativistic electron beams instead of lasers. In 1988 Rosenzweig et al. demonstrated 
experimentally for the first time the interaction of a trailing electron beam with the wakefield 
excited by a preceding drive electron beam [148]. They fed a 21 MeV drive beam with 2 - 3 nC 
charge which was 7 ps long into a 33 cm long plasma with a neutral density of 10^^ cm'^ and 
observed 50 keV energy shifts of the witness beam at 5 MV/m gradients. Those experiments 
were later extended in the mildly nonlinear regime in 1989 [14, 149] (by allowing the plasma 



^ Single laser pulses were not intense enough at the time to excite those wakes, thus the beating between 
two laser pulses was utilized. 



density n^ to become smaller than the beam density n^), as was analyzed theoretically in 1987 
also by Rosenzweig [146]. 

By 1991 the so-called blowout or bubble regime of the plasma wakefield accelerators 
was identified through simulations [13, 147]. In this highly nonlinear regime the drive beam is 
much denser than the plasma, n^ »^o' ^"^^ ^^ ^ result the region behind the driver is void of 
electrons leaving a uniform density ion column [15, 113, 114, 115]. A witness beam placed in 
that region was found to experience simultaneously uniform accelerating and focusing forces 
along its transverse size, an extremely desirable situation that is not possible with linear drivers. 

The first evidence of accelerated electrons using a laser-driven plasma were reported by 
Kitagawa et al. in 1992 [99] by trapping electrons in a beat-wave setup. Soon after, in 1993, in a 
series of very well monitored experiments, Clayton et al. showed the acceleration of externally 
injected electrons in a laser beat-wave plasma with 0.7 GV/m gradients [37], which was later 
followed by acceleration of trapped electrons up to 28 MeV by sampling 2.8 GV/m gradients, as 
reported by Everett et al. in 1994 [55]. By 1995 Modena et al. [122] reported the use of an 
intense 5x10^^ W/cm^ laser to drive the plasma to wave breaking by inducing the Raman 
Forward Scattering instability, thus producing the most energetic electrons from a laser-driven 
accelerator at the time (44 MeV)\ Still, the laser pulse needed to be many times longer than the 
plasma wavelength for these schemes to work. 

As lasers became more powerful in the early 2000s reaching multi - TW levels, it 
became possible to excite plasma waves in the blowout regime [112] by using a single short 
pulse (with length comparable to the plasma wavelength). In 2002, Malka et al. published the 
first such results [119] by accelerating trapped electrons to 200 MeV using 1 Joule energy 30 fs 



The LWFA gradients where much higher than the PWFA gradients at the time due to the difference in 
the driver intensities: typically 10 W/cm for laser beams, 10 W/cm for particle beams. 



long TitaniumiSapphire laser pulses onto 2x10^^ cm"^ density gas jets, albeit with 100% electron 
energy spread at the exit. Soon it was clear from the first 3D LWFA simulations [117, 145, 170] 
that in this regime very low energy spread bunches could be produced; it wasn't long before 
three groups independently reached this "sweet spot" in 2004 [58, 63, 120] by trapping and 
accelerating 0.5 nC bunches to over 100 MeV with less than 10% energy spread. More than 20 
groups worldwide have since repeated those results, and in 2006 the output energy was 
increased to 1 GeV by propagating similar laser pulses inside a cm-scale capillary plasma by 
Leemans et al. [106, 127]. This last record-breaking result for laser-driven accelerators also 
brought forth the main issue with these schemes, the fact that laser pulses tend not to 
propagate focused over long distances due to diffraction, and it is challenging to extend these 
acceleration techniques to meter-scale plasmas that are required for particle-physics-quality 
energies. 

Relativistic charged particle beam drivers on the other hand are "stiffer" and do not 
suffer from diffraction as severely^ While the laser-driven acceleration was flourishing, the 
PWFA community had set forth to demonstrate that plasma accelerators can work over meter- 
scales by using the 50 GeV SLAC electron beam as a driver in the ultrarelativistic blowout 
regime^ [73]. This collaboration produced a wide array of wonderful experimental results in the 
2000s, such as the transverse beam-envelope dynamics by Clayton et al. in 2002 [35], the study 
of positron focusing by Hogan et al. [75] and the 56 MV/m positron acceleration in positron- 
driven wakefields by Blue et al. [19] in 2003, the 150 MV/m acceleration of the electron beam 
by Muggli et al. in 2004 [123], the first multi-GeV energy gain for any plasma accelerator by 



^ The effective relativistic parameter / for a laser \s -cd^/ cd ~ 10 for a 1 pinn laser at 10^^ cm"^ density, 

while the SLAC beam has ^=^10^ 

^ This was initially proposed by John Dawson. 
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Hogan et al. in 2005 [74], and the production of ultra-short trapped electron bunches by Oz et 
al. in 2007 [136], culminating with the doubling of the energy for some of the electrons of the 
beam from 42 GeV to 84 GeV in a little less than one meter of plasma as reported by Blumenfeld 
et al., also in 2007 [18, 20]. Although the energy spread of those electrons was an undesirable 
100%, this latest remarkable result proved that over 50GV/m accelerating gradients can be 
maintained in long plasmas. 

1-3 Present Challenges for PWFAs and Scope of the Dissertation 

In this section we will analyze the upcoming challenges that plasma wakefield accelerators have 
to surpass in the next 10 - 20 years in order to realize a future plasma-based lepton collider, and 
how the work described in this dissertation is placed along that roadmap. 

1.3.1 Future Collider Design Requirements 

We initially need to briefly summarize the requirements that a future lepton (electron-positron) 
collider must satisfy in order to produce interesting new physics [3, 175], since a plasma 
accelerator must be designed according to them [178]. First, the particle energy needs to be at 
least 250GeV\ For comparison the previous largest electron-positron collider, CERN's Large 
Electron-Positron (LEP), accelerated particles up to 105 GeV. 

Equally important is the luminosity L , or the number of collisions per cross-sectional 
area, defined as 

L = f^ (4) 

A 

Here / is the repetition rate of the collisions (per turn for circular colliders), N is the number 
of particles per collision in each beam, and A is the area of each beam. The Tevatron at 



^ For a total of 500 GeV per collision. 
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Fermilab has enough energy to create the Higgs Boson, but its luminosity is too low to provide a 
satisfying measurement, on the order of 10^^ cm"^s"\ This is why LHC, CLIC and ILC luminosities 
are targeted towards the 10^^ cm'^s"^ value. Furthermore, in order to produce high-quality data 
the beam quality also needs to be very high, which translates to 0.1 - 10% energy spread and 
emmitance of 0.01 - 0.1 mm-mrad. 

In order for this design to be constructible at a fraction of the present design cost 
estimates for CLIC and ILC, currently at $10 billion per TeV, an average accelerating gradient 
larger than 100 MV/m is required along with a wall-plug efficiency (the fraction of electrical 
power that translates to beam power) of 10 - 20%. In comparison, SLAC operated at 40 MW and 
CLIC & ILC will require a constant supply of 200 - 300 MW of AC power to operate at those 
efficiencies. 

1.3.2 Roadmap to a Plasma-Based Collider 

Presently the most promising technology for a plasma accelerator appears to be an electron- 
beam driven PWFA as opposed to a laser-driven LWFA, for two primary reasons. First, lasers are 
limited by diffraction and it is not straightforward to extend the laser-plasma interaction length 
from 1 cm (which is the case now [106]) to several meters, required to achieve multi-GeV 
particle energies. Second, despite today's PW-scale power, lasers do not provide enough energy: 
a required 250 GeV beam with 6x10^ electrons (1 nC charge) contains 250 Joules of energy, as 
opposed to 1 Joule of energy that most present-day LWFA lasers can offer at a very low 
repetition rate. In contrast, the 50 GeV SLAC electron beam provides about 100 Joules of very- 
low-entropy energy hundreds of times every second. 
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PWFAs have already demonstrated that energy doubling at 50 GV/m gradients over a 
meter-long plasma is achievable. With that starting point, the following major milestones will 
need to be reached in the next few years in order to realize a future plasma-based collider: 

1. High-energy acceleration of a trailing electron beam with high charge and low 
energy spread 

2. High-energy acceleration of a trailing positron beam 

3. Wake excitation from a high-efficiency driver 

4. Multiply the energy gain either by staging or by increasing the transformer ratio^ 

All 4 of the above milestones may depend on the use of multiple bunches, which is the main 
focus of this dissertation. 

Regarding item 1, the accelerated electrons observed in PWFA experiments to date have 
originated either from the same bunch that drove the wakefield or from unshaped witness 
bunches that were much longer than a plasma wavelength, thus sampling a large phase of the 
wake which increased the energy spread and spoiled the beam emmitance. Eventually a witness 
bunch that is much shorter than the plasma wavelength will have to be placed appropriately and 
accelerated with high-gradients. Such an experiment was performed at the Accelerator Test 
Facility (ATF) of the Brookhaven National Lab (BNL), and the main results are presented in 
Chapter 4 here (also published in [85]). Specifically, it is shown that a trailing bunch which is 
shorter than the plasma wavelength is accelerated under 150 MV/m loaded gradients for the 
first time. This result, combined with SLAC-type propagation in a meter-scale plasma experiment 
could lead to the realization of the first milestone. The process would then have to be repeated 
for an injected positron witness bunch (item 2). 



^ For a definition of the transformer ratio, see section 2.3. 
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Items 3 and 4 deal with high efficiency and multiplication of the energy in a single-bunch 
PWFA stage. Coupling 250 Joules of energy (required for a strong 250 GeV witness bunch as we 
showed earlier) with a realistic 10% efficiency using plasma acceleration requires 2,500 Joules of 
input beam energy. This is one order of magnitude more energy than a single electron bunch 
can carry today. As a result, multiple bunches will most likely need to be used to provide the 
necessary energy, either through staging (by placing many single-bunch PWFA cells in series), or 
by using multiple bunches in a single stage [118]. As we will analyze in Chapter 5, the latter 
solution not only could provide the necessary energy but also optimize the beam-to-plasma 
energy transfer efficiency without the need for specially shaped bunches (only their positions 
and charge need to be tuned). 

Furthermore, staging significantly reduces the average energy gradient of the total 
accelerator because it requires relatively long non - PWFA modules to be inserted between the 
different stages to achieve extraction of spent beams and injection of new ones. On the other 
hand, multiple bunch PWFA schemes can multiply the energy of a witness bunch in a single 
stage, thus eliminating the need for staging and preserving the high gradients. Such high 
transformer ratio scenarios will also be analyzed theoretically and through simulations in 
Chapter 5. In addition, some preliminary experimental results using multiple electrons bunches 
to drive wakefields will also be presented in Chapter 6. 

Multiple bunch schemes are favored in the linear or mildly nonlinear regime of plasma 
acceleration because the sensitive nature of the phasing of the bunches in the waves is not 
affected by the nonlinear frequency shifts that occur in the blowout regime. In addition, 
although the transverse focusing properties are not optimized in this regime and the gradients 
may not be as large, the linear regime offers design simplicity and symmetry in the acceleration 
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between electron and positron witness bunches (neglecting for now spot size constraints 
imposed by luminosity requirements). It is also less susceptible to instabilities such as the hose 
instability [44] which may prevent the beam propagation over very long distances otherwise. 

The rest of the dissertation is organized as follows. Chapter 2 analyzes results from the 
PWFA linear theory relevant to understanding the experiments and simulations presented later. 
Chapter 3 provides a series of experimental data regarding the design and diagnosis of plasma 
sources. Chapter 4 demonstrates the high-gradient witness bunch acceleration in the wake 
driven by a preceding bunch. Chapter 5 presents designs and simulations for high-transformer- 
ratio and high-efficiency multibunch PWFA drivers. Chapter 6 demonstrates experimentally 
some preliminary schemes for multibunch PWFA accelerators. Finally, Chapter 7 summarizes the 
results and points out future research directions. 

1,4 Chapter Conclusions 

The following items are concluded from this chapter: 

• Particle accelerators are the main tool for doing modern physics but they have limited 
gradients; fields greater than 150 MV/m cause cavity wall breakdown. 

• Plasma - based accelerators overcome the gradient limitations because of the collective 
motion of the free electrons and provide a promising potential alternative. 

• A future electron-positron collider will require high energy (> 250 GeV), high luminosity 
(10^^ cm'^s"^) and high beam quality (low emmitance and energy spread) achieved over 
relatively short distances. Present plasma accelerator schemes do not solve these issues 
satisfactorily. 
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This dissertation demonstrates high - gradient acceleration of a trailing bunch and 
proposes the use of multiple electron bunches to multiply the energy of short trailing 
beams at a fraction of the conventionally required distances in an efficient manner. 
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In mathematics you don't understand things. You just get used 
to them. 

- Johann von Neumann 



Chapter 2 

PWFA Linear Theory 



In this chapter the basic principles of plasma wakefield theory will be analyzed. The focus will be 
in the linear regime and specifically in those aspects fundamental to the understanding of the 
physical principles of the simulations and the experiments presented the in later chapters of this 
dissertation. Apart from the generated wakefield amplitude we will analyze the transformer 
ratio, energy efficiency, beam loading and also the transverse effects. 

In principle the linear regime is valid when the plasma density perturbation Sn is small 
compared to the neutral plasma density n^, or Sn/n^'^l. Since the density perturbation is 
typically on the order of the peak beam density n^^, this criteria is equivalent to u^qI nQ<^\. In 
practice and through simulations it is found that the predictions of the linear theory are good up 
to ^,o/^o~l [116]- 
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2 . 1 Longitudinal Wakefields 

2.1.1 Basic Equations 

A number of assumptions are made when deriving the plasma wal<efield equations in the two- 
dimensional linear regime. 1) There is azimuthal symmetry, dld(p = 0. 2) The plasma ions 
remain static since their mass M. is much larger than the electron mass m . 3) The second order 
perturbations of all the quantities are much smaller than the first order perturbations. 4) The 
plasma electrons and ions are cold and have uniform neutral density. 5) The beam velocity is 
close to the speed of light c , i.e. v^-cz. 

Next we define the electron beam density distribution as follows 

n^{z,(p,r,t) = n^{^ = Z-ct,r) = n^^^{^)'n^^{r) (5) 

Hence the beam charge density is Pj^ =-enj^ and the beam current /^ =-zecn^ . Here (^ = z-ct 
is a variable in the frame of the moving beam, and ^^i(^) is an amplitude-normalized to 1 

unitless probability distribution such that ^^^(0) = 1. Combining Maxwell's equations [121], 

Newton's second law of motion for the electrons [130] and the equation of continuity for the 
plasma, it can be found (see Appendix A) that the plasma electron density is then n^-n^+n^, 
where n^ = 5n is the density perturbation which follows the harmonic oscillator equation 



n , / 2,^ _ r 2 



de 



+ k^n,=-k^n, (6) 



and its solution as a function of ^ is given by the Green's function integral 

n, (r,|) = -k^n,, (r) | n„ {^')sin(k^ (# - ^'))d^' (7) 

-00 
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Here kj^=2n I X^^ is the plasma wavenumber. The corresponding longitudinal electric field 
perturbation E^[r,^^ can be split into a longitudinal component £^||((^)and a transverse 
component R{r) as £'^(r,^) = £j|(^)-7?(r) which are given by [92, 116, 174] 

(8) 

R{r) =^ j J^' j /2,^ (r')i^o (^;,^/^' + r'^ - ^rr'co^cp'Y'dr' 

2^^'=0 r'=0 

Notice that £'||((^) = Z'((^). For a complete derivation see Appendix A. K^ is the zeroth order 
modified Bessel function of the second kind [4]. R{r) is a unitless function that describes the 
transverse dependence of the wakes and it has the property that /?(r) - n^^ (r) for transversely 

wide beams. The equality is satisfied exactly when R[r) = n^^[r) = \. It is interesting to point 

out that in that ID limit, equation (8) shows that the longitudinal wakefield after the beam is 
proportional to the cosine Fourier transform of the beam density. 

2.1.2 Transverse Component 

Experimentally the most common transverse beam density profile is a Gaussian, defined 

through n^^ (r) = e'"^ ^""^ . We shall examine these profiles closer since they are most commonly 

assumed throughout this dissertation. Using equation (8), the transverse component Rij) of 
the longitudinal wake can be evaluated numerically and is plotted in Figure 2.1 for different 
beam spot sizes. The on-axis value /?(r = 0) as a function of cr^ is plotted also in Figure 2.1 for 
different asymptotic expressions and it is given by the formula [116] 
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Transverse Wake for Gaussian Bunches 
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Figure 2.1: Left: Transverse dependence R{r) of the longitudinal wakefield for transversely Gaussian 
bunches. Right: The on-axis transverse connponent of the longitudinal wakefield 7?(0) . The two 
asynnptotic expressions for extrenne values are also plotted here. 
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^X[0-05797-ln(^^^J] 
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V J 
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k^CT^«\ 



yk^a^ 



^^^,»1 



(9) 



^e-" 



Here E^{h)= I dt is the exponential integral of the first kind. We observe that for beams 

1 ^ 



narrower than ^^cr^-lthe transverse effects become important, specifying the transition 
between ID and 2D regimes. 

2,2 Transverse Wakefields 

In this section we will investigate the transverse wakefield of a beam in the linear regime, which 
for the relativistic beams studied here is equal to Wjj,^^ = F^I e^E^-cB^, and determines 
the transverse forces on the beam [33, 129, 167]. The transverse wakefield can be estimated 
once the longitudinal wakefield W, -E^ is known by using the Panofsky-Wenzel theorem [139] 
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or simply by taking the transverse component of Faraday's law: 

— L = ^ ^ ^^ = \^d^ = Z(^)R'(r) . Equivalently [30, 33, 166] 

dr d(^ ^ dr 

W,{r,^) = ^]n,,{^')sm{k^{^-^'))d^'.^ (10) 

^0% -00 ^^ 



2.2.1 Focusing 

The envelope equation of a transversely Gaussian beam in one dimension (neglecting space 

charge effects which are of the order of 1/^ ) [35, 72, 161] is 



d a^ ^^ s 



^,- + Ka^=^j (11) 

Here C^ is the transverse spot size of the beam along the dimension x, K is the focusing 
strength, and ^^ is the normalized emmitance of the beam. Deviations in the initial spot size of 

the beam are minimized when the beam emmitance is matched such that — ^^ = 0/ oi" 

dx 

s^=y4KcjI (12) 

Inside an ion channel (void of electrons) there is no azimuthal magnetic field [B^=0) and the 
remaining transverse electric field E^ (calculated from Gauss's law) increases linearly with r, 
while the focusing strength is a constant: 

ymc zs^ymc ly 
The matched emmitance is then given by the formula 



^n = y-k,<yl (14) 
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This scenario, which often occurs in the blowout regime of PWFAs, is highly desirable since 
either the beam spot size remains constant (if the emmitance is matched) or it stably oscillates 
around the equilibrium value. 

In other cases though, the focusing strength is not constant and it may depend on the 
longitudinal position along the bunch. It can be proven from first principles that it is then given 
by the more generalized equation [72] 

K{^.) = ^^i^-W,{x)) (15) 

Here W^ is the transverse wakefield acting on the beam in that dimension and 

+00 

x-W^{x)= j x-W^{x)-n^{x)dx (16) 

-00 

In the above expression ^^(-^) is the ID density probability distribution of the beam in that 
dimension^ 

2.2.2 Emmitance Balancing in the Linear Regime 

Specializing now in the linear regime of PWFAs, by replacing the transverse wakefield from 
equation (10) on equation (15) we have 

K{^;a^)= ^ Y 2 \xR'{x)-n^{x)dx (17) 

rmc (7^ _i 

The approximately "matched" or balanced emmitance is then found from equation (12) as 

{^■,a^) = 4^^%{^)-^. (18) 

Y mc 



^N 



^ For a Gaussian beam density ^^(x) = l/v2^o-_^ -expf-x^ /2o-M the constant value of K = k^^l2y is 
retrieved for an ion channel, independent of <y^ . 



22 



X 



U.J 

0.4 
0.3 
0.2 
0.1 

A 




1 1 i 1 
III 


4 




! ! ! 1 


















^ 




















1/ 












^ 


1 








III 



10 



-1 



vO 



10 



10 
ka 

P X 

Figure 2.2: Transverse focusing dependence for a Gaussian distribution, useful for calculating the 
ennnnitance that balances the plasnna's transverse focusing force. 



The function l{o'x) defined in equation (17) can be evaluated numerically for a Gaussian beam 

_^ 
distribution n^(^x) = —^ — e ^"^^ . The minimum value of the emmitance that can prevent a 

single beam from collapsing under the focusing force is found by assuming a transverse 
Wakefield on the order of Z-E^lk^, where E^ is the maximum value of the longitudinal 



Wakefield £',, inside the beam. Then we find that 



^ mccOp 



(19) 
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2 ■ 3 Transformer Ratio 

2.3.1 Physical Interpretation 

The transformer ratio R is defined as the ratio of the maximum accelerating wal<efield behind a 

bunch ^max/ divided by the maximum decelerating wakefield inside the bunch E^^^. It is an 
artificial measure on how efficiently energy is transported from a drive bunch into a witness 
bunch through the plasma oscillations. 

R^^^ (20) 

E~ 

max 

The importance of the transformer ratio can be illustrated in an ideal case as follows. 
Assume that a single point charge Q drives a wake, its total (kinetic) energy being Wq . The 

charge will generate a maximum decelerating wakefield E^^^ and it will lose all its energy after 

traveling a distance L, such that Wq=E^^L. Now suppose that another test charge trails 
behind and gains energy by sampling the wake of the first charge. It can only gain energy over 
the same distance L , because in principle after that distance the first bunch has deposited all 
its energy to the plasma waves and has come to a halt. The total energy AW gained from the 
second charge is 

W 

^W = EI^L = E:^-^ = RW, (21) 

E^ 



max 



Therefore the somewhat surprising result is reached that the total energy gained does not 
depend on the wakefield amplitude, but only on the transformer ratio. 
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2.3.2 Achieving Large Transformer Ratios 

In most plasma acceleration experiments to date the transformer ratio has not been an 
important factor because the interaction between the plasma and the beam happens over 
relatively small distances and the drive beams do not have enough time to deposit a significant 
fraction of their energy. However a realistic future plasma accelerator will have to achieve high 
transformer ratios in order to effectively multiply the energy of an incoming beam over short 
distances. 

According to the fundamental wakefield theorem, the transformer ratio of a single 
symmetric bunch in a linear homogeneous lossless plasma has a maximum value of 2 [10, 152, 
153] (for a simple proof, see Appendix B.l). In order to achieve higher transformer ratios, either 
an asymmetric bunch must be used or more than one driving bunches must be fed into the 
plasma. In the former case, it was found by Bane et al. [9, 32, 104] that a single bunch which has 
a linearly ramped longitudinal density profile (along with a sharp kick in the beginning^) can 
create a constantly retarding wakefield along the bunch length, with the efficiency approaching 
100% and a transformer ratio theoretically up to 2nM , where M is the length of the bunch in 
units of kp (see details in Appendix B.2). However these shaped bunches are not trivial to 

generate using present-day technology [49] and no PWFA experiment has been attempted using 
a single ramped bunch. 

Another promising way to achieve high transformer ratios using regularly shaped beams 
is the use of more than one drive bunches, as initially hinted by Laziev et al. in 1988 [104, 154]. 
Power et al. [144] and Jing et al. [78, 79] more recently have proposed schemes for dielectric 
wakefield accelerators where if the distance between the bunches and the charge in each bunch 



^ The sharp kick brings the wakefield in the optimal value, after which the ramp preserves it. 
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is adjusted appropriately, then the decelerating wakefield experienced by each bunch could be 
identical for all bunches and hence the transformer ratio and the efficiency will be enhanced. 
Those ideas can be extended to plasma accelerators^ and they will be analyzed in detail in 
Chapters. 

2,4 Energy Transfer 

In this section we will consider the driving energy efficiency of a plasma accelerator in the linear 
regime. This is of great interest because the plasma accelerator is essentially an energy 
transformer, extracting the energy from the drive beam(s) and delivering it through the plasma 
waves to a trailing beam. The efficiency is a very important parameter in the design of the 
accelerator as design efficiencies around 20% are highly desirable in order to reduce the cost of 
the accelerator. 

2 .4. 1 Driving Efficiency 

Given an electron beam with a specific initial energy, the driving efficiency is defined as the 
fraction of the total beam energy that has been transferred to the plasma waves at the end of 
the beam-plasma interaction. 

Let us focus on the energy transport between particles and fields using general 
electromagnetic theory. Assuming a current / created by a number of particles that are 
present in an electric field E , then the instantaneous power density transferred from the 
particles to the fields and vice versa is defined by J E (with units of W/m^) . For the case of a 
charged bunch distribution, the total instantaneous power P dissipated at the plasma (in 



^They are applicable to other wakefield accelerator schemes as well, as long as they share the same 
formulation. 
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Watts) is just the integral of the above power density integrated over the whole volume of the 
bunch: 



+00 Ik +00 



P =\\ \ J-EdV = 

+00 271 +00 



r=0 ^=0 ^=-00 

+00 



(22) 

^=-00 r=0 

In order calculate the total energy transferred from the beam to the plasma, equation (22) 

needs to be multiplied by the total time of the interaction, T = L/c (L is the length of the 

interaction). This assumes that the velocities of the particles in the bunch are identical and do 

not change with time; in other words, the bunch is rigid\ 

The energy transfer efficiency of the interaction is now defined as the fraction of the 

total energy of the bunch that has transferred to the plasma after time T . Assuming that there 

are N particles in the bunch, each having the same initial energy W^ ^^Vq =^£'^^L, the 

efficiency r/ is written as 

energy dissipated to the plasma _ T P _ 
total energy in the bunch N • W^ 

(23) 



ri ^■ 



\ n,,{r)R{r)2nrdr\ n,,{^)E,{^)d^ 

r=0 ^00 

NE 

max 

Note that f"^^ =£"0 •/?(0) (assuming the wakefield peaks on-axis). By inspecting the above 
expression we observe that the efficiency is maximized to unity when the electric field is 



^ Otherwise the time-dependent integral \Pdt needs to be evaluated. 
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constant for every point along the bunch distribution, E^^(^(^^ = Eq, with 7?(r) = l. Physically 

this means that every particle inside the bunch is depositing its energy to the plasma at the 
same rate. In that case the efficiency becomes 



Eo j n,^{r)2nrdr\ n,^^{^)d^ 

ri = ^^ --^ = ^^ = 1 (24) 

NE, NE, 

It may not be obvious right away that the constant electric field E^ along the bunch (in the 
enumerator of the above expression) and the electric field that corresponds to the initial energy 
of the bunch E^^^ (in the denominator of the above expression) should be the same. However 
the physical meaning of E^^^ is the electric field required such that a charge with energy 
Wq = cVq = eE^^L will come to a stop after distance L . Vq then is simply the voltage difference 
induced between the start and end points. Of course the symmetrical argument can also be 
applied: It is the electric field required to accelerate a charge e over distance L to an energy 
W =eV =eE L 

2.4.2 Transverse Efficiency 

The efficiency specified in equation (23) will be less than unity for bunches for which /?(r)<l 

for some r. For transversely Gaussian bunches, ^tii^)-^"^ ^""^ ' ^^^^ transverse driving 
efficiency which is defined as 

+00 

n.-^ (25) 

I n,^{r)R{0)rdr 
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Figure 2.3: Transverse efficiency of a Gaussian bunch (beann width ^Ink g^) and of a flat-top bunch 



(beann width k^a ). 



is plotted in Figure 2.3 as a function of the bunch width k^G^ . 7?(0) is the maximum value of 

/?(r), which here occurs on-axis. 

The transverse efficiency increases with decreasing cr^ because the wakefield variation across 
the beam is less for narrower beams. The efficiency is compared with the transverse efficiency 
of a transversely flat-top bunch (with an infinitely sharp cutoff at r = k^a ). We observe that for 

narrow bunches [k^G^ <1 or k^aKl) the efficiency is independent of the exact bunch shape 
since the wakefield always extends out to a distance 1/^^. However for wide bunches the 
efficiency of a Gaussian bunch is 50%, unlike the wide sharp flat-top bunch which approaches 
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100%. This can be explained by invol<ing physical arguments when realizing that for wide beams 
(kp(7^ »1), the transverse component of the wakefield R{i") and the transverse beam profile 

^bi{^) w'" overlap very closely since the expelled plasma electrons stay close to the beam ( 

/l^<^(Tj. Hence for wide beams R[r)^nj^^[r) = e~' ^^""^ with /?(0)^1 and in this case 

-j +00 ^ 

equation (25) will read rj^- \ [e '^""^ ) iTrrdr = 1/2 . Wide flat-top bunches consist a 



r r=0 



special case where nbi{^)^{^)-^bi{^)-^bi{^) ^^^ therefore yield transverse efficiency of 
7/^ = 1 (equivalent to setting 7?(r) = 1 which is the ID limit). 

2.4.3 Efficiency vs Transformer Ratio 

The concepts of driving efficiency and transformer ratio are very closely related. In this section 
we provide an example to illustrate the similarities and differences between these two 
parameters. Assume a single longitudinally square bunch centered at ^^ with a longitudinal 
width W containing N^ particles that has a beam density profile 



^b{^^^) 



N. 



Inaw 



-reel 



^#-#,^ 



w 



(26) 



Here the function rect[^lw) is 1 for |^|<w/2, V^ for |(^| = w/2 and otherwise. Using 
equation (8), the electric field generated by this bunch in the ID limit is 



^0% 



+ 



V W y 



^0% 



^kw^ 



V 2 , 



cos[it^(^-^o)]-C/(^-^o-H^/2) 



(27) 
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Figure 2.4: Transformer ratio and driving efficiency of a single square bunch as a function of the 
nornnalized bunch width k w . The results are connpared with the transfornner ratio and efficiency of a 

Gaussian bunch with nornnalized width 42nk a . The transfornner ratio for the Gaussian is defined at lo. 



Here U is the Heaviside step function. The driving efficiency is calculated from equation (23) 
and is equal to 

l-cos(/:^w) 



ri(k^w) = - 



k^w 



(28) 



The transformer ratio is simply from equations (27) and (20) 



R{k^w): 



{k^w/2) 



cos 



2 sin 



y 2 , 



kw< — 

' 2 



n 
kw> — 

' 2 



(29) 



The values of the longitudinal efficiency 7/ and the transformer ratio R are plotted in 
Figure 2.4 as a function of the normalized width of the bunch k^w . For comparison, the 



31 



corresponding parameters for a Gaussian bunch of the same number of particles (achieved 
when v2;r(T^ = w ) are also plotted in the same graph\ 

It is observed that although the two parameters (77 and R ) are optimized around the 
same rough range of values for k^w , they peak at slightly different points. In order to analyze 
this difference, we can compare the situations where the two above quantities are maximized, 
namely kpW = 7r (which maximizes the transformer ratio at R = 2 and then the efficiency is 

r/ = 2/ 71 = 63.6%), and kjyW = 0J427r (which maximizes the efficiency at 77 = 72.5% and the 
transformer ratio then is R = l.S4). 

In the maximum transformer ratio case, kpW = 7r , the maximum wake inside the bunch 

's ^max ^ 3nd the maximum accelerating wake after the bunch is exactly 

^max ^ — ^-^—2 = 2E^^ . Assuming that each particle in the bunch has an initial energy Wq , then 

S.TT 

the interaction length is limited by the length over which the fastest decelerating particle loses 
all its energy: 

L = --^ = ^^ (30) 

The energy left behind (transferred to the plasma) is proportional to the square of the 
accelerating field amplitude and the length of the interaction (C is a constant): 

W,^a.«. =ciE:J-L = C- 4^^W„ (31) 



^ Strictly, the transformer ratio of a Gaussian bunch is always equal to 1 (since the bunch extends to 
infinity). Here we evaluate the transformer ratio at a distance la from the center. On the other hand the 
definition of the efficiency incorporates the bunch profile and does not suffer from this technicality. 
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In comparison, in the maximum efficiency case, k w = 0.1427r , the maximum wake inside the 



€H W 

bunch is £"'' = -135 E~^^ and the maximum accelerating wake after the bunch is 



en^^w ^ . fO .14271 
-2sin - 



-1.24£''^ . So in this case the decelerating wake is 35% higher. 



"^'^ ^o0.742;r ^ 2 
and the accelerating wake is 24% larger after the bunch. The interaction length now is shorter: 

W 
V = — '— = 0J4L (32) 

eE'- 

max 

Finally, the total energy transferred to the plasma is 

Kas.a=CiEZ^f-L' = LU.W^,^^^ (33) 

So 14% more energy will be transferred to the plasma in the maximum efficiency case, even 
though the transformer ratio is actually smaller. The lower R is counterbalanced by the higher 
wake, even if the interaction length is shorter. This is because the energy transferred is 
proportional to the square of the wakefield, but only linearly proportional to the interaction 
length (which is inversely proportional to the decelerating wakefield). 

2.5 Beam Loading 

In this section we will investigate the process of energy transfer from the plasma wave to a 
witness bunch. In order for the bunch to gain energy it needs to be phased properly at the 
accelerating crest of the plasma wakefield. The witness bunch will also suffer energy spread if 
the longitudinal wakefield is not uniform along its longitudinal and transverse dimensions. 
Finally, an energy transfer efficiency from the plasma to the bunch electrons is associated with 
this beam loading. The question is then what is the optimal set of parameters for the witness 
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bunch (charge, width, phasing) to optimize the acceleration process and produce an accelerating 
bunch with low energy spread and high efficiency. 

2.5.1 Longitudinal Beam Loading 

Assume a external longitudinal plasma wakefield that has the form (on-axis) 

EU^) = EoCOs(k^^)R,{0) (34) 

Here ^o(O) '^ given by equation (9). At this point it is not relevant how this wakefield was 
created; it may have been generated by a single bunch or by a series of bunches, adding in or 
out of phase. After the drivers though, the wake that is left behind can be fully described with 
the above expression. 

For simplicity let us assume a square witness bunch with A^^ particles and length W 

phased at a position (p^ =k^(^^ after the location of the external wakefield maximum <^ = 0, of 
the form 



A^., 



iTra^w 



-red 



'^-L 



w 



2< 



(35) 



The total wakefield inside and after the bunch (on-axis) can be found by superimposing the 
external field from equation (34) to the field generated by the witness bunch from equation (27) 
. The beam loading efficiency then is found to be 
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(36) 
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Figure 2.5: Longitudinal beam loading efficiency (left) and fractional longitudinal energy spread (right) of a 
square bunch accelerated by an external wakefield as a function of the normalized bunch length k w , 

the phase (p^ of the bunch in the wake, and the number of particles in the bunch (related to the ratio 

Here sine (x) = sin (x)/x is the unnormalized sine funetion and £"0^ = ^ = — ^^— is the 

longitudinal wakefield amplitude per normalized unit length generated by the witness buneh. 
The maximum effieieney is aehieved when the buneh is phased at the most flat portion of the 
external wave (so that different partieles experienee as similar as possible aeeelerating gradient), 
i.e. (p^=kp^^=n . In that ease and for very short bunehes, /:^w^0, then in the ID limit ( 

/?Q (0) = /?^ (0) = 1 ) the expression derived by Katsouleas et al. in 1987 [92, 174] is retrieved: 



2 — —^ 

E , 

^0 J 



(37) 



On the other hand, inereased effieieney oeeurs at the expense of larger energy spread for this 

s E 

unshaped buneh. 100% effieieney happens when Eq^=Eq (or N^=-^—^27r(j^ partieles are 

e 

loaded into the wake) and means that no wakefield is left behind the buneh, therefore the 
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Wakefield along the bunch will vary from E^ (in the head) to (in the tail) resulting in 100% 
energy spread. The fractional energy spread /^ for the square bunch is given by the formula 



fw = = Pq sine 



^k^w^ 



V 2 , 



(38) 



In Figure 2.5 the beam loading efficiency and the energy spread for a square witness bunch are 
plotted in the ID limit as a function of the bunch parameters: the bunch width k^w , its phase in 
the external wake cp^, and its self-wakefield £"0^ (which is proportional to the number of 
particles in the bunch) normalized to the external wakefield amplitude E^. We observe that 
although the efficiency is substantially lower when the witness bunch is placed before the peak 
of the wake, its charge can be adjusted such that the energy spread is minimized (at some 
expense of the accelerating gradient). 

2.5.2 Reverse Ramped Bunch 

The tradeoff between efficiency and energy spread can be surpassed if a longitudinally specially 
shaped bunch is utilized, as described in [89, 92, 174]. The efficiency is maximized with no 
longitudinal energy spread if all the particles on-axis feel the same accelerating wakefield. Using 
equation (8) and assuming an external wakefield of the form in equation (34), by requiring the 
total field inside a bunch with its head phased k^(^^ relative to the external wake to be constant, 
the optimal longitudinal beam distribution is found 

while the phase is adjusted between 7r/2<k^(^^ <n . The (constant) accelerating gradient in 
this case is E"^ =-£'q/?q(0)cos(/:^4o) ' ^^^ maximum length of the bunch is 
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Properties of triangular witness bunch 




120 150 

Relative Phase [degrees] 

Figure 2.6: The efficiency, accelerating field and the relative nunnber of particles accelerated for a rannped 
witness bunch in a sinusoidal wakefield as a function of its position relative to the external wakefield. The 
longitudinal energy spread is zero for this specially shaped bunch. 



kpL = -tan(kp(^], the longitudinal efficiency is 7/^ =sin^(fc^4w) ^^^ ^^^ ^o^^' number of 

sin' (/:/o) s^Ittct^E^R^ (O) 
particles in the bunch is N^=-Nq -. -r , with A^q ~ ^^^ ^ 

2cos(/:/o) e 

transversely Gaussian profile. Although now there is no energy spread, the efficiency (and the 
number of particles) are increased at the expense of the accelerating gradient. Those tradeoffs 
are depicted in Figure 2.6. 

2.5.3 Transverse Beam Loading 

So far we have ignored transverse effects in the beam loading process. Now suppose that the 
same square witness bunch defined in equation (35) is found inside a 2D external axi-symmetric 
wakefield 
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EA^^r) = E,cos(k^^)R,{r) (40) 

The Wakefield generated from that witness bunch (after the bunch) has the general form (see 
also equation (27)) 

E„{^,r) = E^cos(k^^ + g,„)R„{r) (41) 

Here (p^ is the phase shift of the witness wake relative to the wake of the external wakefield. 
The total energy in the wakefield before the witness bunch over one period of the wave is (for 
wide beams where the other field components can be neglected) 



27ilk„ 



2 +^ 

W,,, = I Inrdr \ d^El cos^ (k^4)R^ (r) = ^^o J rdrRl (r) (42) 

The total wakefield after the bunch is given by the superposition of (40) and (41): 

= ^0^0 ( n/i + P' + 2P cos (p^ sin (^/ + (p^^,^^ ) (43) 

'^^ EMr) 

After calculating the energy W^^^^^ in the wake left behind by the witness bunch similar to 
equation (42), the beam-loading efficiency of the system can be written as 

+00 

l{\ + p^+2pcos(p^,)R^{r)rdr 

%=i--f^ = i-^^ r. (44) 

lR^{r)rdr 



Some interesting special cases can be discussed here. First, if E"^ =0 then % =0 (no energy is 
transferred from the wake to the bunch). Second, if the bunch that excited the wake and the 
witness bunch have the same radial dependence (i.e., /?Q(r) = /?^(r)), then the transverse 
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components of the wakes overlap perfectly in space (albeit with different strengths) and the ID 
result of equation (36) is retrieved, yielding the efficiency 



r . ^ 



^0 



^0 J 



(45) 



Negative beam loading efficiency means that the bunch transfers energy to the plasma instead 
of the other way. In order for the witness bunch to absorb energy, the phase of the witness 
bunch must be selected such that cos^^ <0. The maximum beam loading efficiency is obtained 
when (p^ =71 , in which case the external wave and the bunch wake are exactly out of phase and 
cancel out. 

To illustrate the principles of beam loading, we can plot in Figure 2.7 the efficiency 
according to equation (44) when the external wake is of the form (40) with amplitude E^ and 

transverse dependence ^o(^) corresponding to a transversely Gaussian bunch with spot size 
(Jq. The witness bunch is placed at the accelerating phase, (p^=n . The efficiency is plotted 
versus the Wakefield amplitude ratio, E^l E^,^n\\^xq. E^ is the amplitude of the sinusoidal wake 
that the witness bunch alone leaves behind it (proportional to the number of particles in the 
witness bunch). We also vary in the same plot the transverse Gaussian spot size of the witness 
bunch, G^, by plotting curves for different a^ /cTq, for three different values of ^^cTq- Note that 
the cr^ / cTq = 1 curves are the same as the (p^^n curve in Figure 2.5. 

We observe that when the drive and witness bunches have the same transverse width ( 
ct^/cTq =1) then the beam loading is 100% when E^=E^. This of course happens with 100% 
energy spread since the two waves overlap out of phase everywhere in space. For narrower 
witness bunches, g^I g^<\, the efficiency is reduced and denser bunches (higher E^ ) need to 
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be utilized to compensate for the narrower witness wakes. After the efficiency peal<s for each 
case it eventually becomes negative since the wake energy left behind is mostly that of the 

witness bunch. Note that simply by changing a^ will also alter the witness wake E"^ by l/o"^ - 

see equation (35). For example, for the case of k^^G^ = 2.4 (top panel), if initially cr^ / cTq = 1 and 

E^/ Eq =0.1 (20% beam loading) then by making the spot size 10 times smaller (a^^ I a^ =0.1) 

will also cause E^ I E^ =10, yielding an efficiency of ^28%. 
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Figure 2.7: Beam loading for a transversely Gaussian witness bunch of spot size a^ in the presence of an 
external wave of amplitude E^ created by transversely Gaussian bunches of spot size cr^ . The wake 
amplitude left behind by the witness bunch alone is £ . The different colored curves correspond to 
narrower witness bunches. The curves for 3 different values of k^a^ are plotted in this figure. 
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It is advantageous to use narrower drive bunches (bottom panels) since the wake 
energy is more efficiently absorbed. The reason for this is that when driving with a narrow 
bunch then the beam loading is insensitive to the width of the witness bunch since the wake 
always extends out to a plasma skin depth. Furthermore, narrow drive bunches increase the 
driving energy transfer efficiency (see Figure 2.3). The drawback is that the absolute number of 
particles N that can be used in the drive bunch has to reduce as o-^ decreases, since the drive 

beam density ^^q ^N / (Jq cannot be raised arbitrarily without going into the nonlinear regime 
where these derivations break down (for a fixed plasma density). The plasma density could be 
decreased to allow for more particles, but only at the expense of the wakefield amplitude which 

scales as £ oc .Jn^ for fixed n^^ / n^ (see for example equation (27) for the scaling). 

2.5.4 Transverse Energy Spread 

As with the longitudinal beam loading, there is always a tradeoff between the transverse beam 
loading and transverse energy spread. This energy spread is in addition to whatever energy 
spread may be imparted on the bunch by the longitudinal effects and field variation along the 
bunch. The transverse energy spread is caused by the off-axis field variation along the radial 
profile of the bunch. Unlike the transverse beam loading efficiency, it is not independent from 
the longitudinal profile of the bunch. Therefore, we shall examine the transverse energy spread 
on a longitudinally square and transversely Gaussian bunch of the form specified in 
equation (35), which is loaded onto an external wakefield of the form 

£^^(<^,r) = £'QCOs(/:_^<^j/?Q(r). Since the wake left behind by the witness bunch alone is equal 

to E^(^(^,r) = E^cosykp(^ -\-(p^jR^(^r) , the total wakefield inside the witness bunch is a 
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superposition of £'^_^^(^,r) and the wal<efield found by substituting the bunch profile into 
equation (8), and can be written as 



1 E R (r) (46) 



2 sin 



V 2 J 



EoM^) 



The fractional transverse energy spread (defined at tlie radial tail r = 3(7^ of the bunch^) is found 

E ., (f,r = 3cr ) 

^ _ 1 inside \ ^ ' wj_ lAl) 

The above energy spread is plotted in Figure 2.8 as a function of the various parameters of the 
system for (p^=7r. 

Two regions are clearly distinguishable in these graphs. For light beam loading (left 
side), the energy spread is determined by the external wakefield and E.^^.^^ ~/?Q(r) only. For 
heavy beam loading (right side), the energy spread is determined by the transverse variation of 
the bunch wakefield along the bunch, i.e. E.^side ^ ^w{^)- For transversely similar bunches, 
a^=(jQ, the energy spread is independent of the witness bunch charge since the fields overlap 
everywhere in space. We observe that in order to retain a small energy spread while loading a 
significant amount of charge, a witness bunch narrower than the drive wave must be used. 



^As with the transformer ratio, the energy spread for non-finite bunches strictly is always 100%. The 
emmitance (area of the beam particles in phase space) is a much better definition of beam quality, 
although it complicates the issue by including spatial particle information as well. 
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Figure 2.8: Fractional transverse energy spread for a longitudinally square (of width k^w ) and 
transversely Gaussian witness bunch (of spot size g^ ) in the presence of an external wave of annplitude 
E^ created by transversely Gaussian bunches of spot size a^. The wake annplitude left behind by the 
witness bunch alone is £'^^y^^^ = E^ . The differently colored curves correspond to narrower witness 
bunches. The dashed curves correspond to bunch width k^w^nl \ , while the solid curves correspond to 
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2,6 Chapter Conclusions 

The following items are concluded from this chapter: 

• The transformer ratio can be greater than 2 in the linear regime if either non-symmetric 
or multiple bunches are used. 

• The electron bunches suffer from nonlinear focusing forces in the linear regime, which 
can be approximately tuned out by adjusting the emmitance of the individual bunches. 

• The beam to plasma longitudinal energy transfer efficiency is maximized (for single 
bunches) when all particles lose energy at the same rate, i.e. the field is constant under 
the beam. The transverse efficiency is maximized for narrow beams. Regularly shaped 
single bunches have maximum efficiencies around 60 - 70%. 

• Efficiency and transformer ratio are maximized almost simultaneously, however high 
efficiency schemes transfer more energy to the plasma over shorter distances even if 
the transformer ratio is smaller. 

• There are certain tradeoffs between efficiency, energy spread and accelerating gradient 
when the wake is loaded with a trailing beam. The resulting energy spread can be 
reduced either at the expense of the accelerated charge (unless the beam is specially 
shaped), or at the expense of the gradient (by placing the beam earlier than the peak 
accelerating wakefield). 
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Science is facts; just as houses are made of stones, so is science 
made of facts; but a pile of stones are not a house and a 
collection of facts is not necessarily science. 

- Henri Poincare 



Chapter 3 
Plasma Sources 

3-1 Overview 

In order for the plasma accelerators to operate optimally, reliable plasma sources are required. 
These sources must be able to produce high density (10^^ - 10^^ cm'^), low temperature 
(<10eV^) and longitudinally uniform plasmas over relatively long distances (> 1cm) in a 
controllable, repeatable and consistent manner. 

3.1.1 Gas Jets and Metal Vapors 

In some occasions the driver (either laser or electron beam) is intense enough to ionize a neutral 
gas and create the plasma, inside which the wakefields are excited [45]. Such was the case in the 
recent LWFA experiments by Malka et al., Faure et al., Geddes et al., and Mangles et al. [58, 59, 
63, 90, 119], where they managed to systematically produce relatively monoenergetically 
accelerated electron bunches of good quality (low emmitance and energy spread) using gas jets. 



^ What is usually quoted when we refer to a plasma temperature T is, by convention, the corresponding 
energy k^T . 
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In these setups there is typically a tank that holds a gas such as helium or hydrogen in 
very high pressure. This gas is released into vacuum through a valve and a small nozzle. The 
front of the focused laser pulse ionizes the gas, while the rest of the pulse expels the electrons 
and initiates the plasma oscillations. For instance, Faure et al. used a helium gas jet that had a 
nozzle diameter of 2 mm. The 10^^ W/cm^ titanium-doped sapphire laser was focused 800 lum 
away from the nozzle, generating a longitudinally flat plasma density on the order of 10^^ - 
10^^ cm"^ over a 1.5 mm length with 400 lum density ramps on each side. However in those gas 
jet systems the acceleration length is limited by the length of the longitudinal plasma density 
plateau which in turn is limited either by the small nozzle sizes, by the length that can be ionized 
by the laser or by the dephasing length, thus it cannot be straightforwardly extended to multi- 
cm-scale lengths. 

In a similar fashion, a lithium vapor inside a heat-pipe oven [171] was used as a plasma 
source in the SLAC energy doubling PWFA experiment [20]. The SLAC 42 GeV electron bunch is 
50 fs long, carries 1.8x10^° electrons and can be focused down to a 10 lum spot size, thus being 
intense enough (electric field > 5 GV/m) to ionize the lithium gas\ The meter-long plasma 
density (on the order of 10^^ cm'^) was adjusted with better than 10% accuracy by tuning the 
neutral gas pressure and assuming 100% ionization [131]. 

3.1.2 Capillary Discharge Plasmas 

The plasma sources that are the focus in this dissertation are capillary discharges. As opposed to 
gas jets and metal vapors, capillary discharges do not require an external strong driver to ionize 
a neutral gas. Rather, the plasma is produced by applying a high voltage (typically in the several 
kV range) along a capillary tube, which typically have lengths anywhere from a few mm to 



^ Lithium has very low first ionization energy for a metal, equal to 5.4 eV. 
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several cm and transverse lengths up to couple of mm. The plasma electrons originate either 
from the ionization of a gas that has pre-filled the capillary tube (gas-filled capillaries [65, 95, 
159]) or from the evaporation of the capillary wall material (ablative capillaries [143, 158, 164, 
180]). Currents in the kA range may be required to break down the materials, so capillary 
discharges need to be supported by high voltage circuits that deposit enough energy to the 
capillary. 

Capillary discharges have many favorable characteristics for plasma wakefield 
acceleration. They can be used to extend the acceleration lengths, especially in LWFA 
experiments because the typically cylindrical symmetry can create a near-parabolic in the radial 
dimension profile, which favors laser channeling [52, 65, 84, 107, 159]. Capillaries are also great 
plasma sources for weak driving electron beams or lasers that are not intense enough to ionize 
the gas. In addition, the extra parameters introduced, such as the diameter or the initial gas 
pressure, allow for better control of the selected plasma density. Since the discharges may last 
for several jus before they die out and the relativistic beams travel close to the speed of light, a 
beam will transverse the capillary length in a matter of ps. The proper timing between a beam 
and the discharge also allows for a choice in the plasma density. Typically densities on the order 
of 10^"^ - 10^^ cm"^ can be obtained. 

On the other hand, the density produced at a capillary may be sensitive to parameters 
that are not easy to reproduce widely, such as the exact material composition, the shape of the 
discharging electrodes and the internal impedance of the power supply along with the current 
pulse it can produce. Since each experiment usually builds its own capillary and power supply, 
there are difficulties in the reproducibility of the results among the researchers and there seems 
to be a lack of common knowledge in the literature. In addition, due to the complexity and the 
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violent nature of the discharges there exists no solid theory that predicts the produced plasma 
densities with accuracy and only recently a few simulation codes have been published [21, 24, 
25, 109, 181] in an attempt to validate experimental results. 

The goal of this chapter is to present experimental data for capillaries that were taken at 
the Pulsed Power Lab at the University of Southern California and at the Accelerator Test Facility 
at the Brookhaven National Lab. We examine how the various parameters of the capillary 
geometry (size, length and diameter) and discharge characteristics (charging voltage, gas 
pressure and electrode shape) affect the generated plasma density. In addition, the evolution of 
the plasma density with time, which is critical for the experiments described later in this 
dissertation, is also analyzed. 

3-2 Spectral Diagnostics 

Plasmas emit light; the light emitted at different wavelengths is a signature of the plasma and 
carries information about the processes that occur inside it. By collecting that light one can work 
backwards and estimate the plasma parameters such as the density and the temperature. There 
are a number of spectroscopic techniques available for analyzing plasmas, and we will review 
both standard used methods [76] as well as methods used specifically in the experiments 
described later in the dissertation. 

3.2.1 Typical Spectroscopic Methods 

A popular method for diagnosing both the plasma density and temperature is Thompson 
Scattering, which is based on the interaction between an electromagnetic wave and charged 
particles [156]. The diagnostic setup requires an external laser source that is focused into the 
plasma. The electromagnetic field of the laser will cause the plasma electrons to start oscillating. 
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Figure 3.1: Experimental setup layout for the collection of the plasnna light for the ablative discharges. 
Source: Daniil Stolyarov, BNL 



and as the electrons quiver they emit radiation which is strong around the wavelength of the 
incident field (scattered wave). Since the interaction cross section is on the order of 10"^^ this 
radiation is much weaker than the incident radiation. As a result a powerful laser needs to be 
used (such as the Ruby at A = 0.694 lum) and the scattered radiation needs to be collected and 
imaged carefully. 

More relevant to capillary discharges, a widely implemented method for diagnosing the 
plasma density is laser interferometry [150, 158, 159]. In these setups typically a Mach - Zender 
interferometer is built with one of its arms passing through the plasma. The plasma's index of 
refraction will depend on its density and will cause a phase shift in the incoming laser pulse. This 
method has the advantage that it provides full time-dependent evolution of the density during 
the plasma discharge; its main drawback is that the laser light may be deflected or absorbed by 
the plasma and pollute the interference signal, thus the density has to be extracted from the 
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Figure 3.2: High-density time-integrated hydrogen plasnna discharge spectrunn. The absorption line around 
590 nnn is attributed to the nnolybdenunn ends of the electrodes. 



number of fringe shifts with extra care. For example, using a He-Ne laser over 3 mm of plasma a 
full phase shift of 360° will occur for every 1.2x10^^ cm"^ change in the plasma density^ Another 
advantage of this method is that when a 2D spatial detector is used to record the interferometer 
signal, then information about the transverse profile of the plasma can be also gained [47, 65]. 

3.2.2 Stark Broadening of Hydrogen Lines 

The principle diagnostic utilized in this work is Stark Broadening on hydrogen plasmas [6, 22, 26, 
53, 54, 61, 62, 64, 172, 173], which emit light largely in the optical wavelengths. The main 
physical principle is that transitions of electrons between atomic levels that would normally 
produce a spectral line almost at a single wavelength, now occur between perturbed levels and 
as a result the spectral line broadens [2]. The widening of the lines is caused by the electric 
fields^ generated by the free plasma electrons, therefore in general a larger plasma density will 



^ The phase shift over a plasnna length L with index of refraction rj \s roughly A(^-k^^^^^L?j , and the 
corresponding plasnna density n^ -(4;r^QmA^)/(^^/i^^^^^ -Lj , assunning that rj^ = -^l-nj n^ -\-nJ2n^. 

^The Stark Effect is the electric field analog of the Zeennan Effect, where a nnagnetic field is applied to 
perturb the energy levels. 
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Figure 3.3: Photograph of the hydrogen capillary discharge at USC. Source: Jessica Hao Chen. 

cause a larger widening and a bigger Stark broadening [140, 141]. A sample spectrum is shown 
in Figure 3.2, and a typical hydrogen plasma discharge in Figure 3.3. 

Three different spectral techniques were used. First, the Stark broadening of the Balmer 
Hct (656.3 nm) and Hp (486.1 nm) lines: The width of the lines yields information on the density 
and the temperature of the plasma\ Second, the line intensity ratio between those two lines 
can provide information on the temperature. Finally, the temperature can be estimated from 
the background light that is emitted from hydrogen^ [80]. 

A typical experimental setup is shown in Figure 3.1 [29, 164]. The plasma light was 
collected from the side of the capillary through a series of lenses, it was fed into an optical 
spectrometer with a resolution of approximately 1 nm, and the resulting spectrum was imaged 
onto a CCD camera. In other setups one or more optical fibers were inserted transversely 
directly into the center of the capillary body to capture the plasma light at the center of the 
capillary. 

The plasma spectrum is post-processed and the FWHM linewidths A/l^/^ of the Balmer 
lines are measured by fitting the theoretically-predicted Lorentzian shapes under the lines (most 
often Hct). The plasma electron temperature T is measured either from the ratio of the Balmer 



^ The Balmer a (red) and (3 (blue-green) lines originate from electron transitions from the third (n=3) and 
fourth (n=4) atomic shells, respectively, to the second (n=2) shell. The a line gives to the famous Orion 
Nebula its reddish color. 



^ This light originates mostly from bremsstrahlung radiation. 
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Figure 3.4: The plasma density as a function of the linewidth of the Ha line for different tennperatures. The 
plasnna density is in units of cnn' and the linewidth in units of nnn. 

line intensities [16] or from the line to background ratio under one line [66]. Then, using the 
theory developed by Griem [66, 67, 68, 69, 70, 94, 160], the plasma density is estimated with a 
30% accuracy [16] using the nonlinear equation [7, 67] 



^[cm"'] = 8.02x10^' 



A^,[A] 



(48) 



Here a^^ is a parameter that depends on the plasma density, the temperature and the spectral 
line observed. The above formula is based on a number of assumptions on the state and 
behavior of the plasma, such as local thermal equilibrium (LTE) for the plasma and validity of the 
quasistatic approximation, which assumes that the spectral broadening is dominated by the 
effect of the electrons and that the ions remain static during that time. The theory seems to be 
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reliable for plasma densities up to 10^^ cm"^ and temperatures up to 5 eV. In Figure 3.4 the 
dependence of the plasma density on the measured linewidth for different temperatures is 
plotted for the case of the Ha line. 

3-3 Dependence of the Plasma Density on the Capillary 
Parameters 

The goal of this section is to investigate the dependence of the plasma density and temperature 
on the various parameters controlling a capillary discharge. These data are important because 
they allow design flexibility when a certain plasma source is to be designed for a plasma 
Wakefield experiment. Most of the work described in this section was performed at the Pulsed 
Power Lab of the University of Southern California using transparent glass capillaries filled with 
Hydrogen gas. Unless otherwise noted, the charging voltage was set at 60 kV. The plasma light 
was collected by imaging the light onto a lens-to-fiber coupler, which was then fed through an 
optical fiber to an optical spectrometer with resolution of 1 nm. Every possible measure was 
taken so that only one parameter was varied at any given time, everything else remaining 
constant^ 

3.3.1 The Effect of Neutral Gas Pressure 

In this section we present a comparison of the spectra recorded for different neutral hydrogen 
pressures, up to 2 atm. The data are presented in Figure 3.5, where the spectra are normalized 
to the peak of the Balmer Ha line. 



^ Unfortunately the energy transferred from the power supply to the plasma cannot be independently 
controlled externally since it depends on the impedance of the plasma, which in turn depends on the 
discharge conditions. 
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Figure 3.5: Normalized spectra recorded for increasing neutral hydrogen pressures. The right peak is the 
Ha line and the left peak is the Hp line. 



It is evident from this figure that higher plasma densities (wider linewidths) and stronger 
light emission can be obtained for higher pressures. The linewidths of both lines increase, 
suggesting an increase in density. This is also accompanied by a change in the ratio of the two 
line intensities since we observe that the Hp/Ha ratio increases with pressure, suggesting an 
increase in the plasma temperature with pressure, according to the line-ratio technique that 
was mentioned previously. Quantitative measurements of the density and temperature are 
presented in the following sections. 

For low pressures and densities the line shape of the Ha line remains mostly Lorentzian, 
which strengthens the validity of the density calculations. For larger pressures the shape is not 
purely Lorentzian; however the right side of the Ha line is still Lorentzian enough and in those 
cases the twice of the right-side half width at half maximum was used to estimate the density. 



55 



Plasma Density for Various Capillary Lengths 
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Figure 3.6: The plasnna density for different capillary lengths as a function of pressure. ''New pulser'' refers 
to shots nnade after sonne slight nnodifications to the power supply unit. The charging voltage was 60 kV 
and the capillary dianneter 0.6 nnnn. 



3.3.2 The Effect of the Capillary Length 

Figure 3.6 presents a set of data that depicts the plasma density dependence for different 
capillary lengths. Each data point is an average of 3 - 5 consecutive discharge shots with the 
error bar (min - max) being on the order of 10% (not shown). The inaccuracy due to the density 
extraction from linewidth and temperature uncertainty is higher than this shot to shot variation. 
It was observed that in general the density increases with pressure and initially 
increases with the capillary length, with a maximum density of 9.5x10^^ cm'^ at a length of 12 
mm and pressure of 2 atm. For larger lengths (> 12 mm), the density drops and capillaries that 
had more than 30 mm of length did not always breakdown for fixed voltage settings due to the 
Paschen curve limitations. 
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Plasma Temperature for Various Capillary Lengths 
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Figure 3.7: The plasma temperature for various capillary lengths. ''New pulser'' refers to shots made after 
some slight modifications to the power supply unit. The charging voltage was 60 kV and the capillary 
diameter 0.6 mm. 



A possible explanation for the observed increase in density vs. increasing length could 
be the following. Since the discharge is based on an avalanche effect, where each electron 
suffers a number of collisions along its path from one electrode to the other, a lengthier 
capillary would imply that each individual electron is involved in more collisions. In a short 
capillary the electrons diffuse out quickly before colliding with enough other particles to 
increase the plasma density. Then, after some optimal length is crossed, it is tougher to ionize 
the same number of particles per unit volume because the acceleration and speed of the 
electrons die out before they reach from one electrode to the other. 

Figure 3.7 shows the temperature dependence for the same capillary lengths as the 
previous graph. The temperature was extracted from line to background intensity ratio of the Ha 
line. It is interesting to observe that the temperature also increases with length, and the 
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Plasma Density vs Gas Pressure 
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Figure 3.8: The plasma density as a function of the neutral gas pressure for three different charging 
voltages. The capillary length is 20 nnnn and the dianneter 0.6 nnnn. 



maximum temperature is reached at 16 mm, while the maximum density was reached at 
12 mm. In fact, the absolute linewidth was larger for the 16 mm capillary, but according to these 
results the large linewidth was related to a large temperature increase as opposed to a plasma 
density increase. Since these plasma densities are close to lO^^cm ^ the temperature affects the 
linewidth more strongly as compared to low densities where the temperature effect is not as 
important. 

3.3.3 The Effect of the Charging Voltage 

Figure 3.8 shows plasma density vs. pressure for different charging voltages, 54 kV, 66 kV and 
78 kV, for a 20 mm long capillary with 0.6 mm diameter. 

The plasma density for a given pressure is slightly higher when the charging voltage is 
increased from 54 kV to 66 kV, but there is no significant increase when the voltage is further 
increased to 78 kV (except at low pressures). The extra energy that is deposited into the plasma 
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in the latter case is most lil<ely distributed into the kinetic energy of the plasma electrons. This 
conclusion is drawn because the plasma temperature was measured to increase from an 
average of 2 eV to an average of 4 eV as the voltage was increased from 66 kV to 78 kV (not 
shown here). Thus further increasing the charging voltage in this case is not expected to yield 
analogous increased plasma density. 

It should be noted here that the energy that is deposited to the plasma during a 
discharge is not controlled solely from the external power supply. The plasma has its own 
impedance which is variable depending on the initial pressure, the plasma density and the 
plasma temperature. The energy transfer would be maximized if the internal impedance of the 
power supply was adjusted to match the resulting plasma impedance. A mismatched circuit is 
identified by the large number of oscillations in the current that flows in the system as a 
function of time. 

3.3.4 The Effect of the Capillary Tube Diameter 

Figure 3.9 shows plasma density vs. pressure for 3 different capillary diameters (0.6 mm, 1.2 mm 
and 3 mm). The capillaries were 20 mm long and the charging voltage was set to 60 kV. 

The data shows that narrower capillaries greatly enhance the plasma density. When the 
diameter decreases 2.5 times from 3 mm down to 1.2 mm the density increases 3.25 times, 
from 0.8x10^^ cm'^ to 2.6x10^^ cm'^. Then, when the diameter is decreased further by a factor 
of 2, the density jumps to 7x10^^ cm'^ which is an increase by 2.7 times. This is probably due to 
the enhanced current density. Assuming that the current supplied is approximately same over 
those discharges, ideally the current density should be inversely proportional to the square of 
the diameter of the capillary. The enhancement observed in density is not linear to the current 
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Plasma Density for various capillary tube 

diameters 
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Figure 3.9: The dependence of the plasma density as a function of pressure for 3 different capillary tube 
diameters. The capillary lengths are 20 mm and the charging voltage 60 kV. 



density but it is still very strong, suggesting that in order to achieve high one must minimize the 
capillary diameter (while preserving the breakdown). 

Additionally, when a 0.3 mm diameter capillary of the same length was tested, it was 
found that the spectrum was radically different with a very large spectral background and no 
clear Ha line visible. This was attributed to the very high temperature of the plasma caused by 
the mismatch in the impedances between the high-voltage pulser and the plasma. Only when a 
50 load was added in series to the capillary did the temperature drop and the Ha line 
reappeared (see also the data in the next section). 

3.3.5 The Effect of the Hollow Electrode Inner Diameter 

In this setup we tested two pairs of electrodes, one pair with 0.6 mm hole diameter (equal to 
the diameter of the capillary in most cases) and one pair with no holes at all (flat). The purpose 
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of this experiment was to test for any hollow cathode effects that may enhance the plasma 
density in the capillary. The capillary that was tested was 16 mm long. We observed that the 
plasma density is relatively insensitive to the electrode shape, as the differences never exceeded 
the shot-to-shot statistical variation. Temperature measurements similarly showed that the 
plasma temperature is also not affected significantly by the choice of the electrode. 

3-4 Variation of the Plasma Density along the Length of the 
Capillary 

For the resonant experiments described later in this dissertation, a major requirement for the 
plasma source is that the plasma remains uniform at a constant density over the length of the 
capillary. The tolerance in the uniformity of the plasma is determined by the number of bunches 
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Figure 3.11: Plasma density along the capillary axis for a 60 kV discharge in a 17 nnnn long capillary at 
1 atnn of pressure. The error bars indicate 1-a variations over 10 discharges. 



N that drive the wakefields, and scales roughly as 1/ N . For example, for A^ = 5-10 drive 
bunches (typical for the experiments under consideration) the plasma density needs to be 
constant within 10 - 20%, otherwise the excited wakefields from each bunch will be out of 
phase and cancel each other before the bunches leave the plasma. Therefore, it was important 
to establish whether the capillary discharges can produce longitudinally uniform plasmas. 

In order to measure the longitudinal dependence of the plasma density an optical fiber 
was mounted in a 2 dimensional translation stage and it was placed directly against the outer 
side of the glass capillary walls. The time-averaged plasma light was collected in successive 
discharges by tuning the position of the fiber along the length of the capillary. 

It was measured that for the given 17 mm capillary shown in Figure 3.11, the plasma 
density over a length of 8 mm near the center of the capillary did not vary more than 10% from 
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the average value. However near the cathode side the density differs by a factor of two 
compared to the center of the capillary. This suggests that the resonant interaction length is 
shorter than the quoted length for the capillary. The plasma temperature behaves similarly 
along the capillary, having a flat region near the center and dropping fast near the electrodes. 
Similar results that show reduction in the plasma density and temperature near the electrodes 
were reported by Levin et al. in [108]. 

3-5 Variation of the Plasma Density with Time 

The time-integrated measurements presented in the previous two sections highlight the effects 
of the various capillary and discharge parameters, yet they are expected to yield smaller plasma 
densities than the actual ones because they are not sensitive to high densities that occur for 
only a short period of time - those peaks are averaged out. Therefore, since in many cases the 
exact plasma density is selected by tuning the arrival of the electron beam with respect to the 
discharge, it is important to resolve the plasma density in time. 

The time-dependence of the plasma density was investigated in the course of the 
multibunch experiments at ATF. We tested a 10 mm long hydrogen-filled capillary with 1 mm 
inner diameter using the following technique. The light during the discharge was collected by an 
optical fiber inserted directly into the capillary center, which was then fed into an optical 
spectrometer equipped with a 1,200 gr/mm grating. The entrance and exit slits were adjusted 
such that a 0.5 - 1 nm spectral window could be imaged onto a Photo Multiplier Tube and yield 
the signal recorded at that wavelength range as a function of time. By varying the central 
wavelength of that window, the profile of the Ha line versus time could be reconstructed after 
successive shots. By fitting a Lorentzian distribution to each set of wavelengths for every point in 
time the plasma density could then be extracted. 
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Figure 3.12: Plasma density as a function of tinne for different discharge paranneters. The plasnna 
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noisy at later tinnes where the light ennitted fronn the plasnna is weak. 



The success of this method relies on the reproducibility of the discharge among the 
shots (which needs to be monitored by collecting the light emitted in all wavelengths), and it is 
limited by the resolution of the spectrometer and the Signal-to-Noise ratio which decreases with 
narrower spectral windows. Figure 3.12 shows typical results for the evolution of the plasma 
density for different discharge parameters. The maximum density was obtained for 100 Torr of 
neutral gas pressure and 20 kV of charging voltage (two successive sets of measurements are 
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shown here). The peak density dropped by a factor of 2 when the charging voltage was reduced 
to half its initial value, 10 kV. 

For 30 Torr of gas pressure, the peak voltage is even lower, near the 10^^ cm"^ range. 
Several hundred nanoseconds after the peak, the density is seen to drop linearly which indicates 
exponential diffusion. The plasma current typically peaks 50 - 100 ns before the density and 
vanishes 200 - 400 ns after that peak. 

In those sets of data the spectral window was 1 nm in the 100 Torr cases and 0.5 nm in 
the 30 Torr cases. This change allows smaller densities to be measured (narrower linewidths), 
yet the method becomes unreliable earlier in time due to the decrease in the light recorded. It 
can be improved though if a higher-resolution spectrometer is utilized while the optical 
sensitivity is preserved. For the spectra shown in this graph, the maximum single-to-noise ratio 
was 10, which implies that the density is known with 20% accuracy or better^ at densities down 
to 10'' cml 

3-6 Chapter Conclusions 

The following items are concluded from this chapter: 

• Capillary discharges are a favorable option for PWFA experiments that require pre- 
ionized plasmas. The density can be tuned through the initial neutral gas pressure, the 
capillary geometry or by tuning the arrival of the beam with respect to the discharge. 

• The Stark Broadening of the Hydrogen line Ha is a good diagnostic for plasma densities 
in the range 10^' cm'^ - 10^^ cm'^. The plasma temperature also needs to be known at 
the high densities by measuring the line-to-continuum spectral ratio. 



' Higher densities than 10^' cm'^ are known more accurately and lower densities are known less 
accurately, since higher densities produce more light. 
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• The density increases with increasing gas pressure, albeit less than linearly proportional. 
It becomes harder to ionize the plasma unless more energy is deposited. 

• The maximum plasma density initially increases with the capillary length up to an 
optimal of 12 mm length and a density of 9.5x10^^ cm"^ for 2 atm of pressure. The 
density decreases for longer capillaries and eventually it becomes hard to even break 
down. 

• Increased charging voltage increases the amount of light emitted but this is attributed to 
temperature increase and not due to density increase. 

• The density increases between linearly and quadratically with decreased capillary tube 
diameter due to higher current density. 

• The hollow cathode effect is insignificant at the pressure range measured here 
(>100Torr). 

• The plasma density varies less than 10% longitudinally near the center of the capillary 
(for about half the total length), but it decreases by as much as a factor of 2 near the 
electrodes. 

• As a function of time the density has a peak shortly after the discharge current peak and 
then eventually starts dropping exponentially, probably due to diffusion as the plasma 
recombines at the capillary walls. 



66 



The things we know best are the things we haven't been taught. 

- Marquis de Vauvenargues 



Chapter 4 

High-gradient Acceleration of a 

Trailing Electron Bunch 



Plasma waves can sustain extremely large electric fields that are orders-of-magnitude larger 
than those in conventional radio-frequency accelerators, which are limited by vacuum 
breakdown to accelerating gradients of up to 150 MV/m [23]. Such large amplitude electron 
density waves, or wakes, can be excited in plasmas by a laser pulse (laser wakefield acceleration 
- LWFA), or a relativistic particle beam (plasma wakefield acceleration - PWFA). Recent LWFA 
experiments demonstrated quasi-monoenergetic acceleration of self-trapped plasma electrons 
[58, 63, 106, 120]. Further scaling of LWFA to higher energies, by using higher laser power but 
larger spot sizes and lower density plasmas, will probably require injecting relativistic electron 
bunches into a plasma wave rather than starting with plasma electrons at zero energy (in order 
to maintain small emmitances). So far, this approach to potentially monoenergetic particle 
acceleration using plasmas has not been explored. 
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Experiments using PWFA methods similarly face the challenge of producing low-energy- 
spread acceleration of an injected relativistic particle bunch. In previous PWFA experiments (see 
for example Refs. [12, 20, 179]) a single relativistic electron-bunch both drove the wake and 
provided the electrons to be accelerated. Using this scheme at SLAC, a record-high energy gain 
of 42 GeV over 85 cm of plasma was demonstrated [20], albeit with an undesirable -100% 
spread of the electron energy spectrum. To realize a future collider, such as one incorporating 
the PWFA afterburner concept [105], a well-defined bunch is required that is suitably phased on 
the plasma wake of a preceding drive bunch thereby to achieve high efficiency and a small 
energy spread (~0.1% which is typical for conventional accelerators). In the following sections a 
double-bunch PWFA experiment is detailed which demonstrates controllable high-gradient 
acceleration of a witness relativistic electron bunch injected into a plasma wave. 

4,1 Previous Work 

Earlier double-bunch PWFA experiments [14, 100, 128, 134, 148] utilized relatively long, 
picosecond electron bunches in low-density (10^^ cm'^) plasmas; wakefields of only up to 
~4 MV/m were inferred. This experiment differed substantially from these studies. First, the 
driver and witness bunches had subpicosecond lengths (cr^ - 100 fs), and both were shorter 
than the plasma wavelength. Consequently, the energy shift of both bunches could be directly 
observed rather than mathematically extracted^ as was required in earlier works. Second, the 
shorter bunch lengths and the higher plasma densities employed (up to 10^^ cm'^) resulted in 
generated wakefield amplitudes that were two orders-of-magnitude larger compared to those 
studies. 



^ The wakefield in some of these works was extracted by evaluating the shifts of the bunch centroid after 
background signal subtraction. 
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Figure 4.1: Layout of the double-bunch PWFA experimental setup. Figure by W. Kimura. 

In this experiment, the plasma density could be chosen such that energy lost by the drive 
electron bunch was transferred through the plasma wave to the second bunch, which, in turn, 
gained energy with minimum energy spread. The drive bunch lost about 1.0 MeV over 6 mm 
propagation in a plasma of 10^^ cm'^ density, and the witness bunch, delayed by 500 fs, gained 
-0.9 MeV corresponding to an average loaded accelerating gradient of -150 MV/m, while the 
beam loading efficiency was estimated to be approximately 77%. The measured energy gain and 
loss agree well with 2D linear theory calculations. This experiment was the first to generate and 
directly probe large plasma accelerating gradients (> 100 MV/m) utilizing a trailing electron 
bunch [85]. 

4,2 Experimental Setup and Diagnostics 

The experiment was performed at Brookhaven National Laboratory's (BNL's) Accelerator Test 
Facility (ATF). A photocathode rf gun followed by a conventional 2.856 GHz (S-band) accelerator 
produced a 1.5 ps-long (rms), 500 pC, 60 MeV single electron bunch [28] that was compressed 
and split into two distinct (in time and energy) subpicosecond bunches after traveling through a 
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Figure 4.2: Example of three different double bunch energy spectra taken nnany nninutes apart, indicating 
the stability of the beann break-up. 



chicane compressor and "dog-leg" dipoles downstream from the linac (see Figure 4.1). These 
two bunches were focused transversely to cr^ -100 jum at the entrance of a 10^^ - 10^^ cm"^ 
density plasma produced by an ablative capillary discharge [83]. A magnetic spectrometer at the 
end of the beamline recorded the energy change imparted to the bunches by the plasma. 

4.2.1 Coherent Transition Radiation Diagnostics 

Production of subpicosecond driver-witness electron bunch pairs separated by a distance 
appropriate for PWFA experiments in high-density plasmas (10^^ - 10^^ cm'^) is challenging 
[155]. The breakup of the electron beam into two distinct bunches at ATF was attributed to the 
combination of a nonlinear energy chirp introduced by the linac, over-compression of the bunch 
in the chicane, and coherent synchrotron radiation effects in the chicane and the dog-leg dipole 
magnets. 

Therefore, although the breakup of the bunch is repeatable and consistent [98, 165], 
control of its characteristics is limited. The two bunches were separated by approximately 1.8 
MeV in energy and had a typical full-width-at-half-maximum (FWHM) of 0.4 MeV, thus allowing 
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Figure 4.3: Experimental setup of the CTR interferometry diagnostic for the double-bunch experinnent. 
Figure by W. Kinnura. 

direct observation of their energy shifts caused by their interaction with the plasma. Figure 4.2 

shows three double-bunch energy spectra recorded many minutes apart, demonstrating the 

stability of the break up. 

In addition, coherent transition radiation (CTR) interferometry [125, 157] was used to 

diagnose the two bunches in time (see details in Appendix C.4). The experimental setup for this 

diagnostic is shown in Figure 4.3. Figure 4.4 shows the CTR signal recorded as the position of the 

mirror in one of the interferometer arms was varied, for the cases of a single bunch and both 

bunches. The dips observed in the signal are attributed to the frequency limitations of the 

systems which act as a high-pass filter and neglect frequencies below a certain threshold^ 

Assuming that each bunch has a Gaussian shape, by fitting the observed signal traces to the 

model it was found that the high-energy and the low-energy bunches, respectively, were 

roughly 45 jum (150 fs) and 27 jum (90 fs) long (rms), and were separated by Az=l50jum 

(500 fs). Beam-position monitors and a Faraday cup also indicated that the total charge was 



^ This frequency cutoff was included in the CTR model. 
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Figure 4.4: CTR interferometry diagnostic for the single bunch (left) and for both bunches (right). The blue 
circles indicate experinnental data, while the solid red line indicates the nnodel fit. Figure by W. Kinnura. 



preserved during the break up, with the low-energy bunch having about 60% (180 pC) of the 
charge of the high-energy bunch (300 pC). Finally, observations of the interaction between the 
bunches and the plasma showed that the high-energy bunch lost energy in the plasma 
independently of the presence of the low-energy bunch, thus implying that the former preceded 
the latter, which acted as the witness bunch. 

4.2.2 Plasma Density Diagnostics 

The plasma source was a 16.5 kV pulsed electrical discharge through a readily available 6 mm 
long polypropylene capillary with a 1 mm inner diameter. It should be noted that there is no 
fundamental limitation in using a longer, centimeter-scale, capillary to reach a larger energy 
gain. The breakdown in the vacuum ablates the capillary walls creating a carbon-hydrogen 
plasma. The discharge light was collected and guided by an optical fiber into a spectrometer that 
measured the hydrogen Balmer Ha linewidth, and the plasma density was derived from the well- 
tabulated Stark broadening [7] (see section 3.2). An intensified time-gated camera collected 
spectral measurements with a resolution of 300 ns. 
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Figure 4.5: Plasma density diagnostic for the double-bunch experinnent. The Stark broadening of the Ha 
line was used to identify the density in the first |is (solid line), after which the density was extrapolated 
assunning exponential diffusion. The two dashed lines indicate the density dependence for two different 
tinne constants. 



Those measurements are plotted in Figure 4.5 for the first \xs, after which the density 
was extrapolated due to insufficient plasma light reaching the camera. The data show that the 
plasma reaches a maximum density of approximately 5x10^^ cm ^ after which it is assumed to 
exponentially decay due to diffusion through the capillary openings [21, 179], according to the 

formula n(t) = nQe ^ . The diffusion constant chosen was estimated to be equal to T = 03jus 
for the data sets in this experiment, which corresponds to a nominal diffusion rate of one order- 
of-magnitude per 1.5 jus. Therefore, the required plasma density, n^, could be selected by 

tuning the time-delay of the e-beam after the discharge started. The peak electric field of the 
discharge was about 3 MV/m (16.5 kV over 6 mm), which was much smaller than the typical 



73 



plasma wakefields in this experiment (> 80MV/m), and did not affect the energy spectra of the 
electron bunches significantly. 

4,3 Plasma Interaction Results 

4.3.1 Double Bunch Interaction 

The peak beam density, n^, of the two compressed bunches was approximately 10^^ cm ^ so 

placing the physics of the beam-plasma interaction mostly in the linear overdense regime ( 
^b^^p) fo'' the range of plasma densities used in these experiments [116]. In this regime, each 

bunch independently creates a wakefield with a spatial period equal to the plasma wavelength. 
While the drive bunch only loses energy due to its own wakefield, the witness bunch samples 
the superposition of the two fields, and loses or gains energy depending on its relative phase in 
the wakefield driven by the first bunch. This phase was controlled by tuning the plasma density 
(and hence, the plasma wavelength). Specifically, energy gain was expected at plasma densities 
such that the witness bunch samples the second half-period of the plasma wave driven by the 



first bunch, i.e., when approximately — ^ < Az < yl ; here, X = with co = \ . 

In order to model the behavior of the witness bunch, a 2D numerical program was 
developed to calculate the linear wakefields driven by two Gaussian-shaped bunches [86] with 
the aforementioned experimental parameters (see Appendix D on ThemOsiris code). The 
validity of the modeling was verified by comparing the findings with the fully-explicit 2D particle- 
in-cell code OSIRIS [71]. The amplitudes of the wakefields simulated by both methods agreed 
within 10%, thus justifying the use of the thousand-times-faster numerical calculation. 
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Figure 4.6: Experimental and sinnulated energy spectra of the double-bunch beann after the 6 nnnn-long 
capillary discharge at 4x10^^ cnn'^ plasnna density (left colunnn), and at 1x10^^ cnn'^ density (right colunnn): 
a,f) raw energy spectrunn without plasnna; b,g) raw energy spectrunn with plasnna on; c,h) experinnental 
energy profiles; d,i) sinnulated energy profiles; e,k) sinnulated plasnna wakefield and position of the 
bunches inside the wake. 
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Comparison of the double-bunch interaction at two different densities is shown in 
Figure 4.6. The left-hand side column presents experimental energy spectra of the e-beam with 
the plasma discharge turned off [Figure 4.6(a)] and as measured after the 4x10^^ cm'^ density 
plasma [Figure 4.6(b)], together with plots of the area-normalized line profiles integrated 
vertically (across the beam's profile) in Figure 4.6(c). At this plasma density, both bunches 
clearly lost energy, since the plasma wavelength is longer than the separation between the 
bunches (\ =529jum>2Az-300jum). The drive bunch lost ~0.5 MeV, and the witness bunch 

lost -0.8 MeV, corresponding to average wakefield amplitudes (over the 6 mm capillary length) 
of 83 MV/m, and 133 MV/m. Figure 4.6(d) plots the simulations of the predicted energy spectra 
before and after the plasma interaction; they confirm that both bunches should lose energy. 
Figure 4.6(e) shows simulated combined plasma wakefield and longitudinal phasing of the two 
bunches in this field. 

The right-hand side column in Figure 4.6 shows the experimental energy spectra 
recorded when the plasma density was increased to IxlO^^cm ^ therefore reducing the plasma 
wavelength to/l^ =334jUm-2Az . The witness bunch in this case mainly gained ~0.9 MeV in 

energy [Figure 4.6(g-h)]. The simulation illustrates that the observed peak in the energy 
spectrum around 59 MeV resulted from the superposition of both accelerated electrons from 
the witness bunch and decelerated electrons from the drive bunch (which only lost energy) 
[Figure 4.6(i-k)]. Some residual charge of the drive bunch was recorded around 60.5 MeV, 
probably reflecting the non-Gaussian initial energy distribution in the experiment. 
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Figure 4.7: Experimental and sinnulated energy spectra after the 6 nnnn-long capillary discharge at 1x10^^ 
cnn'^ plasnna density for the double-bunch beann (left colunnn), and for the witness bunch only (right 
colunnn): a,d) experinnental energy profiles; b,e) sinnulated energy profiles; c,f) sinnulated plasnna 
Wakefield and position of the bunches inside the wake. 



4.3.2 Single Bunch Interaction 

In order to confirm that the witness bunch was indeed affected by the drive bunch, the higher- 
energy bunch was partially blocked by closing a slit located at the dispersion plane inside the 
dog-leg. Side-by-side comparison of the double-bunch and single-bunch interaction is shown in 
Figure 4.7. When the plasma density remained at 1x10^^ cm ^ the witness bunch had an average 
energy loss of ~1 MeV [Figure 4.7(d)]. Since the witness bunch lost 1 MeV due to its own wake 
and gained 0.9 MeV in the presence of the drive bunch, we conclude that the net energy shift 
due to the drive bunch was 1.9 MeV, corresponding to a -315 MV/m unloaded accelerating 
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Figure 4.8: Experimental data points for the energy shift of the witness electron-bunch centroid for a 
range of plasnna densities. The solid curve represents 2D nunnerical calculations for the centroid energy 
shift. 



Wakefield amplitude driven by the first bunch. This value agrees with the numerical calculation 



150 
and therefore the beam-loading efficiency is estimated approximately rj^ -1 



■ 11% for 



this interaction^ 

4.3.3 Plasma Density Scan 

Figure 4.8 depicts the measured energy shift of the witness-bunch centroid as the plasma 
density was scanned from 1x10^^ cm'^ up to 4x10^^ cm"^. Comparison of the recorded shifts with 
the 2D numerical calculation of the corresponding predicted energy shifts demonstrates good 
agreement with the simulation for the chosen plasma density decay profile. The experimental 
points show that the witness bunch progressively lost more energy as the plasma density was 



^ That would be the fraction of the energy transferred if the bunches were let to propagate over long 
distance while preserving their wakes. See section 2.5. 
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increased from its minimum value, provided that the plasma wavelength remained longer than 
the bunch separation ( /l^ > Az ) and both bunches reside at the decelerating phase of the wake 

[see Figure 4.6(e)]. As the density was increased further near the 10^^ cm'^ range, the plasma 
wavelength decreased and the witness bunch now sampled the second, accelerating half of the 
plasma wave period [/l^ -2Az; see also Figure 4.7(c)], and therefore it gained energy while 

simultaneously loading the wake since its charge was comparable to that of the drive bunch. 
Energy loss was observed again at an even higher plasma density (> 4x10^^ cm'^) where the 
witness bunch sampled the decelerating first half of the second period of the excited plasma 
wake (Xp<Az). 

On the low plasma density side the wakefield eventually drops when the transverse 
beam size of the bunches (a^ ~100//m) becomes shorter than a plasma wavelength because 
the total charge in the bunches is not enough to excite a plasma oscillation with the 
corresponding large plasma wavelength. On the high plasma density side, the wakefield also 
drops because the longitudinal size of the bunches ( a^ ~ 30/um) becomes long when compared 
to the plasma wavelength and thus the wakefields cancel out inside the bunches. For the plasma 
densities in between these two extreme cases, the wakefield excitation is optimal and its 
amplitude is determined by the relative phase at which the wakefields excited from each bunch 
add together. 

In summary, it was demonstrated that by adjusting the plasma density tunable energy 
gain or energy loss of a short ( O"^ < /l^ ) witness bunch that samples the wakefield generated by 

a leading subpicosecond drive bunch is possible. The measured energy gain of 0.9 MeV over a 
6 mm long, 10^^ cm'^ density plasma corresponds to an unloaded average accelerating gradient 
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of 315 MV/m. The loaded gradient measured was 150 MV/m and the extrapolated beam loading 
efficiency of 77%. This work demonstrates that short-pulse accelerated trailing electron beams 
can be produced with a double-bunch PWFA scheme, and hence, the method constitutes a 
promising step forward towards producing monoenergetic bunches in next-generation ultra- 
high gradient plasma accelerators. 

4,4 Chapter Conclusions 

The following items are concluded from this chapter: 

• Two -100 fs long 500 pC bunches could be generated at ATF, distinct in both energy and 
time (as diagnosed through CTR). Both bunches are shorter than half plasma wavelength 
for densities 10^^ cm"^ and below. 

• By tuning the plasma density, the second bunch could be made to sample the 
accelerating phase of the wakefield driven by the first bunch. Loaded accelerating 
gradients of 150 MV/m were measured for the first time for a trailing short bunch. 

• When the drive bunch was blocked the witness bunch lost 1 MeV energy at the same 
density it gained 1 MeV before, thus implying a 315 MV/m unloaded accelerating 
gradient and an estimated 77% beam loading efficiency. 

• The energy shift of the witness bunch as a function of the plasma density (selected 
through the timed arrival of the beam with respect to the plasma discharge) agrees well 
with a simple 2D linear theory model. 
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For every complex problem there is an answer that is clear, 
simple, and wrong. 

- Henry Louis Mencken 



Chapter 5 

Multibunch Schemes and 

Simulations 

In this chapter the theoretical framework of designing a multibunch PWFA accelerator will be 
analyzed. The difference compared to previous chapters is that now the wakefield is driven by 
more than one bunches (a witness bunch is still required to sample the generated wakefield). 
When using multiple bunches, by controlling their position and charge the accelerator can be 
tuned to either maximize the wakefield, the transformer ratio or the efficiency of the system. 
Specifically of interest for future PWFA-based colliders are schemes where the energy of the 
incoming witness beam is multiplied. In section 5.1 the various schemes will be analyzed in the 
linear regime, and particle-in-cell simulations using OSIRIS will be shown in section 5.2. In 
section 5.3 we will simulate an example of multiple drive bunches in the nonlinear regime. 

5-1 Multibunch Schemes 

In this section the basic principles of multibunch PWFA will be analyzed. Specifically, we will 
explore how, by adjusting the bunch spacing and the charge per bunch, that the plasma 
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accelerator can be tuned to maximize the wakefield, the transformer ratio or the efficiency. In 
the models we will utilize transversely Gaussian and longitudinally square bunches of the same 
width (with the exception of the witness bunch which may be shorter). Although longitudinally 
Gaussian shapes are more realistic, the square bunches do not have significantly different 
behavior (see Figure 2.4) yet they are much more tractable analytically and provide quick and 
easy physics insight in the process. 

So let us assume a series of M drive bunches with the same width W. The m — th 
bunch has N^ particles, is placed at a position ^^ and has a transverse spot size a^ . This 
corresponds to the following beam density profile: 






(49) 



The function rect has a value of 1 for |<^-<^^|< w, y^ for |<^-(^J = w and otherwise. The total 

longitudinal electric field is simply the superposition of the fields created by the individual 
bunches, which was given in equation (27), and is equal to 



£,(r,#) = +^X^sin[^,(#-#,+w/2)].r.c? 
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The transverse wakefield can be found from equation (10): 
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(50) 



Xn,sin[^4|-|J].f/(#-#^.-w/2).- (51) 
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In the next few sections we will explicitly describe what can be achieved with different setups. 

5.1.1 Enhancing the Wakefield 

In order to create the maximum possible wakefield, the individual wakefields from each bunch 
need to be added in-phase. This occurs when the bunches are separated by one plasma 
wavelength^ placing each bunch in the decelerating phase of the wake created by previous 
bunches. This means that 



L={m-^K 



The wakefield left after the bunches (^>^^ +w/2) is then equal to (from (50)) 



p f k w^ ^ 

E,{r,^) = + ^2sm -^ cos (fc/)- 7? (r).£n,. 

^0% V ^ y i=l 



(52) 



(53) 



It is interesting to observe in this case that the number of particles per bunch does not affect the 
wakefield, as it simply adds together and the total wake is proportional to the total number of 
particles in all bunches. The focusing force at the center of bunch m (<^ = <^^ =(m-l)/L^) is 



equal to 



W.(^^J = - 



^0% 



'^kpW^ ^-1 



Z^ 



V ^ y 



j=i 



dR 

dr 



(54) 



We observe that in addition to its own focusing force, each bunch will experience the force due 
to all the previous bunches. Thus, later bunches are focused stronger. 

The tradeoff for the enhanced wakefield in this scheme is the reduction in the 
transformer ratio, which approaches unity as the number of bunches increases. For identical 
bunches (N^=Nq) with transformer ratio R^ , the total transformer ratio R^^^ is equal to 



^ More general, the separation can be any integer number of plasma wavelengths. 
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Figure 5.1: Example of a multibunch plasma accelerator in a maximum wakefield setup. The bunches are 
identical with k^w = 7r and are separated by one plasma wavelength apart. A witness bunch 180° out of 

phase samples the accelerating wakefield. In this example A^ = 250jum and a^ = lOOjum . 
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An example of a PWFA under such a setup is presented in Figure 5.1 where the on-axis beam 
density and wakefield of a 5+1 bunch system (5x30 pC drive bunches and 1 witness bunch) is 
shown. In this example kpW = 7r , and the bunches are separated by 250 jU/r? at 1.8x10^^ cm'^ 

plasma density (/l^ =250 jum). A relevant experiment that creates the bunches through IFEL 

using a laser is analyzed later in section 6.1. 

The reason why the transformer ratio and the efficiency are not enhanced is because 
the last bunch feels the almost full decelerating field, while the early bunches feel a smaller 
decelerating field. In this example the transformer ratio is R = l.l and the efficiency is only 
13% (as defined in section 2.4.1). In order to increase the transformer ratio each bunch needs 
to feel the exact same decelerating field, and this is explored in the next section. 



84 



5.1.2 Enhancing the Transformer Ratio 

In order to maximize the transformer ratio, each bunch needs to be placed in the accelerating 
phase of the wake created by the preceding bunches, and its charge needs to be increased such 
that the total decelerating wakefield under each bunch is the same for every bunch. This 
scheme was suggested first by Laziev et al. [104] and more recently by Power et al. for a 
dielectric wakefield accelerator [144] and was demonstrated experimentally in a proof of 
principle experiment in 2006 by Jing et al. using two drive bunches [78, 79]. 
The total wakefield under the m-th bunch is given by equation (50) 



eA^^^) = 



^0^/7 



«*.sin[^^(^-4+w/2)] + 2sin ^ T% cos[k ^ (^ - ^ .)] 



\ ^ J 



7=1 



■R{r) 



(56) 



^0^/7 



«Msin[^,(^-4+w/2)].i?(r) 



The last equality is set such that the wakefield under each bunch is equal to the wakefield under 
the first bunch (shifted by the location of the bunch center). By requiring the above expression 
to hold true for every ^ , a solution for the number of particles N^ for the m-th bunch and 

its center ^^ can be found. For example, for an arbitrary bunch width k^w , a solution for the 

number of particles and the location for the second bunch is 



N,^N, 



l + 2sin 



V 2 , 



57t kw 



(57) 



^^=^p^i+ — + ^ 



Under an arbitrary bunch width the bunches are not equidistant and their charge may not scale 
in a simple way - those values need to be evaluated numerically from equation (56). 

The most interesting case though is the one that maximizes the transformer ratio of the 
single bunch, i.e. k w = 7r for which R^=2. In this case the wake around each bunch is 
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symmetric with respect to its center and the equations are greatly simplified. In that special 

scenario it is found that the number of particles per bunch and their location scale in a simple 

way: 

N^=N,{2m-\) 

(58) 

The decelerating wakefield under each bunch is now the same, and the envelope of the beam 
density of the bunches is a linear ramp, much like the continuous single-bunch linear ramp that 
is analyzed in Appendix B.2 . Equation (58) implies that the bunches must be placed 1.5 plasma 

wavelengths apart from each other^ and their charge must scale as A^^ x [l : 3 : 5 : ...] . 

The total wakefield left after the drive bunches in the case where equations (58) is 
satisfied for the k^w^n case is 

£, (r,^) = M . ^(-if -' cos(^/) . R{r) (59) 

Note that the phase of the wakefield resets after each drive bunch. Using equation (56), the 
transformer ratio of the system is then 

R^^^=2M=M'R, (60) 

The last equation summarizes the ultimate advantage of this scheme, since the transformer 
ratio of the system is proportional to the number of bunches used to drive the wakefield. 

The transverse wakefield at the center of the m-th bunch is in this scenario equal to 
(from equation (51)) 



^ They can also be any half-integer number of plasma wavelengths apart. 
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Figure 5.2: Example of a multibunch plasma accelerator in a maximum transformer ratio setup. The 
bunches are identical with k^w = n and are separated by 1.5 plasma wavelengths apart. The total charge 
is 500 pC and in each bunch is increased linearly. A witness bunch equidistant to the other bunches 
samples the accelerating Wakefield. In this example X^ = 250//m and cr^ =lOOjum. 



^.(^^J = -^{(--ir+sin(^^(^-^j)].^ 



(61) 



Here \(^ -(^^\<w/2. It is interesting to observe that although the system is designed such that 

the longitudinal wakefield is identical under each bunch, as equation (56) suggests, on the other 
hand the transverse wakefield is not the same for each bunch but it scales at their center 
quadratically with the number of drive bunches instead. This is a consequence of the fact that 
the total charge fed into the system scales also quadratically with the number of bunches, i.e. 

A^^^^ = M 'N^. Ass result, later bunches experience stronger focusing force. 

Figure 5.2 illustrates an example where a ramped bunch train of 4 drive bunches is fed 
into a 1.8x10^^ cm'^ plasma. The bunches are separated by 1.5 plasma wavelengths apart and 

have a total charge of 500 pC, while their charges scale as 31pCx[l: 3:5:7]. The transformer 



87 



ratio is 7.9 and the ID driving efficiency of the system is 64%, equal to the efficiency of a single 
bunch. 

Several things need to be pointed out in this figure. Although the total wakefield 
increases linearly with the number of bunches, similar to the maximum wakefield scheme of 
section 5.1.1, this occurs here by supplying quadratically more charge (as opposed to linearly 
more). The total wakefield amplitude per total unit charge put into the system scales inversely 
proportional with the number of bunches in the ramped bunch train scheme described here 
[87]. In addition, from an energy perspective, each particle is depositing into the plasma as 
much energy as the particles in the first bunch do, and the rest of its energy is being transferred 
to the following bunch in the train to prevent it from decelerating faster. In this way all the 
particles deposit their energy at the same rate (on average along the bunch), and the last bunch 
in the train releases this net energy into the plasma. 

Note that the efficiency of the system is limited by the efficiency of a single bunch, 
which is achieved when the charge in each bunch is adjusted properly such that the wakefield 
under each bunch is the same. Even so, the variation of the wakefield inside each individual 
bunch sets an upper bound for the efficiency. This limit can be overcome (at the slight expense 
of the transformer ratio) if the positions and number of particles per bunch are fine-tuned, as 
we will show in the following section. 

5.1.3 Enhancing the Efficiency 

Gaussian, square and other realistically-shaped bunches have maximum single bunch 
efficiencies around 60% - 70%, depending on their width relative to the plasma density. This is 
the maximum fraction of their energy that can be transferred to the plasma. Even if the 
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Figure 5.3: Example of a multibunch plasma accelerator in a maximum efficiency setup. The bunches are 
identical with kpW = 0.567r and are separated by 1.28 plasma wavelengths apart. The charge in each 

bunch is 17.5pCx[l: 2.70:5.20:8.22j. a witness bunch samples the accelerating wakefield. In this 
example A = 250 jum and a^ = 100//m . The driving efficiency of this system is 84% and the transformer 
ratio is 5.14. 



transformer ratio is maximized, as was explored in section 5.1.2, the efficiency is never larger 
than the single bunch efficiency for those setups. 

It is possible however using multiple drive bunches to obtain a total system efficiency 
that is larger than the individual bunch efficiency, without a need to specially shape the 
bunches. This is achieved by placing the bunches at those phases of the wakefields such that the 
total wakefield under each bunch is most flat (thereby depositing energy more efficiently), even 
if the wakefield is not exactly the same between bunches. Notice in Figure 5.2 for example, 
there are particles inside each bunch that experience an almost zero wakefield (near the edges), 
while other particles near the center of the bunches experience the full decelerating field. This 
can be avoided if the bunches are placed earlier in the wakefield such that they sample both a 
portion of the decelerating phase of the wakefield and a portion of the accelerating phase as 
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well, with the number of particles per bunch is adjusted accordingly. Such an example using 4 
drive bunches is shown in Figure 5.3. 

In this scenario the bunches are equidistant and separated by 1.28 plasma wavelengths 

apart, and the total charge is 300 pC and scales as 17.5/7Cx[l:2.70:5.20:8.22] in each bunch. 

Those optimal values are found numerically by running an optimization routine for the efficiency 
as a function of the bunch separation and the bunch charges. Further enhancements (few %) 
can be achieved if the restriction for equidistant bunches is lifted. The ID driving efficiency of 
the plasma accelerator shown here is 84%, while the transformer ratio is 5.14 (instead of the 
theoretical maximum of 8 for 4 bunches). 

It is interesting to observe the physics of this interaction. In contrast to the maximum 
transformer ratio scheme of the previous section, the wake under each bunch is not identical 
anymore. The first bunch experiences a weaker decelerating wakefield for example. However 
this is allowed here because from an energy perspective the first bunch is not as important as 
the later bunches that carry more energy and affect the efficiency more significantly. In essence, 
the first bunch jump-starts the wake such that the second bunch will sample a more flat portion 
of it. Shorter bunches increase the efficiency further as they sample narrower portions of the 
sinusoidal wakes. 

5-2 Meter-scale Simulations of a Multibunch Accelerator 

So far in this chapter we have been limited to the initial wakes generated by multibunch 
electron drivers. These waves are valid at time t = 0, before the beam starts to evolve inside the 
plasma. In reality though the plasma will generate transverse forces that will act on the beam 
and unless the beam is matched the longitudinal wakefield structure will be disrupted after a 
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few mm of propagation in the plasma. In the blowout regime the beam is more stable along the 
plasma propagation because the lack of electrons creates a radially linear focusing force to 
which the beam emmitance or thermal spread can be matched, however those longitudinal 
wakefields excited cannot be efficient for multibunch drivers (see section 5.3). It is the purpose 
of this section to show that in the linear regime, even though (ironically) the transverse forces 
are not linear, a high-transformer-ratio wakefield structure can be maintained over long 
distances. The scenario analyzed in section 5.1.2 under the conditions of equation (58) will be 
explored. 

5.2.1 Multibunch Emmitance Balancing 

We seek to find a "matched" or balanced emmitance for each drive bunch such that they remain 
as stable as possible and the wakefield is not altered significantly as the bunches propagate in 
the plasma. In the linear regime examined here, and assuming that the bunches have 
transversely Gaussian shape, the matched emmitance along the bunch was analyzed in section 

2.2.2 and is given by equation (18), where the function 2(^) is the longitudinal component of 

the transverse wakefield which was found in equation (61). The matched emmitance along the 
longitudinal beam profile of the m-th bunch is then 

^.(^) = #J-[('^-l)'+sin(^^(#-4))].77(^c7^ (62) 

Under realistic terms, one might attempt to match in an average way the emmitance at the 
center of each bunch, ^ = ^^ : 

s,{t) = {m-iy4^l^:J7{^CT^ (63) 
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The value of l{o'r) '^ given for transversely Gaussian bunches in Figure 2.2. It should be noted 
that the value given from the above formula is approximate, as the bunch will eventually evolve 
and the assumptions under which this formula was derived (e.g. transversely Gaussian bunches) 
will not hold. Still, it serves as a good starting point. Note that all such transverse nonlinear 
focusing effects can be avoided if a hollow plasma channel is utilized [34], thus avoiding the 
need to adjust the emmitances. 

5.2.2 Simulation Results 

The parameters of the multibunch accelerator that is modeled in 2D OSIRIS are as follows. The 
plasma has a neutral density of 2x10^^ cm"^ and a length of 50 cm. M =4 drive bunches have 
longitudinal widths of kpW = 7r (37.5 lum), transverse spot sizes of k^a^ =2.4 (28.6 lum), initial 

energy of ;i^ = 2xl0 (102.25 GeV), and are separated by 1.5 plasma wavelengths or 112.5 jum. 
Their charge scales as 62/7Cx[l:3:5:7] for a total of ~1 nC of drive charge, which corresponds 
to normalized beam densities of n^^ / n^ =0.01x[l:3:5:7], thus placing the interaction into the 
linear regime. The transverse emmitances of the drive bunches where tuned according to 
equation (63) to the values £^ I cr^ =22:22:45:67 \ A short and narrow test witness bunch 

with 5 pC charge, width k^w^^nlA (9.4jU/r7), transverse spot size kp(j^=0.2 (2A jum) and 

same energy is placed another 1.5 plasma wavelengths after the last drive bunch, at the peak of 
the accelerating phase of the wake ( n^^ I %= 0.3 for that bunch). 



^ Realistic thermal spread was added to the first bunch to prevent its tail from collapsing under its own 
focusing force. 
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Figure 5.4: Simulation of 4 drive bunches in the linear reginne at the entrance of the plasnna. Top left 
panel: On-axis electron beann density (green area), longitudinal on-axis electric field (blue line), on-axis 
beann density (black line) and on-axis plasnna density (gray line). Top right panel: 2D dependence of the 
electron beann density (in units of n^ ). The white line shows the on-axis longitudinal position of the drive 
bunches, and the gray line (on the left) the sunnnned transverse profile. Bottonn left panel: Beann phase 
space p-z (in units of mc). The blue line shows the longitudinal position of the bunches. Bottonn right 
panel: The 2D dependence of the focusing force of the plasnna. 



The results of the simulation at the entrance of the plasma are shown in Figure 5.4. 
Using now the theoretical background that was detailed in the previous chapters of the 
dissertation we can make estimations for the system and compare them with the simulations. 
The theoretically predicted transformer ratio is 8 (from equation (60)), the driving efficiency is 
^ = ^11 -^1=64% -55% = 35% (from equation (28) and Figure 2.3) and the wakefield amplitude 

after the bunches is equal to 2.8 GV/m (from equation (59)). In the simulation the transformer 
ratio is 7.2, and the on-axis wake amplitude 2.9 GV/m. We also observe that the focusing force 
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Figure 5.5: Simulation of 4 drive bunches in the linear reginne after a 43 cnn long plasnna. Top left panel: 
On-axis electron beann density (green area), longitudinal on-axis electric field (blue line), on-axis beann 
density (black line) and on-axis plasnna density (gray line). Top right panel: 2D dependence of the electron 
beann density (in units of n^ ). The white line shows the on-axis longitudinal position of the drive bunches, 
and the gray line (on the left) the sunnnned transverse profile. Bottonn left panel: Beann phase space p-z (in 
units of MC). The blue line shows the longitudinal position of the bunches. Bottonn right panel: The 2D 
dependence of the focusing force of the plasnna. 



(in the bottom left panel) increases with each bunch and also varies within each of them, as 
predicted by equation (61). Overall the agreement between theory and simulation is very good. 
The simulation results after 43 cm of plasma are shown in Figure 5.5. The wakefield and 
the transformer ratio oscillate during the length of the plasma but on average they remain 
stable, probably a result of the chosen emmitances for the bunches. In the snapshot shown here 
the transformer ratio is 5.3 and the wakefield amplitude is 3.5 GV/m, since the bunches have 
slightly focused (especially near the tails). The most intriguing result of course is shown in the 
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bottom left panel which depicts the energy of the bunches, where we observe that the witness 
bunch has gained 3GeV / mxO.5m-l.5GeV , while the drive bunches have lost only 
1.5/ 6 - 0.25GeV each, thus demonstrating the high-transformer ratio scheme. 

The beam loading was kept light in this simulation in order to preserve the energy 
spread of the witness bunch. For the given witness bunch parameters, we can calculate the ratio 

E, n,,-R,„, 0.01-8 
Here E^ and E^ are the longitudinal components of the wake amplitudes left behind by the 
witness bunch and drive bunches alone, respectively (defined in equations (40) and (41)), and 
r^ =2sin(/:^w^ /2j & R^^^=2'M are their corresponding transformer ratios. Then, using 

Figure 2.7 we can estimate the total beam loading efficiency to be 77^ - 10% (since a^/ a^-O.! 
). The longitudinal energy spread is estimated through the ratio 

£^(0) n„R,^-lt,(0) 0.018 0.8 
The parameters shown here are defined in section 2.5.1. From Figure 2.5, the 3-o longitudinal 
energy spread is estimated f^\ -25% and the transverse energy spread (from Figure 2.8) is 
negligible since the witness bunch was made 10 times narrower than the drive bunches. 

5-3 Multiple Drive Bunches in the Blowout Regime 

In this section we will explore the possibility of using multiple drive bunches to enhance the 
transformer ratio and the efficiency in the highly nonlinear regime, where the bunches are 
strong enough to blowout all the electrons near them. Attempting to extend the techniques 
used in the linear regime, we place each subsequent drive bunch in the accelerating phase of 
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Figure 5.6: Sinnulation of 3 drive bunches in the blowout reginne. Top left panel: On-axis electron beann 
density (green area), longitudinal on-axis electric field (blue line), on-axis beann density (black line) and on- 
axis plasnna density (gray line). Top right panel: 2D dependence of the longitudinal electron field. The 
white line shows the longitudinal position of the drive bunches. Bottonn left panel: Plasnna density real 
space (in units of n^ ). The white line shows the longitudinal position of the bunches, and the gray line (on 
the left) shows the transverse profile of the bunches. The red line shows the plasnna density on-axis. 
Bottonn right panel: The 2D dependence of the focusing force of the plasnna. 



the Wakefield driven by the previous bunch, while at the same time ramping its charge so that 
the total Wakefield in the center of the bunch is the same for all bunches, and the bunches 
deposit their energy at the same rate. 

This scheme is examined through OSIRIS simulations only, since it is an analytically 
intractable scenario. The plasma density is set to n^ =2.0x10^^ cm'^. There are 3 drive bunches 
that have bi-Gaussian shapes with 6.7 jU/r? spot sizes (in either dimension), and their charge 
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scales as 135:405:675 pC. This corresponds to normalized beam densities of 
^bo / ^0 = 2.5 : 7.5 : 12.5 . The second bunch is placed 45 jU/r? after the first, and the third bunch is 
placed 55 jU/r? after the second bunch. This non-equidistant placement was found empirically 
and it is necessary because of the frequency shift that increases with higher drive charges in the 
nonlinear regime, which results from the electrons being blown out at larger radii from the axis, 
thus increasing the wavelength of the oscillation. The results from the simulation after 2 cm of 
propagation in the plasma are shown in Figure 5.6. 

We observe (top left panel) that although the on-axis wakefield at the center of each 
bunch is the same (blue line), it varies greatly along the length of each bunch. This is because in 
the nonlinear regime, the wakefield from highly-accelerating becomes highly-decelerating inside 
each bunch, resulting in the front of each bunch gaining energy and the tail of each bunch 
depositing energy. This sharp dichotomy is clearly observed when looking how the 2D r-z 
dependence of the longitudinal wakefield (top right panel) varies inside each bunch, thus 
reducing the driving efficiency and the transformer ratio which is only 1.5 in this example. This is 
in stark contrast with the linear regime where the wakefield inside each bunch is always linear, 
i.e. a section of a sinusoidal function, which at worse has a slope of k^ \ 

The bottom left panel of Figure 5.6 shows the plasma density real space r-z. The 
blowout region is clearly visible, and it is interesting to notice the electrons that are expelled 
away from the first bunch, which are then returning to the axis but they are blown out again at a 
greater radius by the second bunch. When the ions still attract them back, they are blown away 
for a third time by the third (and strongest) bunch. Only after that last bunch do the electrons 
have a chance to return back to the axis creating a density spike and a peak electric field of 



^ This occurs near the zero-crossing of the wakefield. 
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30 GV/m. Finally, the bottom right panel shows the 2D transverse focusing wakefield along the 
bunch dimensions. 

We conclude that because of the nonlinearities in the blowout regime the wakefield 
inside each drive bunch cannot be made flat and thus the efficiency and the transformer ratio 
are not enhanced in the same way they do in the linear regime using simply unshaped bunches. 
The most probable workaround to resolve this is to utilize shaped bunches, such as a single long 
triangular bunch or multiple properly shaped drivers. Those drivers will most likely require a 
delta-function-like front in order to bring the wakefield from accelerating to decelerating inside 
the rest of the bunch. However, for regularly shaped bunches (e.g. Gaussians, flat-tops) the 
linear regime is more favorable. Finally, a combination of the two regimes where the driving of 
the wakefield occurs in the linear regime and the beam loading in the mildly nonlinear regime 
may also be optimal. 

5-4 Chapter Conclusions 

The following items are concluded from this chapter: 

• By placing multiple drive bunches separated by an integer number of plasma 
wavelengths apart the wakefields add resonantly and the amplitude is maximized and 
increases linearly with the total charge in the beam. The transformer ratio and efficiency 
are minimized. 

• By placing multiple drive bunches at the accelerating phase of the wakefield driven by 
preceding bunches, then by adjusting the charge in each bunch the wakefield under 
each bunch can be made equal, thus increasing the transformer ratio beyond two which 
now scales linearly with the number of bunches. Single bunch efficiencies are obtained. 
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By fine-tuning the position and charge in each bunch the system can be more efficient 
than a single bunch system without shaping the bunches, at a slight expense of the 
transformer ratio. An example with efficiency 85% and transformer ratio 5.1 is 
presented using 4 drive bunches. 

By linearly increasing the emmitance on each subsequent bunch, 4 bunches are tuned to 
balance the transverse focusing force of the plasma and propagate over 43 cm of 
2x10^^ cm"^ density plasma while maintaining a transformer ratio of 6 in a 2D OSIRIS 
simulation. 

These schemes cannot simply be extended into the nonlinear regime because the wakes 
change very abruptly and each drive bunch experiences both energy loss and energy 
gain. This may be circumvented by shaping the bunches. 
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In theory, there is no difference between theory and practice; 
in practice, there is. 

- Chucl< Reid 



Chapter 6 

Multibunch Experiments 

The use of multiple drive electron bunches to excite plasma waves has been limited to only two 
notable sets of experiments, mainly due to the unavailability of multiple ultra-short drivers. 
Berezin et al. in the 1970s have reported the use of a series of 2 MeV electron bunches in 
10^^ cm'^ density 1 m long plasma and observed gradients of 0.25 MV/m [17, 176]. Similar 
experiments were conducted in the KEK 500 MeV linac in Japan, where 6 drive bunches were 
fed into aim long 5x10^^ cm'^ plasma and up to 20 MV/m gradients were recorded [126, 132, 
133]. In both cases there was very poor agreement between the results and the theoretical 
predictions. The bunches were separated by the linac RF electric field, which corresponds to 
several hundred picoseconds period. The goal of those experiments was to excite resonantly 
wakefields in the plasma. Such strategies can have interestingly high gradients only when high- 
density plasmas are used, which is not possible using the above mentioned technique to 
generate the bunches. 

In this chapter, we will investigate two methods to excite large wakefield resonantly. 
The first uses ~100 ultra-short bunches in 10^^ cm'^ plasmas created through the use of the 
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Inverse Free Electron Laser (IFEL) bunching technique. This method will be analyzed in section 
6.1. The second method for creating subpicosecond bunches uses a mask technique to 
selectively block portions of the initial beam profile. This method offers the advantage of 
creating adjustable charge per bunch and adjustable spacing between the bunches (along with 
the creation of a witness bunch), thus bypassing the IFEL method limit which is restricted to 
bunch separations fixed to the laser wavelength only. Beam diagnostics and beam-plasma 
interactions results using this method will be presented in sections 6.2 and 6.3, respectively. 

6-1 Resonant Wakefield Excitation with IFEL 

The experimental setups described in this section explore the possibility of creating GV/m 
wakefields by modulating a single electron beam through IFEL and then feeding the created 
bunches into a high-density plasma [88]. The motivation and beam-plasma simulations along 
with diagnostics of the microbunched beam will be presented here. The experiments took place 
at the Accelerator Test Facility (ATF) of the Brookhaven National Lab. 

6.1.1 Motivation and Theoretical Background 

Given an electron beam with a fixed charge and transverse dimensions, what is the optimal 
beam distribution that will excite the strongest wakefield in the plasma? Since the plasma 
operates as a simple harmonic oscillator in the linear regime, the wakefield amplitude will be 
maximized when the plasma is driven to its resonance by a periodic force with a frequency equal 
to the natural oscillation frequency of the system, i.e. the plasma frequency. This suggests that 
the use of a series of microbunches separated by a distance equal to the plasma wavelength X^ 

to be driven into such a plasma. For M Gaussian drivers, the beam density would be then given 
by 
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Figure 6.1: The plasma wakefield amplitude response of a non-bunched beam (top) and of a 
corresponding microbunched beam (bottom) of equal charge. The bunch separation is 10.6 |im and the 
microbunch length 1 pim FWHM. 
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Here A. is the normalized amplitude of each microbunch and ^j = j^p for resonant interaction. 

The response of the plasma to such a driver is shown in Figure 6.1, an example modeled after 
the beam parameters at ATF. 

The initial broad Gaussian-like beam distribution is shown on the top left panel, along 
with its plasma response in the top right panel for various plasma densities. For a 1x10^"^ cm'^ 
peak beam density, the plasma wakefield excited is around 100 MV/m. In the bottom left panel, 
the same total charge of the beam is distributed into a series of microbunches, separated by a 
distance 10.6 jum. The plasma response on the bottom right panel shows some new interesting 
features. Although the beam envelope response remains at the low densities (since the plasma 
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wavelength is too large to resolve the microbunches), when the plasma density is such that the 
plasma wavelength matches the separation of the bunches (here this happens at a 1x10^^ cm'^ 
density) then each bunch creates a wakefield which adds constructively and drives the total 
Wakefield in a resonant manner. The amplitude is now almost an order of magnitude larger, 
more than 800 MV/m, limited by the finite size of the microbunches and the residual charge 
between the peaks. 

Also visible are the harmonics of the interaction, appearing at densities 4x10^^ cm ^ 
9x10^^ cm"^ etc, where the plasma wavelength is an integer multiple of the resonant wavelength 
and there is now a microbunch placed every two, three or more periods of the plasma wave. 
The decreasing amplitudes of those resonances reflect the finite width of the microbunches and 
are determined by the Fourier transform of the microbunch Gaussian-like shape (a series of 
delta functions would create harmonic responses of equal amplitude). 

6.1.2 Experimental Setup and Simulations 

The idea of resonantly generating an enhanced wakefield is sought after experimentally at ATF\ 
There are two main parts in achieving the goal, 1) create equally spaced microbunches and 2) 
generate the corresponding resonant plasma density. A proposed experimental setup is shown 
in Figure 6.2. 

The process for creating the microbunches is the Inverse Free Electron Laser (IFEL) , 
whereby an unmodulated electron beam and a long laser pulse are fed co-propagating into a 
magnetic wiggler [27, 101, 137]. The longitudinally sinusoidal-dependent static magnetic field 
inside the wiggler causes a perpendicular force on the beam electrons and initiates electron 
motion in the transverse plane. Using this transverse motion the electron beam can be coupled 



^ This idea was initially suggested by Nan Ben-Zvi. 
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Figure 6.2: Experinnental setup of resonant nnultibunch experinnent at ATF. The unnnodulated electron 
beann is inserted co-propagating with a CO2 laser in the IFEL wiggler nnicrobuncher. The output nnodulated 
beann is then fed into a high-density plasnna, and the energy change innparted onto the beann is innaged 
with a spectronneter on a phosphor screen. 



to the laser pulse (by tuning the wiggler field period) such that some electrons feel the 
decelerating electric field of the laser and some other the accelerating field [38]. This picture of 
half-acceleration and half-deceleration repeats over every laser period along the beam length 
and results in a velocity modulation of the beam electrons [162]. After the exit of the wiggler the 
beam is left to propagate in vacuum over some specific distance, after which the low-velocity 
electrons have caught up with the high-velocity electrons thus creating microbunches, 
separated at the laser's wavelength [111]. In this experimental setup a 10.6 [um CO2 laser at 
50 MW power is used [96, 97, 163]. 

The required optimal bunching distance s^ for an electron beam with a given / is 
determined by the condition that electrons compress but never outrun each other, i.e. their 
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Figure 6.3: Dependence of the relative resonant wakefield annplitude on the nnicrobunch longitudinal spot 
size for a plasnna with a 10.6 [xm wavelength. 



relative drift distance within a laser period A^ is less than a quarter of the plasma wavelength. 
This is described by the formula 
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Here q is the fractional initial energy modulation imparted on the electrons at the exit of the 
wiggler, determined by the laser electric field. The approximate expression above holds true for 
small modulations and large /'s. For the 45 MeV ATF electron beam (/-lOO) the required 
drift distance for q = l% modulation is estimated at SQ-2.7m for the given laser parameters. 

As we will present in the next section, such microbunches have been successfully 
created and measured at ATF. The relative dependence of generated plasma wakefield 
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Solution of the 1 D nonlinear plasma equation 
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Figure 6.4: Simulated plasma wakefield amplitude response on-axis after the beam has fully entered the 
plasma. 1) Multibunched beam density 2) ThemOsiris calculated wakefield 3) OSIRIS calculated wakefield. 
For a transverse spot size 35 jum (shown here) and 500 pC of total charge, the generated wakefield 
amplitude is 7 GV/m. The plasma density is 1x10 cm' . 



amplitude on the Gaussian microbunch quality (i.e., their width) is shown in Figure 6.3 for a 
10.6 jU/T7 plasma wavelength and bunch separation distance. 

The simulated wakefield response on-axis for the multibunched ATF beam is shown in 
Figure 6.4 when the plasma density is exactly on resonance, 1.0x10^^ cm'^. The wakefield 
amplitude versus time is calculated using the fully explicit code OSIRIS, as well as the much 
faster quasi-2D code ThemOsiris, and they are shown to agree well. The wakefield amplitude 
resonantly excited from the 500 pC 35 jum spot size beam is 7 GV/m in this scenario. Since the 
individual wake from each bunch is added constructively, only the total charge that is fed into 
the plasma is important and not the charge distribution per bunch. Therefore, the wakefield 
amplitude at a given time is proportional to the integral of the charge up to that time. 
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Figure 6.5: Energy spectra of the 45 MeV microbunched beam for a plasma density exactly on resonance 
(left) and 3% off resonance (right). The initial energy spread is 0.1 MeV. The bottom of the figure shows 
the relative phasing between the bunches and the wakefield near the tail of the beam (the charge is too 
low to enhance the wakefield visibly). 

It is interesting to look at the electron energy spectra after a 15 mm long plasma 
interaction, both on and off resonance. This is shown in Figure 6.5, for a case exactly on 
resonance and for a case with a plasma density 3% off resonance, along with the relative phase 
of the wakefield with respect to the beam microbunches. 

We observe that the energy shifts drop sharply even slightly off resonance, as is 
apparent also from Figure 6.1: the resonance width is inversely proportional to the number of 
microbunches. For M =100 or more bunches seen here, the resonance width is on the order of 
l/M-1%, and the higher wakefield is obtained at the cost of a narrower resonance. In 
addition, in Figure 6.5 (left panel) we observe that the accelerating wakefield occurs only in 
between the bunches, thus the accelerated electrons are the residual particles outside the 
microbunch peaks (roughly 10% of the total charge using this IFEL scheme). In the slightly off 
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resonance case though (right panel), the energy gain is much larger and comparable to the 
energy loss since the wakefield now partially cancels out in each bunch, resulting in larger 
number of accelerating particles compared to the on-resonance case. This may be a welcome 
tradeoff for the experimenter. 

Ultimately, the crucial item in such an experiment is to obtain a reliable high plasma 
density with very high accuracy (1% or better). This is not straightforward since in typical plasma 
sources such as capillary discharges (see Chapter 3) the density varies by several orders of 
magnitude in microsecond timescales. 

6.1.3 Diagnostics of the Multibunched Electron Beam 

In order to diagnose the quality of the microbunching of the IFEL modulated electron beam, we 
collected the CTR emitted when the microbunches traverse two 1 [um thick titanium foils [177]. 
CTR is emitted because the electric fields of the electrons in the beam displace violently the free 
electrons in the metal surfaces, which in turn radiate due to the acceleration they suffer. 

In the experiment, the first metal foil was placed perpendicular to the beam direction of 
propagation and served the purpose of blocking the CO2 radiation used to drive the IFEL, since it 
could interfere with the CTR signal measurement. The second foil was placed at 45° with 
respect to the direction of propagation and emitted radiation through a window out of the 
beamline. The sum of the radiation emitted from both foils was collected. 

In contrast to the CTR interferometry technique used in the experiment described in 
Chapter 4 (see section 4.2.1 and Appendix C.4), the microbunch spacing here was too short to 
be resolvable using interferometry. The spectrum of the CTR radiation was recorded instead, 
since it contains information about the longitudinal current profile of the electron beam. 
Specifically, the on-axis spectrum is proportional to the amplitude-squared of the Fourier 
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Figure 6.6: The spectrum of the CTR emitted when the microbunched electron beam passes through a 
metal foil. 



transform of the current profile of the beam (see Appendix C.3). For a beam density defined in 
equation (66), the bunch form factor (which determines the CTR spectrum) is given by the 
formula 



ah^ sin^ -ah^ cos^ 



f{k;0) = e-^^''"''e 



j=i 



Ae 



+ik cos 0(^ J 



(68) 



Here k = 2n I X is the wavenumber of the radiation, is the radiation angle with respect to 
the radiation propagation direction \ a^ is the radial spot size of the beam, cj^ is the 

longitudinal spot size of each microbunch, M is the total number of microbunches, A^ are 
their relative charge amplitudes, and Zj = jX^ are their longitudinal locations. 



^ The detector has a surface area large enough to capture the radiation emitted at most of the angles. 
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Figure 6.7: Ratio of energy collected in each harnnonic of the CTR spectrunn for various nnicrobunch 
widths. The solid colored curves correspond to the theoretical predictions. The arrows on the left depict 
the experinnentally recorded data. The three sets of data agree at the region around a^ -0.1 /um. These 
ratios include the response of the detector and the transnnission of the filters. 



The geometry of the electron beam is uniquely mapped into its CTR spectrum, which is 
shown on-axis in Figure 6.6 (for the beam profile shown in Figure 6.4). The low frequency (long 
wavelength) radiation on the right corresponds to the envelope of the beam. For other 
wavelengths, the radiation in general adds out-of-phase except at the harmonics of the 
separation wavelength X^ =10.6//m between the microbunches, X^=\l m, m = l,2,... . The 
existence of radiation at each harmonic is therefore the result of the periodicity of the 
microbunch train. Finally and most importantly, the amplitudes of these harmonics are 
modulated by the Fourier transform of each microbunch, which is assumed to be a Gaussian 
profile. Therefore, the measurement of the amplitude ratio between these harmonics provides 
a direct way of estimating the microbunches' width. 
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In order to collect the CTR, a helium-cooled bolometer detector sensitive between 3 jum 
- 20 lum wavelengths is used. Furthermore, narrow Gaussian band-pass filters with roughly 0.5 
lum FWHM transmission were used to isolate the radiation of the first, second and third laser 
frequency harmonics, in successive measurements. After recording about 100 events at a 
roughly constant IFEL laser pulse energy (within a factor of 2) at all three harmonics, we 
calculated the three possible ratios of the signals. The data range for each ratio is shown on the 
vertical axis in Figure 6.7. 

Using equation (68) for small angles (^-0) and after accounting for the response of 
the detector and the filters, we calculated the expected theoretical ratio between the energy 
radiated under each of the harmonics as a function of the microbunches' widths. These three 
predicted ratios are also plotted in Figure 6.7. The data show that the measured ratios indicate a 
microbunch width around a^ =0.7±0.1//m, in good agreement with the predicted value of 0.5 
[um. Although the microbunch width was inferred using CTR before [111], this is its first direct 
measurement that also utilizes the information from different harmonics. 

In order to confirm the 10.6 jum separation between the microbunches, the IFEL drive 
laser was tuned to 10.2 jum, while still using the same narrowband 10.6 jum filter to detect the 
CTR. In that case the signal recorded was at least 100 times less than when the IFEL was driven 
at 10.6 jum and very close to the noise level of the detector, thus confirming the periodicity of 
the bunching at the laser wavelength only. 

6-2 Generation of Microbunch Trains by Masking Chirped Beams 

In this section we will describe a method for creating subpicosecond microbunches with 
adjustable spacing, which are required for most of the experiments described in this chapter. 
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Previous multibunch PWFA experiments relied on the different bunches to be created by very 
long beams being placed and bunched over different RF buckets in the linac [17, 79, 135]. Since 
those fields are in the microwave range, very low plasma densities had to be used (10^^ cm'^ - 
10^^ cm'^) in order for the plasma wavelength to be comparable with the bunch spacing, thus 
reducing the accelerating gradient significantly (which depends on the square root of the plasma 
density - see equation (2)). A future PWFA collider will require GeV-scale wakefields, which can 
only be achieved at high densities (10^^ cm'^ and above). The plasma wavelength at those 
densities is on the order of 1 ps or less, and as a result a multibunch PWFA accelerator requires 
bunches separated at subpicosecond-scale spacing, a challenging process [155]. 

One method to create such microbunches is outlined in section 6.1.2 using velocity 
bunching after an IFEL interaction. However, in this case the spacing is not adjustable because it 
is fixed by the laser's wavelength. The method illustrated in the following section 6.2.1 uses a 
metal mask in a high-dispersion plane of the accelerator to create microbunches and does away 
of those limitations. CTR diagnostics of the output microbunches will be then presented from 
proof of principle experiments in section 6.2.2. 

6.2.1 Mask Method 

The principle of generating a train of microbunches is shown schematically in Figure 6.8. The 
electron beam enters the "dog-leg" portion of the beamline having a correlated energy spread 
along its length. In the case shown here, the front of the beam ("F") has a lower energy than the 
back of the beam ("B"). The first dipole disperses the different energies of the beam in space, 
and thus the energy-time correlation becomes also energy-space correlation. Subsequently, the 
metal mask that consists of a periodic set of metal wires is placed along the path of the beam, 
thus spoiling the emmitance of the particles that hit upon the solid parts of the mask which are 
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Figure 6.8: Left panel: Simplified schematic of the microbunch generation technique using the mask. The 
beam with a correlated energy spread enters the dog-leg on the left, is dispersed in space, goes through 
the mask, and then is brought back to energy-time correlation. The letters "F' and "B" stand for the front 
and back of the beam, respectively. Right panel: microbunch structure created by the wire mask in both 
time and energy. Source: P. Muggli. 



then lost along the beam transport line. The energy-space microbunches then become energy- 
time microbunches when passing through a second dipole that removes the energy-space 
correlation of the beam. 

The advantage of this technique is that the mask is fixed in space and therefore the 
microbunches are extremely stable, affected only by the incoming beam variations. The initial 
correlated energy spread can be introduced by accelerating the beam off the crest of the linac 
RF wave. The number of bunches can be selected by tuning the width of a limiting energy slit 
aperture, thus selectively blocking some parts of the initial beam. The width of the bunches can 
be selected by tuning the width of the wires on the mask. Note that a fraction d / D of the 
initial beam charge is lost after hitting the mask, where D is the separation between the wires 
on the mask and d their individual width. In addition, the shadow that is cast by the mask is 
crisp as long as the finite transverse size of the beam at the mask, o-^, which is determined by 
its emmitance s^ , is much smaller than the individual wire width: 
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<^.=.F^«d (69) 



Here j3^ is the beam beta function at the location of the mask, s^ is the normalized emmitance 
of the beam, and / is its relativistic parameter. 

The distance between the microbunches is selected through the fractional energy 
spread imparted on the beam, AE/ E^, with E^ being the mean beam energy, for a given 

dispersion of the beam at the mask ^j^^^ask ■ Specifically, the longitudinal length of the beam at 
the exit of the dog-leg is o-'^=o-^±R^^AE/ E^, where R^^ is the dog-leg longitudinal dispersion 
function and cr^ is the beam length before the dog-leg. The transverse size L^ of the beam at 

the mask is a combination of its intrinsic beam size cr^ = ^P^Sj^ I y and the dispersion width 

added by the magnets, (J^ =?]^^^j^AE / E^. The mask is placed at a high-dispersion and low j3^ 

plane along the beamline such that o"^ <^ cr^ and thus 

K=4<^1+<^1 - ^, (70) 

Therefore the number of microbunches is N^=g^I D and thus the bunch separation is given 
by the formula 
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Although the mask here is periodic and the bunches are equidistant, a non-periodic mask with 
properly designed solid parts can be used to generate bunches at any position and with different 
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Figure 6.9: Left panel: Microbunches generated with tine nnasl< nnetliod dispersed in energy, for two cases 
with different nunnber of bunches. Right panel: Auto-correlation tinne-integrated CTR traces as a function 
of the interferonneter arnns path length. N=3 bunches with 1.4% energy spread and N=4 bunches with 
3.4% energy spread are shown. Source: P. Muggli. 
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mount of charge in each bunch. This is ideal for the multibunch setups explored in this chapter, 
especially for the creation of out-of-phase witness bunches. 

6.2.2 Beam Diagnostics 

In order to experimentally verify the generation of microbunches, CTR interferometry was 
employed in proof-of-principle setups at ATF [124]. This is similar to the technique described in 
section 4.2.1 that was used to diagnose the bunches in the double-bunch experiment, and is 
analyzed more in Appendix C.4. The CTR signal that is emitted when the bunches pass though a 
copper mirror is collected and sent through an interferometer that auto-correlates the signal in 
time. This auto-correlation time-integrated trace is recorded for different path lengths inside the 
interferometer by a liquid helium cooled silicon bolometer detector. A train of N equidistant 
bunches is expected to yield a symmetric trace with 2N -h 1 peaks. 



115 



For the experimental setup at ATF, the mean beam energy was E^^S^MeV , the 
dispersion of the beam at the mask was Tl^^^ask =1-4^, the initial envelope bunch length was 
cr^ =1650 //m (or 5.5 ps full width at half maximum), and the dog-leg dispersion was 
R^^ =-4cm. The total beam charge was 500 pC , the typical energy spread on around 1% and 
the beam emmitance was s^ =2mm-mrad . The mask was made of circular steel wires of 
diameter d =500 jum and separation period of D = 1270 jum. For this beam cr^ -140 jum and 
condition (69) is moderately satisfied, while (70) is easily satisfied. 

The left panel on Figure 6.9 shows the energy distribution profile of the beam after 
dispersed at the spectrometer at the end of the beamline, for AE/ E^ =0.5% . By changing the 
width of the limiting slit either 6 or 7 bunches could be selected. The modulation of the number 
of particles per bunch reflects the initial profile of the incoming beam. 

The right panel on Figure 6.9 shows two CTR auto-correlation traces for two incoming 
bunch parameters with energy spreads 1.5% & 3.4% and A^ = 3 & N = 4 bunches, respectively. 
The number of peaks in the traces is 5 and 7, in agreement with the 2N -\-l expectation value. 
The distance between the bunches is measured from the distance between the peaks in the 
traces (shown in the inset panel), and corresponds to 434±22jum and 216±60jum , 
respectively. Note that because only a portion of the CTR spectrum could be collected (low 
frequencies), the auto-correlation traces are not ideal and they have values smaller than the 
average. However this does not affect the bunch separation measurement. The above measured 
bunch periods are correctly predicted by equation (71) considering that the full width of the 
beam is a few times the quoted cr^ value. 
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Figure 6.10: Electron energy spectrunn of the bunched beann with 3 drive bunches and one strong witness 
bunch (out of phase) as recorded at the energy spectronneter at the end of the beannline. This is a 
favorable setup for a PWFA experinnent. The horizontal axis is increasing energy to the right, while the 
vertical axis contains beann transverse profile infornnation. 

Finally, an example shot of the energy spectrum of the beam with 3 drive bunches and 
one witness bunch out of phase is shown in Figure 6.10. The rest of the bunches have been 
blocked by limiting energy slit. Notice that the separation of the drive bunches in energy /SE is 
fixed by the period of the wires D for a given beam energy E^ through AE^E^-D / rj^^^^^. 

6-3 Beam-Plasma Interactions Using the Mask 

A series of experiments at the ATF have the goal of demonstrating proof-of-principle PWFA 
accelerators using multiple drive bunches. We will present here some simulations and 
experimental results specific to the ATF parameters, using the bunch generation technique 
described in the previous section. In the set presented here the drive bunches and the plasma 
density are adjusted such that the accelerator works in the maximum wakefield setup (see 
section 5.1.1) since we wish to first detect the maximum possible effect of the plasma on the 
bunches, and thereafter increase the transformer ratio in a future second set of experiments. 
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Figure 6.11: Plasma wakefield amplitude response as a function of the plasma density for a bunched and 
non-bunched square beam. The bunches are separated by one plasma wavelength and their width is half 
a plasma wavelength. 



6.3.1 Predicted Response 

The predicted plasma response as a function of the plasma density is shown in Figure 6.11. In 
the no-mask case the single 300 pC 4 ps long ATF electron beam is fed into the plasma, and the 
maximum wakefield amplitude is 60 MV/m assuming cr^ =100/im . The envelope response 
appears in the 10^^ - 10^^ cm'^ density range, while the slight beating visible at higher densities 
relate to the fact that the model of the beam was a square bunch with sharp cutoffs at the 
edges. In the case where the wire mask was blocking sections of the beam to create the 
bunches, the total charge is 150 pC and is assumed to be distributed over 5 bunches which have 
width equal to 125 \xrr] (half the plasma wavelength). The bunches are separated by one plasma 
wavelength apart. It is observed that the wakefields are very low for most densities, expect for 
the resonant density at 1.8x10^^ cm'^ (where the wakes of all 5 bunches add in-phase) where the 
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Figure 6.12: Energy spectrum of 5x30 pC drive bunches separated by 250 [um in space and by 0.25 IVIeV in 
energy as a function of the plasnna density after 10 nnnn of plasnna propagation. In densities above and 
below resonance the bunches suffer snnall energy spread, while at the resonance of 1.8x10 cnn' the 
Wakefield is nnaxinnunn and the energy spread beconnes large. The last drive bunches sannples the highest 
Wakefield. 



wake is expected to peak at 140 MV/m. Weaker wakefield peaks also appear at the plasma 
frequency harmonics. It is important to notice that the resonance width is around 20%\ which is 
the accuracy with which the plasma density must be known in order to operate near it. 

Since the bunches are also spaced apart in energy, each bunch samples a different 
phase of the wakefield and thus it will shift in energy. Those expected shifts in energy as a 
function of plasma density are simulated in Figure 6.12 using ThemOsiris after 10 mm of 
propagation in the plasma. The bunches are assumed equidistant in energy between 60 MeV 
and 61.5 MeV, separated by 0.25 MeV and their individual energy spread is 0.05 MeV. The first 
drive bunch has the higher energy, and is only experiences its own wake, thus the energy shift is 
relatively small. Later bunches experience their own wakefield plus the wakefield created from 



^ Inversely proportional to the nunnber of bunches. 
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Figure 6.13: Recorded energy spectra of the bunched (right panel) and non-bunched (left panel) before 

and after a 6 nnnn long plasnna at a low density. 



all the previous bunches. At the resonant density of 1.8x10^^ cm"^ the wakefield is maximum and 
this shows up as a strong spreading in energy, since the bunches are not much shorter than the 
plasma wavelength and sample half of its phase. In higher densities the bunches suffer small 
energy spreads. 

6.3.2 Experimental Data 

The experimentally recorded energy spectra of 7 bunches fed into a 6 mm plasma produced by a 
25 kV capillary discharge at 150 Torr gas pressure (see details on the plasma source in Chapter 
3) are shown in Figure 6.13 at a low density of 7x10^^ cm ^ selected by delaying the arrival of 
the beam with respect to the discharge. On the left panel the effect of the plasma on the full 
300 pC - 400 pC beam was recorded, showing energy loss gradient of 35 MV/m compared to 
the plasma off case (as deduced from the energy loss for the given plasma length). In this setup 
the beam length is between 0.5 mm - 1 mm and has a transverse spot size close to 100 nm. 
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Figure 6.14: Simulated wakefield and energy spectra for the experinnental data of the beann plasnna 
interaction of Figure 6.13. The top panel shows the 7 bunches with 200 pC total charge fed into a 
7x10 cnn' density 6 nnnn long plasnna. The wakefield under each bunch deduced fronn the energy shifts is 
shown with red dots. The bottonn panel shows the sinnulated energy spectra before and after the plasnna. 



When the mask is inserted the microbunches are now fed into the plasma (right panel), 
with a measured period (through CTR diagnostics) between 100 lum - 200 lum and width equal 
to half their period. The first drive bunch has the highest energy, around 59.4 MeV. The total 
charge now is reduced to roughly 100 pC - 200 pC, which results in a smaller average energy 
loss of 22 MV/m. This is still envelope interaction at the low density, as demonstrated by the 
increasingly larger energy loss of the later drive bunches. The simulated excited wakefield and 
energy spectra for the bunches is shown in Figure 6.14 and compared with the experimentally 
recorded energy shifts, demonstrating good agreement with the data. 

Operation near the resonance is more challenging because the plasma density must be 
stable within -20% for the bunches shown here over the desired plasma length. Plasma 
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instabilities larger than these values, or non-uniform longitudinally plasmas will result in 
decreased wakefield. As was shown in Chapter 3, for capillary discharges the density may vary 
by 3 - 4 orders of magnitude over a few microseconds, and the beam needs to be timed to 
arrive within this 20% window around the resonance. If the plasma density varies as 

^p{t) = n^^^e ""^^ , and the enhanced wakefield is to be observed in a range of densities An 
around a resonant density Hq, then the time scale A^ over which the resonance will evolve is 
given by 



At = -ln 

a 



^, \An^ 

, 2fZo y 



An . 



\ An 

(72) 

a nr. 



''O 



For example, if there is a AnI n^ =20% wide resonance and the density drops exponentially 

with a coefficient a - 1.9 jUS , then A^ - IQns . This is the time over which the density must be 
scanned in order for the resonance to be detected. Finally, instabilities in the incoming beam 
energy distribution may also affect the location of the resonance as the bunch period may drift 
during a run (see equation (71)). 

6-4 Chapter Conclusions 

The following items are concluded from this chapter: 

• Using the IFEL technique the 5.5 ps long ATF beam can be modulated into -100 
microbunches separated by 10.6 jU/r?, which is predicted to resonantly excite a 1 GV/m 
wakefield when fed into a 10^^ cm'^ density plasma. 
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By recording the spectral energy under each harmonic of the CTR radiation when the 
beam traverses a metal foil the length of the individual microbunches was measured 
close to OJ ±0.1jum. 

A method for generating microbunch trains of tunable separation is presented by 
placing a metal wire mask at a high energy dispersion plane along the beam path and 
selectively blocking parts of the beam. 3-7 bunches with 200 |im - 400 |im separation 
were measured using CTR. 

Using the mask technique a 140 MV/m wakefield is predicted to be excited if the plasma 
density is matched to the bunch separation. This is 2 - 5 times stronger than any non- 
resonant wake, with or without the mask. The mask pattern can also be adjusted to 
create bunches for a maximum transformer ratio scheme. 

Wakefields around 20 MV/m where measured at low densities where the plasma 
responds to the envelope of the bunches, by observing the bunch energy shifts which 
increase with later bunches. 
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Write a wise saying and your name will live forever. 

-Anonymous 



Chapter 7 
Conclusions 

7,1 Summary of the Results 

The main goal of this dissertation was to explore the use of multiple bunches as a tool for 
enhancing the performance of plasma wakefield accelerators. This was done through a unique 
combination of theory, simulations and experiments. In Chapter 1 we offered an overview of the 
role of particle accelerators in modern physics and we explained the need for new technologies 
based on harnessing the fields in waves driven in plasmas. In Chapter 2 the main aspects of 
plasma wakefield theory were rigorously laid out, and some aspects relevant to the transformer 
ratio and efficiency (both driving and beam-loading) of plasma accelerators were examined in 
more detail. 

The need for plasma sources was addressed in Chapter 3. The main characteristics of 
capillary discharges were analyzed and a large volume of data was presented, based on many 
months of experimental work in ATF at BNL and in the Pulsed Power Lab at USC. We found that 
densities between 10^^ cm'^ - 10^^ cm'^ can be generated stably using capillary plasma sources, 
tunable through gas pressure, tube diameter, tube length and also by varying the time-delay of 
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the arrival of the electron beam during the discharge. The spectral diagnostic of Stark 
Broadening was found to be an excellent method for measuring the plasma density as long as 
the light spectra during the discharge are relatively clean such that the spectral lines can be 
identified. 

In Chapter 4 we presented experimental evidence for the acceleration of a trailing 
witness electron bunch in the wakefields driven by a preceding bunch in a plasma. For the first 
time the accelerating gradient for such a scenario was above 100 MV/m while the bunch 
remains shorter than the plasma wavelength, thus being a promising technology for accelerating 
bunches to high energies while preserving the beam quality. It is most likely that any future 
plasma-based collider will require the acceleration of a trailing bunch in a similar fashion. 

Chapter 5 laid the theoretical framework for creating high-transformer-ratio wakefields 
in a plasma accelerator through the use of multiple drive electron bunches. We find that it is 
possible to multiply the energy of an incoming electron beam when the position of the bunches 
and the charge per bunch are adjusted properly, as long as the beam can be maintained stable 
over long plasma lengths. The second half of this chapter demonstrates through simulations 
that indeed a witness beam can propagate over meter-scale plasmas if its emmitance of each 
bunch is properly matched to the increasing focusing force that each subsequent bunch 
experiences. Other than those elements the bunches need not be specially shaped, which 
implies that these ideas could be tested straightforwardly experimentally. 

Finally, in Chapter 6 we presented some recent experimental results and diagnostics 
related to multibunch experiments. Two different techniques for generating multiple bunches 
were examined, either through IFEL bunching or by using the mask technique to selectively 
block portions of the incoming beam. Although some charge is inadvertently lost with this 
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method, it is found to offer great design flexibility, such as the easy creation of a witness 
electron bunch. Preliminary experimental results of the interaction of multiple bunches with 
plasmas are also presented in this chapter. 

7.2 Future Work 

After the exciting possibilities that multiple bunches open up for plasma accelerators, a variety 
of approaches require experimental testing. The most immediate one would be the acceleration 
of a trailing witness bunch using SLAC-scale gradients (GV/m) to very high energies. This would 
be a direct extension of both the energy doubling experiment [20] (which had a single bunch 
and high energy spread) and the trailing bunch acceleration experiment presented here in 
Chapter 4 (where the gradient was large but the energy gain was limited by the plasma length). 
Furthermore, such experiments need to also be conducted for positrons, the equally important 
other half of a possible PWFA-based future collider. 

In addition, there is a lot of room for improvement in the demonstration of resonant 
plasma wakefields using multiple drive bunches. The wakefields need to be measured 
thoroughly as a function of the plasma density and the number of bunches in the drive beam. 
Subsequently, it is extremely important to demonstrate experimentally that multiple drive 
bunches can provide wakes with high transformer ratio. So far, only a small enhancement of 
30% has been verified in dielectric accelerators [79]. Using the mask technique presented here, 
the required bunches can be designed. Using the plasma source suggestions also presented 
here, the required plasma source could also be built. ATF could provide the first proof-of- 
principle results at low energies, but eventually these schemes need to be tested in GeV-scale 
machines. 
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Simulations can also help greatly in examining the possibility of using multiple drive 
bunches over very long plasmas. There is still the question whether plasma instabilities, such as 
the electron hosing instability, will hinder the acceleration of witness beams. The electron 
hosing instability arises when the beam accelerates in pure ion channels in the blowout regime 
of PWFA, triggered by small transverse motions of the tail of the beam. This may not be an issue 
if the accelerator operates in the linear regime (as examined in this work) since there is no ion 
channel, yet this needs to be verified. 

The work here can be extended in more detail in the nonlinear regime. Although using 
unshaped bunches does not seem to work directly in the nonlinear regime, simulations can 
indicate whether properly shaped bunches (such as a step function followed by a linear ramp) 
can lead to transformer ratio enhancements in the blowout regime as well. Such bunches, 
although more technically challenging to create, would take advantage of the nice focusing 
properties of the blowout regime and from the higher acceleration gradients (thus reducing the 
required plasma length, if also the high transformer ratio can be maintained). 

Once high beam loading efficiency and high transformer ratio beams are generated, it 
will be important to generate them at high quality, starting with low energy spread. Currently 
this seems to be feasible by shaping the witness bunch in a reverse-ramped fashion that causes 
a constant wakefield under the bunch, thus reducing the energy spread. Another idea that has 
not been explored yet is to achieve that using the technique that is popular with conventional 
accelerators, which is to accelerate the beam on different slopes of the wakefields. By tuning 
the phase the beam can sample either a negative slope or a positive slope of the accelerating 
wave, thus on average compensating for the field variation across the beam. A similar effect 
could also be achieved in a plasma wakefield accelerator by tuning the phase of the witness 
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bunch. For example, one could implement a transverse magnetic field so that the beam acquires 
a transverse velocity that results in its dephasing with respect to the wake - this was suggested 
initially in [91] to actually avoid dephasing. Alternatively, plasma sections with different plasma 
densities could be placed in an alternating fashion such that the plasma wavelength is different 
between sections and thus the beam could sample in a similarly alternating fashion two 
different crests of the plasma wake, averaging out the wakefield variation. 

7.3 Epilogue 

We have come a long way since the ancient Greek philosopher Thales of Miletus wrote that by 
rubbing animal fur onto amber^a particular attraction was observed, thus providing the very 
first understanding of the electromagnetic forces. Had he been around Geneva today to see the 
underground tunnel of the Large Hadron Collider, without a doubt he would have been awed by 
the intricate ways that we have devised to manipulate these forces in order to probe the 
fundamental mysteries of Nature. Though the means may have changed, the spirit for discovery 
has remained unaltered through the ages; it has also inspired the writing of this dissertation in 
an attempt to allow plasma accelerators to aid the future of particle acceleration. 

This dissertation will end in the same way it started, with a quote from Richard 
Feynman, this time taken from his "Lectures in Physics" to the Caltech undergraduate students 
in 1962. It is probably the shortest piece of text that I have ever seen to summarize all of the 
intricate interconnections of science so beautifully. 

"A poet once said: the whole universe is in a glass of wine. We will probably never know 
in what sense he said that, for poets do not write to be understood. But it is true that if we look 
in glass of wine closely enough we see the entire universe. There are the things of physics: the 



^ Amber is ''nAEKipov'' in Greek, hence the name of the electron. 
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twisting liquid which evaporates depending on the wind and weather, the reflections in the 
glass, and our imagination adds the atoms. The glass is a distillation of the earth's rocks, and in 
its composition we see the secrets of the universe's age, and the evolution of the stars. What 
strange arrays of chemicals are in the wine? How did they come to be? There are the ferments, 
the enzymes, the substrates, and the products. There in wine is found the great generalization: 
all life is fermentation. Nobody can discover the chemistry of wine without discovering the 
cause of much disease. How vivid is the claret, pressing its existence into the consciousness that 
watches it! And if in our small minds, for some convenience, divide this glass of wine, this 
universe, into parts - physics, biology, geology, astronomy, psychology, and so on - remember 
that nature does not know it! So let us put it all back together, not forgetting ultimately what it 
is for. Let us give one more final pleasure: drinl< it and forget it allV 
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Appendix A 

PWFA Linear Theory 

A.l Impulse Response of a Plasma 

In this section the electric field that is excited when a relativistic charged beam enters a plasma 
is derived. The analysis applies to linear theory, where the beam density n^, is much smaller that 

the plasma density n^, nf^^n^. Also, we will only deal with two dimensions in a system will 

d ^ 
cylindrical symmetry, so there is no azimuthal dependence and — = . 

d<l) 

We start by writing Maxwell's equations in a homogeneous linear medium: 

dB 

1 dE 
P 



VxE = - 

dt 

^^^^^^■'^^c^ar (73) 



V£ 



^0 



Tl M 

Since in this regime the ion motion is negligible, — ^ «: — - , we can assume static ions and only 



^p ^e 



write down Newton's second law of motion for the electrons: 

dp, 



= -eE-evxB (74) 

dt 



Finally the equation of continuity for the electrons is included: 



143 



^ + V.(«^vJ = (75) 

The electron momentum, the charge density and the current density are defined as follows for 
an electron beam: 

1 

P = Pheam + Pions + P electrons = "^^Z, + ^^i " ^^e (^6) 

^ ~ ^ beam "^ ^electrons ~ ~^^b^b ~ ^^e^e 

There are 10 unknowns in the above 12 equations: E,B,v^,n^ . There are 2 more equations than 
unl<nown variables because Gauss's laws for the electric and magnetic fields can be derived 
from the other two Maxwell's equations when they are combined with the continuity (or charge 
conservation) equation. 

The next step is to linearize the equations, assuming that the quantities that are 
examined are just small perturbations over the initial state values (at ^ = ). The unknowns can 
be written as follows: 



E 





+ E, 


B 


-5, 


+A 


t 





+ Vi 


«e 


-Mo 


+ «! 


«,. 


= «0 





(77) 




There are several items to note here. The plasma starts from a neutral state, so the initial 
density of electrons is equal to the initial density of ions. Also, since the plasma is assumed 
homogeneous there is no initial net electric or magnetic field. In addition, the plasma electrons 
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start out motionless (cold) so there is no initial velocity, and the ions are static so there is no 
density perturbation assigned to them. Lastly, the beam is assumed to move close to the speed 
to light in the z-direction, i.e. v^ - c£ . The equations now become in the nonrelativistic limit: 

' dt 

- 1 9£ 
c at 

V • £j = — (-^«(, + ^«o ~ ^% ~ ^"i ) 

^0 (78) 



m 






dt 



By ignoring the second order terms (the products of two first order terms) because these first 
order perturbations are assumed to be very small compared to the initial values (n^^n^, etc), 
the equations become 

' dt 

~ I dE 

c dt 

^0 ^0 (79) 

V-5i=0 

dt m^ 

dt 
The first goal now is just to solve for the density perturbation. Taking the time derivative of the 
continuity equation and then replacing from the equation of motion: 



145 



r^ + ^.V • ^- = ^ r^ ^ V 'K=0 80 



The electric field can be replaced from Gauss's law to yield: 



1 I ^^2 _ ^2 



con. = -con. (81) 



2 ' ^p'n ~ "^p'^b 



dt 



A,2 The Density Perturbation Response 

Let's now assume as an input beam a delta function, a single electron at a specific point R in 
time and space moving at a relativistic speed v^ - cz into the plasma: 

n,{r,za)=5(R-v,t)=5{7)5{zz-ctz) = 5{7)5{^) 

1 (82) 

since the beam is dependent not simply on the time variable t but on the variable z-ct , we 
expect all the parameters in the system to by finally dependent only on ^ = z-ct . This means 

d d 

that the time derivative — can be replaced with -c — . 

dt d^ 

Equation (81) now becomes 



de 



-kln,=-klS{7)S{^) (83) 



0) 

Where ^ =— ^ is the plasma wavenumber. In order to solve this differential equation the 



' c 



Laplace transform is applied (for the (^ variable): 

s^L{n,}-sn,(0)-n,'(0)+klL{n,}=-klS{7)^ 

, V 

k^ (84) 

5 +Kp 
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This is the density perturbation impulse response of the plasma. Here f/(^) is the Heaviside 
step function (zero for ^<0, unity for ^>0, and equal to 34 for ^ = 0). The total density 
perturbation for a general beam density fz^(r,^) = fz^^(r)fz^||(^) is found by integrating the 
above result (Green's function) over the whole space: 

+00 +00 

n,{r,^)--k^jd^'l2;rr'dr\,{r')n,,{^')S{fym(k^{^-^'))u{^-^')^ 

-00 

+00 +00 



ntl{r) 



n, (r,^) = -k^n,, (r) j n„ (^')sin (k^ (# - ^'))d^' = 

— oo 

= -^Z(#)n,,(r) 



(85) 



Where in the last step we have defined the wakefield function (for reasons that will become 
apparent later) as 



e '- 



^i^>^i''^i^i^>HM^-^')y^' (86) 



^OS -00 



A,3 The Electric Field Perturbation Response 

Now that the density perturbation is known, the response of the electric can be calculated. First, 
take the curl of Faraday's law: 

Vx(VxE^) = (VxBj) ^ V^£i-V(V-£i) = — (Vx5i) (87) 

Now insert the curl of the magnetic field from Faraday's law and the divof the electric field from 
Gauss's law: 
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V'^. -V 



e e 



V ^0 ^0 y 



d_ 
dt 



c dt J 



Sr, Sr, dt dt c dt 



(88) 



•V^E 



"0 ''0 



' C^ dt' 



As before, since everything is a function of ^ = z-c?, the derivatives with respect to time can 

d d 

be replaced as >-c — . Then 

dt d^ 



V'E 



d% 

' ' d4' 



S^.,..2..,dn, e^^^_^^^^ 



d£, 



d4 s, 



(89) 



*^0 



The only unknown here (other that the electric field) is the electron velocity perturbation v^ , 

dv. e - dv. e - 
which can be simply replaced from Newton's law of motion: — - = E^ or — - = E^ . Then 

dt m 5f cm 

e ^ e 



^ d'E, 



e - 



dm e - e - 



V^E^ Y^ecn^ju^ E^ ^ec^ Id^z — Vn^ Vn^ 



d^ 



cm 



d^ s. 






df " 



^dn, e 



E,^+ec fx^z—- Vn, Vn, 

d^ ^0 ^0 



(90) 



Now consider that the V operator can be written as 



" aV 


_1 d 
r dr 


< d\ 

r — 


d' 


^ .d .d 

dr dz 









(91) 



However the variable z is not independent but we have written it as a function of ^ . Therefore 

d d 
since E,=z-ct we can just replace > — . The operators become 

dz dE, 
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v'=v'+- 



1 d 



r ^\ 



r- 



d 



' d'e rdA drj d^e (92) 

V=r — +z — 
dr d£, 

Now the quantities on the right hand side become 

^"' ^ -% |:[<5(r )sin(^/)[/ (^)] - z\ A[^(F)sin(^/)c/ (^)] = 

Now substitute bacl< at the electric field wave equation and separate the two vector 
components, one for the radial direction and one for the longitudinal: 

(94) 

The last equation for the longitudinal electric field can be easily solved in the £, dimension but 
it also requires to be solved in the r dimension, which will yield the field off-axis. The Green's 
function response G{^\7^ to an impulse function 5{r^ in cylindrical coordinates is the 
solution to the equation 

(v^-/:;)G(r;F) = ^(r-F) (95) 

This solution can written as 

G{f-rr') = ~K,{k;\r-r"^ (96) 
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Kq is the modified Bessel function of zeroth order. The longitudinal electric field can be written 



as 



Eu {r,4) = +:^KKo {k^r)cos{k^^)u (^) (97) 

This is the longitudinal electric field response of the plasma to an impulse function. It can be 
divided into two components, a parallel one and a transverse one, namely E^^(^^) and R{r) 

such that ^1, (r,^) = E,, (^)/?(r) . This yields 



kl 



(98) 



R{r) = ^K,(k^r) 



It is interesting to observe that while the plasma density perturbation starts from and it is 
proportional to a sinusoid, the electric field perturbation is actually discontinuous at the starting 
point as it has a value of Vi of its maximum value and then it evolves as a cosine, being always 90 
degrees out of phase compared to the plasma density wave. 

Now assume a general beam density distribution of the type 

n,{r,(p,^) = n,^X^)n,^{r) (99) 

The perpendicular beam density is assumed unitless with peak amplitude of 1, and therefore the 
SI units of the longitudinal beam density are m'^. Then the total electric field is found simply by 
integrating Green's function over the whole beam distribution in space [92, 116]: 

(100) 
R{r,<p) =^ I dcp' J n,,{r')K,{k;\r-r'\)r'dr' 
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Here \r -r\ = Jr^ +r^ -2rrcos((p-(p') . Also notice that £"11 (^) = Z'(^) . For rotationally 
symmetric bunches where ^bi{^'^9')-^bi{^') ' by using the property that 



I /(cosx)Jx= I /(cosx)Jx, one can show that /?(r) reducesto 



1^2 271 +(X) 



R[r) =^ i d(p' j n,^{r')KAk^4r^ + r'^ -2rr'cos(p'ydr' 



(p'=0 r'=0 



■K ] n,,{r')l,(k/)K,(k^r)r'dr' + kl \ n,,{r')l,(k^r)K,(k/)r'dr' (101) 

r'=0 r'=r 

+CX) 



Here r^=min{r,r'} and r^ = max {r,r'} . 

A.4 Bi-Gaussian Bunches 

Since the Gaussian is the most common type of beam density shape, it shall be analyzed 

explicitly. For a transversely Gaussian bunch profile n^^{r) = e ^'^'^ , the on-axis perpendicular 
component of the longitudinal wakefield is 



R{^). 



klcTl[omi91-\n(k^cT^)\ 



f l2 2\ 



\ J 



e 2 r 






V J 



k^(T^ <0.5 



\/k^a^ 



^^o-,»l 



(102) 



+G0 -t +C0 -bt 

e , re 



Here T(0,b)= \ — dt= \ dt = E^{b). A longitudinally Gaussian bunch profile is of the form 



%\{^) = ^b^e 



la] 


_ ^ :k 




/O \3/2 2 



(103) 
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The total number of particles in this bunch is 



N^\n,{n^)dV=\n,,{r)27rrdr j n,,{^)d^ 



(104) 



Near the center of the Gaussian bunch the field needs to be evaluated numerically. Near the tail 
of the bunch though it can be calculated analytically: 



r \ 



^,(#»^j- 



v^^^oy 



'2N^ 



V^r / 



COS 



M 



(105) 



Or, in normalized units 



f5li^^=^ 



r ^ 



mccD„ 






k^a^e 2 



cos 



(^.^) 



(106) 



Then the total on-axis wakefield left behind a bi-Gaussian electron bunch is 



Or, in normalized units 



f \ 
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v4^^oy 



aiki kiai 
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2 



cos 
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(107) 



(108) 



The absolute wakefield amplitude is increased as cr, nI' and cr^ n' and N't (for a constant 
plasma density). Finally, let's investigate the dependence of the plasma wake amplitude on the 
plasma density, or equivalently the plasma wavenumber k^ . There are two questions to be 

answered. First, given all other parameters constant, what is the optimal plasma density? 
Second, what is the wake that can be obtained for that plasma density? 
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A,5 Optimal Plasma Density for bi-Gaussian Electron Drivers 

In the ID limit, one can straightforwardly find the plasma density for which (106) is maximized 
for a fixed density perturbation n^^ln^: it is k^G^=^2 . The maximum normalized wake 
amplitude then is given by 



eE-{^»CT^) _^^^ 



f \ 

riu, 



mcco„ 



(109) 



V ^0 y 



If the bunch shaped is fixed instead, equation (105) indicates that the absolute wakefield 
amplitude increases as the plasma density decreases, or k^ i . This will be true as long as the 

ID linear regime condition is satisfied, i.e. nj^Q^riQ and k^a^^l. Therefore, while in that 

regime, the optimum plasma density is the minimum. 

When the 2D effects are included, in order to get the plasma density that optimizes the 
wake amplitude for a fixed bi-Gaussian bunch, the derivative of the wake amplitude expression 
(107) with respect to k^ needs to be taken and set to zero: 



_d_ 
dk 



P iv^^^oy 



cr'kl kW. ri2 2\ 



■^z "p "p 



Nkle 
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■2E, 
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^K-l^ 



\ J 



+ (7'XE, 



^K-l^ 



\ J 



-cTlk% 



^K-l^ 



\ J 



= 



\ J 



In this last step we used the expression for the derivative of the exponential integral. 



d_ 
dx 



E^ (x) = -Eq (x) , where E^ (x) = I ^ ""^dt = — . The result can be written as 



k>l 



'K-l^ 



^K-l^ 



V y 



(110) 
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One would have to solve this nonlinear equation numerically in order to derive the optimum 
plasma density. However, it is found that in most cases the exact solution can be approximated 
simply by 

klcrl 



:1 ^ k^cj^^^ 



(111) 



This is the ID result. It is true for example for small aspect ratio bunches where 
k (j^ ^k (7^ ^>^-^l and the second term on the right hand side of (110) can be ignored. 

When this aspect ratio is larger than about 0.1, the exact optimal k^ differs from the simplified 
optimal kp by a large factor; however, the approximate maximum wake amplitude never differs 
more than 20% from the exact value even if the approximate expression k^ci^ - v2 is always 
used [116]. This is due to the slowly-varying dependence of the wake amplitude on k^cj^ . 

The Wakefield amplitude for the semi-optimal condition k^cr^ = v2 (using eq.(108)) is 



£--(r = 0,^»aJ-0.74 
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(112) 



In normalized units this reads 



eE£{r = 0^l^aJ^^^ 






mccD„ 



V ^0 y 



2 +^ ^ ^2^. 



^g "^r 



V ^. / 



(113) 



This indicates that other than the slowly-varying monotonically decreasing term in the brackets, 
the optimum absolute wakefield amplitude will scale as 1/cr, if the plasma density is adjusted 
accordingly for a given a . 
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The above results for bi-Gaussian drivers in the linear regime can be summarized as 



follows: 



eE^ 



mcco^ 
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(114) 



And in SI units: 
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(115) 



The convergence of the above three expressions are compared in the following Figure. 
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Figure A.l: Comparisons of the maximunn wakefield expressions as a function of the beann aspect ratio for 
bi-Gaussian bunches, when the plasnna density is adjusted such that k^a^ = ^|2 . 
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Appendix B 

Upper Limits to the Transformer 

Ratio 

B.l Symmetric Bunches 

It can be proved under very general assumptions that the upper limit to the transformer ratio 
for a symmetric bunch inside a single mode lossless medium is equal to 2 [10]. Specifically for 
the plasma wakefield case one can use the results for a linear wakefield excitation. Recollect 
that the longitudinal wakefield excited from a drive bunch distribution in the ID linear regime is 
(from eq.(lOO)) 

^0 -00 

The wakefield at a point ^ inside the bunch with distribution ^^(<^) is 

^0 -00 
^0 -00 ^0 -00 

Assuming a symmetric bunch around ^ = 0, at the center of the bunch the wakefield is 

E;{0) = -]n,,{^')cosk^^'d^' (116) 
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If the Wakefield at the center of the bunch is also the peak decelerating wakefield then 
^max ~^z{^)- 'f ^1^'^ '^ ^^^ ^^^ ^^^^' ^^^^ '^ ^'" ^^ larger than this value and Le^^^ > £^"(0) . 
However in any case Le^^^ > £"(0) . Also, the wakefield at a point c^ after the bunch is 
(assuming that the bunch distribution is localized around ^=0 without extending too far): 

+00 

^0 -00 

+00 +00 

= —cos/:/ j n,|| {^') cos k^^'d^' + —sin/:/ j n,,, (^')sin/:/'J^' 

^0 -00 ^0 -00 

Assuming now that the bunch distribution is symmetric (and the integral is done symmetrically 
around the bunch), then the second integral in the above expression is because an odd 
function is integrated symmetrically (the bunch is an even function multiplied by the sinusoidal 
which is an odd function, resulting in an odd function). Then the maximum value is reached 
when cos^^^ = -l, yielding a result for the maximum accelerating wakefield: 

+00 

e 



^0 -00 

= -2- j n„ {^')cosk^^'d^' = 2\e: (0)1 < 2\e-^ 



R = 



E^ 

max 



<2 (117) 



Therefore we conclude that if the peak of the decelerating wakefield happens to be at the 
center of a symmetric bunch then the transformer ratio peaks at 2; in any other case it will be 
less than 2. 
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B,2 Asymmetric Bunches 

The transformer ratio can be increased beyond 2 if a proper asymmetric bunch is utilized [9, 
104]. The maximum possible transformer ratio (along with the beam to plasma energy transfer 
efficiency) is achieved when the retarding wakefield inside the bunch is constant. The bunch 
distribution that satisfies this condition can be calculated analytically as follows. 



e ' 



+00 77- 

—00 

Here f/(^) is the Heaviside step function. It is assumed that there is a constant wakefield that 

starts at ^ = 0. The longitudinal beam distribution ^^||(^) can be solved with the use of the 

Laplace Transform. By taking the Laplace Transform of both sides and applying the property that 
the Laplace Transform of a convolution integral is just the product of the Laplace Transforms of 
the individual functions that are convolved, we get: 



JjnJr)cos(^^(#-|'))«(^-n< =n^^^(^) 



•pK,(^)} ' -'°^°' 
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5 + ^^ e s 

F F ( 1 ^ 



(119) 



So, mathematically the optimal bunch distribution is linear ramp with a delta function at the 
beginning (required to jump-start the wakefield that remains constant thereafter). The exact 
bunch distribution and retarding wake can be calculated if we specify a linear ramp with ISl^ 
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particles and length k L = M eZ, and a delta function with N^ particles. This can be written as 



(for a transversely Gaussian profile): 



^b{^^^)-- 
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(120) 



The Wakefield excited inside tiiis bunch is 



M^) = 



eN, 



s^TTcrlkll} 
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(121) 



This wal<e is constant if 



^ 
N. 



^ = 2;r'M' 



(122) 



The wal<e excited after the bunch is 



E,{^)- 



eN; 



^/l+(2^Mf 



V ^ 



cos(^^^ + tan"^(2;rM)) (123) 



Therefore, the optimum transformer ratio is then 



E I rM»i 

R = ^ = ^1 + (2;tM) -^ 2nM 



(124) 



This is the ultimate limit for the transformer ratio in the linear regime. 
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Appendix C 

Coherent Transition Radiation from 

an Electron Bunch 

In the following pages we will attempt to calculate theoretically the CTR from an electron bunch 
when going through a metal foil. There are two independent effects that take place in such a 
phenomenon. First, we will calculate the energy radiated from a single electron when it hits a 
metal foil (transition radiation). Second, we will calculate how this radiation is altered by the 
presence of the other electrons nearby (coherent radiation). A more complete description can 
be found in [110]. 

C-1 Transition Radiation from a Single Electron 

The problem of the radiation emitted from an electron approaching normally a perfectly 
conducting boundary is handled using the method of images: the perfect conductor is replaced 
with another electron of equal and opposite charge (i.e., a positron) approaching the first 
electron from the opposite direction with the same initial speed, such that the two particles are 
to collide exactly where the boundary was. A significant assumption here is that the two 
particles decelerate until they meet at the boundary with speed each. Transition radiation will 
be emitted from both electrons as they decelerate to reach standstill. 

The energy radiated from a single charged particle during a collision is given by the 
following formula [77]: 
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d^E 



1 



dcodO. \6n s,x 



•d_ 
dt 



qnx(nx/3] 
\-h-p 



'rico{t-h-r{t)lc) 



dt 



(125) 



The quantity on the left has units of energy per steradian per frequency, in CGS units, ji-v I c 

is the speed of the electron, h is the direction of radiation at frequency co , and r is the 
position of the electron at time t. Since the time T it takes for en electron to cross the 
boundary is small, we let cot <^\ which will set the exponent inside the integral to go to 0. This 
is valid if the radiation of interest is at relatively low frequencies. The integral then becomes 
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Here the positron radiation term was also added. Carrying out the calculation 
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0" sin' e 



4^V(l-/?'cos'^) 



, cos^^n-(-y0) (127) 



So this is the radiation from the transition of the electron approaching perpendicularly a perfect 
conductor. Note that the exact same expression holds true for the case of an electron emerging 
from a perfect conductor into the vacuum. Also, observe that due to the low-frequency 
assumption there is no dependence on frequency in this formula. The peak of the radiation is 
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, , -, located at an angle sin 6^ = — . More general expressions are available for general 
\6n £^13 c py 

cases of oblique incidence [110]. 
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C-2 Coherent Radiation from N electrons 

Let us assume N electrons, distributed in space. Each electron is identified by the number j , so 
j =l,2,...,N . Let the center of this electron bunch to be at the origin of the axes, that the 
observation point is located at R = RR, and that each electron j is located at a distance 
~00 ^^^^ ^^^ center of the bunch. Hence the distance between an electron j and the 
observation point is Xj =R-rj =Xjhj . 

Now, assume that all the electrons are moving at exactly the same speed. This motion can be 
broken down to a set of frequencies. The field generated by an electron j at the frequency co is 

E.ykj, where k = — . The field observed at the observation point which is Xj away will 



propagate with a wavevector kj = khj , so its value at that point will be 

E. {k)e'''''' = E. (^)/*'^<-'^"'' = E. (^)g-*"^"^/*'^"'-^ (128) 

Assuming the far-field approximation, i.e. 7?»r^,then fij-R and E.{k)-E^{k) , which 

means that the electric field amplitude is same as if the electron is located at the center of the 
bunch. The field now becomes 

£ {k)e-''^'^e'^'^^-^ = £ {k)e''''e-^''-'^ (129) 

The total field intensity Ifotdi^)^^ Watts per m^ generated from all the electrons at that 
observation point for that frequency k is 
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CS I ^ 1 2 

Here /^(/^) = — ^£'^(^) is the field intensity generated from a single electron. If all the 

electrons were located in the same point in space, then the phase information would disappear 
and this would yield a factor N^ (it's just like having a charge e-N there, which would create 
an electric field N times stronger and hence intensity N^ times stronger). However the 
distribution of the electrons along the bunch will alterthis simple dependence. 



j=l m=\ 
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(131) 



In order to proceed further, one more assumption needs to be made. Since the number of 
electrons in each bunch is typically large, we could replace the summations over all the 
electrons with integrals over a continuous electron distribution. So, let's assume that the 
number of electrons in a volume element d^r at a distance r from the center of the bunch is 
N'S{r)d^r . Then the function S(r) is just a probability distribution such that \s(r)d^r = l. 

Therefore, the discrete summation can be approximated by the integral (after some short 
calculations) 



t,o.Ak)-h{k) 
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(132) 



f(k;R)^\fe^"^'S{r')d'r' 



f is traditionally noted as the bunch form factor, a quantity that relates the bunch distribution 
with the phases of the electrons inside the bunch. For k^O (low frequencies, high 
wavelengths) the factor is close to its maximum value, 1. In this case If^f^i[k^O)^N^I^(^k). 
Here the electrons act in harmony, and they produce the same intensity as if they all were 
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perfectly located at the center of the bunch. Hence in low frequencies the actual distribution of 
electrons doesn't matter much, because the phase factor is so tiny that it doesn't affect the 
integral. This is the coherent radiation. 

For k^co (high frequencies, low wavelengths) the phase will change with infinitely big 
steps. Even around a small dr we will ideally have all phases adding up in the integral, and these 
will eventually cancel out. Hence at this case where the phases are random the bunch form 

factor will approach and Ifotdi^ ^^)^^ h{^) • Here the electrons act completely at 
random, and they produce a total electric field that is only the square root of the electric field 
they produce when they act in harmony, and hence now the intensity scales linearly with the 
number of electrons and the charge of the bunch. This is the incoherent radiation. 

C-3 The Bunch Form Factor in a more Specific Case 

In order to evaluate the bunch form factor for a general frequency, let 
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This means that we will have a cylindrical coordinate system to describe the transverse 
dependence of the bunch which is assumed Gaussian. Since R-r = p' ^m6co^{(^ - (j)'^ + z co^O , 
then the BFF reads 
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The inside integral can be calculated as follows: 

2? , , <i>-<i>'=w i 

j^..Vsin^cos(^-^)^^, = + j e^''''''^'''''''dii/ = 2nJ,{kp'^me) (135) 
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Then using also the identity \e~"''' J {bt)f^^dt = -e **"' [4], the bunch form factor can 



be written as 



f{k;e) = e 
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i-oo 



(136) 



It is clearly seen that the bunch form factor is directly proportional to the squared amplitude of 
the Fourier transform of the longitudinal distribution of the bunch. Thus, the geometrical 
characteristics of the electron distribution are mapped onto the spectrum of the emitted 
radiation. 

C.4 CTR Interferometry 

Coherent Transition Radiation is emitted when an electron beam hits a foil, and it is emitted 
near a wavelength equal to the dimensions of the beam itself. Since the detectors typically are 
not capable of resolving the picoseconds-long CTR signal in time, interferometry can be used to 
retrieve time-dependent information of the signal. The electric field is added to a slightly 
delayed in time version of the same electric field, and the intensity sum is observed, usually time 
integrated. The output yields information of the shape of the beam, and it is especially useful 
when more than one bunch of electrons are present. 

Let's start with the electric field of the CTR: E(t) . This is a spatially integrated time- 
dependent magnitude of the full field E(^r,(/>,z,t) . At the exit of the interferometer, this field is 
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added to a shifted in time version of the field, E(t-T). Here T is the shift in time, which is 
controlled by adjusting the lengths of the arms of the interferometer. Negative values of T 
imply that the second field comes earlier in time than the first, while positive values correspond 
to an electric field that comes later in time than the original. For a beam that is about 1 mm 
long, the time shift should be controlled in the 1 mm range, which is about 3 ps. At the exit of 
the interferometer the total field is E{t) + E{t-z). 

This field will hit a CTR detector, the same one that was used to measure the plain CTR. The 
detector will yield a time-integrated value of the intensity of the input fields, which will be 
proportional to 

+CO +00 

I{t) = J [E{t) + E(t - T)f dt = 2lQ+2J E{t) ■ E{t - T)dt (137) 

—CO —00 

+CX) 

Where j [£"(0] dt^I^. Hence the output will be a constant background enhanced by the 

—00 

autocorrelation integral of the beam with itself. For a single Gaussian bunch, i.e. 
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E{t) = E^e '^^ (138) 



The recorded interferometer CTR signal is 
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Here /q = EI^ttg^ . It is therefore another Gaussian, only with a standard deviation that is ^/2 
times larger now than the standard deviation of the Gaussian bunch itself. The pattern has a 
minimum at 2/^ and it reaches double this background value at the best case ( r = ) which is 
when the bunch overlaps with itself. By measuring the width of the CTR Gaussian, one can 
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therefore measure the width of the initial bunch. An absolute measurement of the intensity will 
also yield the strength of the bunch. 

Let's assume now that there are two distinct electron bunches that comprise the beam, 
each at different spots along time and with different amplitudes. We have 

E,(t) = E,,e '''- E,(t) = E,,e '"- (140) 

E{t) = E,(t) + E,(t) 

The total signal recorder at the detector will be (after some analytical manipulations) 

(^01 "^02) 
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There are 3 groups of terms in the result. The first group consists of the first 3 terms, and it is 
just the constant background radiation interference. The second group consists of the next 2 
terms, which are the self-correlation terms for each bunch. The third group consists of the last 2 
terms, and they are two satellite lobes, each corresponding to a cross term between the first 
and second bunch and vice versa. It is clear that the location and shape of each satellite lobe can 
provide information for both bunches. 



168 



Appendix D 

ThemOsiris Simulation Code 

The fully explicit particle-in-cell code OSIRIS [71] is widely used for calculating the beam-plasma 
and laser-plasma interactions in a plasma accelerator. However, one of the drawbacks of the 
code is that only a single plasma density can be simulated at any given run. For the multibunch 
ideas and experiments presented in this dissertation though, it is desirable to get a quick answer 
regarding the evolution of the wakefield amplitude with respect to the phase of the bunches for 
different plasma densities. For this reason a simulation code was developed that solves for the 
wakefield amplitude given an input beam profile. The mathematical analysis behind the code is 
briefly presented in Section D.l, while Section D.2 offers an example of an implementation of 
the code in Matlab. This code was used to generate most of the non-OSIRIS simulation figures in 
the work presented in this dissertation. 

D.l Analytical Formulation 

The starting point are Maxwell's equations in a linear medium, explicitly written in equation (73) 
. Instead of assuming linear approximation, we solve the equations in the fully relativistic 
nonlinear ID regime, along with the relativistic Newton's law of motion, equation (74), and the 
continuity condition for electrons, equation (75). Therefore, assuming d/dr = (and thus no 

magnetic field) and dependence on the variable (^ = z-ct, then V = Z'd I d^ and we get 
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The electric field E and the electron velocity v^ can be substituted in the above system of 
equations. After some analytical manipulation, a differential equation for the total electron 
density n^ as a function of the length (^ can be written as follows: 
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This equation is equivalent to the one derived by Rosenzweig in [146]. The electric field E is 
then found from 
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(144) 



The above equations intrinsically include nonlinear effects, such as the sawtooth-shaped 
wakefields and the frequency shifts that emerge in the nonlinear regime. However the 
equations are only valid as long as the plasma can be treated as a fluid, and this does not always 
hold in the blowout regime due to particle trajectory crossings [116]. Still, these equations can 
be used to obtain rough wakefield amplitudes and bunch phasing in the mildly nonlinear regime. 
Note that in the weakly linear regime we have nj nQ<^l and by Taylor-expanding n^ 
around n^ the above equations are approximated by 
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Comparison between OSIRIS and ThemOsiris codes 
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Figure D.l: Comparison between OSIRIS and ThennOsiris codes for the sanne set of input paranneters. 
Three electron bunches with peak beann density of 2x10^^ cnn'^ are fed into a 5x10^^ cnn'^ plasnna. The 
deviation between the codes is less than 5% in ternns of the wakefield evolution, while ThennOsiris is at 
least 1,000 tinnes faster. 
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These are the results of the linear wakefield theory that was presented in Appendix A.l . In 
order to obtain a wakefield amplitude that includes 2D transverse effects, the amplitude found 

through equation (144) can be multiplied (for transversely Gaussian bunches) by /?(0) which is 

defined in equation (9). 

A comparison between the outputs of ThemOsiris and OSIRIS is presented in Figure D. 
An implementation of ThemOsiris in Matlab was written. The input beam electron density 
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consists of 3 bunches, the last one separated by many plasma wavelengths. We observe that as 
benchmarked by OSIRIS the code correctly predicts the wakefield amplitude with less than 5% 
error. Being -1000 times faster than OSIRIS, the code can be used to do quick plasma density 
scans and/or to determine the proper position of the bunches. 

In the second revision of this code, the ability to generate the beam energy spectra at the 
exit of the plasma was added (including intrinsic energy spread). The code is limited to slow 
variation of the initial neutral plasma density (on a scale longer than the beam length) and 
assumes that the beam does not evolve with time. 

D.2 Sample Code 

An example of ThemOsiris code than can run in Matlab is given below. The code outputs the 
wakefield vs. time as well as the predicted energy spectra of given bunches. The functions that 
are called are also explicitly written below. 



% ThemOsiris 

% This is a code that will calculate the phase space and energy gain of 
% 2D bunches into a high density plasma, used for our Brookhaven 
% ATF experiment. . . 

clear all; 

% Useful physical constants 

c = 3e8; 

e = 1.6e-19; 

m = 9.1e-31; 

% Plasma normalization parameters 

lamdaO = 250e-6; % Reference wavelength for normalizations 

k0 = 2*pi/lamda0; % Reference wavevector 

n0 = Ie-6*lamdap2n0(lamda0); % Reference plasma density in cm^-3 

wp = 2*pi*c/lamda0; % Reference plasma frequency in Hz 



% INPUT DECK! ! ! 

% 

% 



% Plasma Density 
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% 

np=1.79el6; % in cm^-3 

% Input parameters for the 2 bunches 
% 

% Energy of the bunches in MeV 
EnergyO = [56 57 58 59 60 61]; 

% Total Charge for drive bunches [pC] 
qt=250; 

% How many drive bunches? 
Nb=4; 

% How much of the charge is *not* blocked? [l=none] 
b=0.5; 

% Charge of each bunch in pC 

% (the drive bunch has half the charge) 

q = b*[qt/Nb*l qt/Nb*3 qt/Nb*5 qt/Nb*7 qt/Nb*0.001 qt/4/2]; 

% Bunch seperation distance in microns 
% (relates to plasma density) 
sep = n021amdap(np)*le6; 

% Position of each bunch in sec 
%mean = [250 250+1. 5*sep 250+3*sep 250+4. 5*sep 1250 1625+125]*le-6/c; 
mean = [250 250+1. 5*sep 250+2*1. 5*sep 250+3*1. 5*sep 250+6*sep 1562.5]*le-6/c; 

% sigma_z of each bunch in sec 

% (such that fraction b goes through) 

sigma_z = b/2*[sep sep sep sep sep sep/4/b]*le-6/c; 

% Spot size of each bunch in microns 
sigma_r = 100*[1 11111]; 



% Define the spatial grid for the simulation 
% 

zmin = 0.0*k0*le-3; % Start of simulation box in mm 

zmax = 20.0*k0*le-3; % End of simulation box in mm 

N = 100; % Grid number of points for simulation box 

z = linspace(zmin,zmax,N); % The grid in propagation direction 

% 

% 

% Solve for the wakefield generated from the above bunch distribution 

% tspan in ps, E in MeV/m, nb in m^-3 

[tspan, E, nb] = multi_gaus_wake(np, mean, sigma_z, sigma_r, q); 

% INITIALIZE PARTICLE POSITION AND MOMENTUM!!! 
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%- 
%- 



NN=length(tspan); 

% Calculate the beam densities seperately 



for i=l 

nbl 
sigma_r(l)) 

nb2 
sigma_r(2)) 

nb3 
sigma_r(3)) 

nb4 
sigma_r(4)) 

nb5 
sigma_r(5)) 

nb6 
sigma_r(6)) 
end 



NN 

i) = squane(tspan(i)*le-12, mean(l), sigma_z(l), 3, q(l)^ 

i) = square(tspan(i)*le-12, mean(2), sigma_z(2), 3, q(2), 

i) = square(tspan(i)*le-12, mean(3), sigma_z(3), 3, q(3), 

i) = square(tspan(i)*le-12, mean(4), sigma_z(4), 3, q(4), 

i) = square(tspan(i)*le-12, mean(5), sigma_z(5), 3, q(5), 

i) = square(tspan(i)*le-12, mean(6), sigma_z(6), 3, q(6). 



% How many particles do you want in the simulation total??? 
Np = 25000*length(q(q>l)); 



Npl = q(l)/sum(q)*Np; 
Np2 = q(2)/sum(q)*Np; 
Np3 = q(3)/sum(q)*Np; 
Np4 = q(4)/sum(q)*Np; 
Np5 = q(5)/sum(q)*Np; 
Np6 = q(6)/sum(q)*Np; 



% Number of particles in each bunch 
% Number of particles in each bunch 
% Number of particles in each bunch 
% Number of particles in each bunch 
% Number of particles in each bunch 
% Number of particles in each bunch 



% Renormalize density so that we have this number of particles 

nbl = nbl.*Npl.*length(nbl)/(tspan(end)-tspan(l))/sum(nbl); 
nb2 = nb2.*Np2.*length(nb2)/(tspan(end)-tspan(l))/sum(nb2); 
nb3 = nb3.*Np3.*length(nb3)/(tspan(end)-tspan(l))/sum(nb3); 
nb4 = nb4.*Np4.*length(nb4)/(tspan(end)-tspan(l))/sum(nb4); 
nb5 = nb5.*Np5.*length(nb5)/(tspan(end)-tspan(l))/sum(nb5); 
nb6 = nb6.*Np6.*length(nb6)/(tspan(end)-tspan(l))/sum(nb6); 



% Create the initial particle location matrix 
% 

dz = tspan(2)-tspan(l); 

zel=[]; ze2=[]; ze3=[]; ze4=[]; ze5=[]; ze6=[]; 

% Initialize the particles location matrix 

for i=2:length(tspan); 
interval dz (dt) 

num_parl = round(nbl(i)*dz); 
particles are in that interval 
num_par2 = round(nb2(i)*dz); 



% Elementary step 

% for every time 
% This is how many 



174 



num_par3 = round(nb3(i)*dz); 
num_par4 = round(nb4(i)*dz); 
num_par5 = round(nb5(i)*dz); 
num_par6 = round(nb6(i)*dz); 



tempzl = linspace(tspan(i-l),tspan(i),num_panl); 
the particles oven that interval 

tempz2 = linspace(tspan(i-l),tspan(i),num_par2); 
tempzB = linspace(tspan(i-l),tspan(i),num_par3); 
tempz4 = linspace(tspan(i-l),tspan(i),num_par4); 
tempzB = linspace(tspan(i-l),tspan(i),num_par5); 
tempz6 = linspace(tspan(i-l),tspan(i),num_par6); 

tempzl = tempzl(2:end); 
overlapping particle 

tempz2 = tempz2(2:end); 
tempzB = tempz3(2:end); 
tempz4 = tempz4(2:end); 
tempzB = tempz5(2:end); 
tempz6 = tempz6(2:end); 



% Evenly distribute 



% Ignore the first 



zel = [zel tempzl]; 
to the location matrix 

ze2 = [ze2 tempz2]; 

ze3 = [ze3 tempzB]; 

ze4 = [ze4 tempz4]; 

ze5 = [ze5 tempzB]; 

ze6 = [ze6 tempz6]; 



% Add those particles 



end 



Npl=length(zel); 
Np2=length(ze2); 
Np3=length(ze3); 
Np4=length(ze4); 
Np5=length(ze5); 
Np6=length(ze6); 



% actual number of particles in beam 1 
% actual number of particles in beam 2 



Np = Npl+Np2+Np3+Np4+Np5+Np6; % Total actual number of 

particles 

ze = [zel ze2 ze3 ze4 ze5 ze6]; % Merging the two beams 

together! 

% Now, the matrix ze contains the positions of the particles with respect 

% to the beggining (head) of the beam 



% Create the initial particle momentum matrix 
% 

s_th=0.1; 

Energyl(l:Npl) = Energy0(l) + s_th*randn(l,Npl); 

beam particles 



% Thermal spread in MeV 
% Energy of the first 
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Energy2(l:Np2) = Energy0(2) + s_th*randn(l,Np2); % Energy of the second 

beam particles 

Energy3(l:Np3) = Energy0(3) + s_th*randn(l,Np3); 

Energy4(l:Np4) = Energy0(4) + s_th*randn(l,Np4); 

Energy5(l:Np5) = Energy0(5) + s_th*randn(l,Np5); 

Energy6(l:Np6) = Energy0(6) + s_th*randn(l,Np6); 

Energy = [Energyl Energy2 Energy3 Energy4 EnergyB Energy6]; 

% Merge the two energy distributions 



% Modify the start and end of simulation box is necessary 
loc_z_start = 1; 
loc_z_final = length(z); 
z_start = z(loc_z_start); 
z_final = z(loc_z_final); 

% Create the matrix that holds the momentums of the particles 

p = zeros(loc_z_final-loc_z_start+l,Np); 

% Initialize the beam electrons at their initial momentum[normalized units] 

p(l,:) = sqrt((l+Energy(:)/0. 511875). ^2-1); 

% Now matrix p holds the momentums of all the Np particles of both bunches 

% Find each particle location in tspan 

par_loc=round(l+((NN-l)/(tspan(end)-tspan(l)))*ze); 

par_loc(par_loc>5000)=5000; % Ignore further away particles 

% This can be done because tspan(i)=1000/(NN-l)*(i-l)=linspace(0,1000,N); 

% 

% 



% Main loop for pushing particles 

% 

for i=loc_z_start:loc_z_final-l % For each step along z 

% Calculate next step for the beam 
dz = z(i+l)-z(i); 

for j=l:Np % For each particle along the beam 
% And now the main push! 
p(i+l-loc_z_start+l, j)= p(i-loc_z_start+l, j) +... 
dz/k0*E(par_loc(j))*le6*e/c^2/m; 
end 
end 

% Translate momentum to energy in MeV 
Energy = 0. 511875. *(sqrt(p.^2+l)-l); 

% Calculate Transformer Ratio 

% 

% Number of driving bunches 
Nd = length(q(q>l))-l; 
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for i=l:Nd 

% Extend of bunch region in psec 

L_low = (mean(i)-sigma_z(i))*lel2; 

L_high = (mean(i)+sigma_z(i))*lel2; 

t_low = length(tspan(tspan<L_low)); % Position in tspan 

t_high = length(tspan(tspan<L_high)); % Position in tspan 

% Maximum decelerating field inside that bunch 

Em(i) = abs(min(E(t_low:t_high))); 
end 

% Maximum decelerating field inside the drive bunches 
Emin = max(Em); 

% Up until the beginning of the drive bunch 

Lb = (mean(end)-sigma_z(end))*lel2; 

tb = length(tspan(tspan<Lb)); % Position in tspan 

% Maximum accelerating field after the drive bunches 

Emax = max(E(t_high:tb)); 

% Transformer Ratio 
R = Emax/Emin; 



% PLOTTING PARAMETERS AND FUNCTIONS 

% 

% 



% Define Bin Centers 

E_min=53; 

E_max=68; 

M=linspace( E_min, E_max, 200) ; 

figure( 'Position', [100 100 800 700]); 

subplot(2,l,l); 

hold on; 

area(tspan, (nbl+nb2+nb3+nb4+nb5)/max(nbl+nb2+nb3+nb4+nb5+nb6)*180, 'facecolor ' , 

'green'); 

area(tspan,nb6/max(nbl+nb2+nb3+nb4+nb5+nb6)*180, 'Facecolor', [1 0.5 0]); 

plot (t span, E, 'linewidth' ,3); 

grid on; 

title([' Wakefield and Beam Density, n_p=' ,num2str(np, '%2.1e' ), 'cm^-^3' ]) 

ylabel( 'Wakefield [MV/m]'); 

xlabel('Time [ps]'); 

axis([0 7 -200 200]); 

legend( 'Drive Beam' , 'Witness' , 'Wakefield' , 'location' , 'southwest' ); 

text ( 5 . 1, 175 , [ ' R= ' , num2st r ( R, ' %1 . 3g ' ) ] ) ; 

enhance_plot; 

[N0,X0] = hist(Energy(l,:),M); 

[N1,X1] = hist(Energy(end,l:Npl),M); 

[N2,X2] = hist(Energy(end,Npl+l:Npl+Np2),M); 

[N3,X3] = hist(Energy(end,Npl+Np2+l:Npl+Np2+Np3),M); 

[N4,X4] = hist(Energy(end,Npl+Np2+Np3+l:Npl+Np2+Np3+Np4),M); 
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[N5,X5] = hist (Energy (end, Npl+Np2+Np3+Np4+l:Npl+Np2+Np3+Np4+Np5),M); 
[N6,X6] = hist(Energy(end,Npl+Np2+Np3+Np4+Np5+l:Npl+Np2+Np3+Np4+Np5+Np6),|vi); 
%[N7,X7] = hist(Energy(end,Npl+Np2+Np3+Np4+Np5+Np6+l:end),M); 
[N ,X ] = hist(Energy(end,:),|vi); 

subplot(2,l,2); 
area(Xl,Nl/Facecolor',[0 1]); 

title([ "Energy Spectrum after ' ,num2str(zmax/k0*1000)/mm of plasma']); 
ylabel('# of sim particles'); 
xlabel( 'eBeam Energy [MeV]'); 
axis([E_min E_max 10000]); 

hold on; 

area(X2,N2/Facecolor',[0 0.85]); 

area(X3,N3/Facecolor' J0 0.60]); 

area(X4,N4/Facecolor' J0 0.45]); 

area(X5,N5/Facecolor' J0 0.30]); 

%area(X6,N6/Facecolor' J0 0.15]); 

area(X6,N6, 'Facecolor' , 'magenta' ); 

%area(X7,N7, 'Facecolor' , 'magenta' ); 

plot(X0,N0, 'black: '); 

plot(X,N, 'r', 'linewidth',3); 

%legend( 'Bunch I'/Bunch 2'/Bunch 3'/Bunch 4'/Bunch 5'/Bunch e'/Witness 

Bunch' , 'Initial Spectrum' , 'Final Spectrum' , 'location' , 'northwest' ); 

enhance_plot; 



% Solves the ID nonlinear plasma differential equation 

% np in cm^-3 

% sigmal, sigma2,in s 

% sigma_r in microns, charge in pC 



function [t, E, nb] = multi_gaus_wake(np, mean, sigma_z, sigma_r, charge) 



c=3e8; 

e=1.6e-19; 

m=9.1e-31; 

wp = 2*pi*c/n021amdap(np); % in s^-1 



tmin = 0; 

tmax = 12; 

NN=5000; 

field 

tspan = linspace(tmin,tmax,NN); 



% start field time in 

% End field time in ps 
% Number of points for 

% in ps 
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tspan2 = tspan*le-12*wp; 
sim units 



% Time grid points in 



options=odeset('reltol',le-3, 'abstol' ,le-10); 



% Solver Options 



[t,v]=ode45(@rhs_multi_gaus,tspan2, [1, 0], options, mean, sigma_z, 
sigma_r, charge, np); 

n = v(:,l); 
density 

E = (v(:,2)./(2*n-l).-(1.5))'; 



sr=sigma_r(l)*le-6*2*pi/n021amdap(np); 
transverse radius 

E = E*sr.^2/2*exp(sr.^2/2).*expint(sr.^2/2); 
transverse effects 



% Electron Beam 
% Electric Field 
% Normalized 
% Account for 



t = t/wp*lel2; 
returned 



% time in ps is 



for i=l:size(t,l) nb(i)=square(t(i)*le-12, mean(l), sigma_z(l), 3, charge(l), 
sigma_r(l)) +. . . 

square(t(i)*le-12, mean(2), sigma_z(2), 3, charge(2), sigma_r(2)); end 
% nb is in m^-3 
t=f; 

Enorm = m*c*wp/e; 

E2 = E*Enorm/le6; % Return E in MeV/m 
E=E2; 



% The function that contains the right hand side of the differential equation 

to be solved 

% np is plasma density in cm^-3 



function out = rhs_multi_gaus(t, y, mean, sigma_z, sigma_r, charge, np) 



n = y(l); 
u = y(2); 

wp = l/n021amdap(np)*2*pi*3e8; 
np = np*le6; 



% plasma frequency in s^-1 
% in m^-3 



% Calculate beam density 
for i=l:length(mean) 

temp(i) = square(t/wp, mean(i), sigma_z(i), 3, charge(i), sigma_r(i)); 
end 
beam = sum(temp); 

outl = u; 
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out2 = 3/(2*n-l)*u^2 - (2*n-l)^(1.5)*(n - 1 + beam/np); 
out = [outl; out2]; 



% Returns the value of the gaussian at point x in sec 

% mean is the mean value in sec 

% Sigma is the Half Width of the square in sec 

% num is the numbers of sigmas on each side after which the gaussian is 

% charge is the charge in pC 

% sigma_r is the rms spot size in microns 

% 

% the output is the beam density at the point x in particles per m^3 

function out = square(x, mean, sigma_z, num, charge, sigma_r) 

charge = charge*le-12; 
sigma_r = sigma_r*le-6; 
e=1.6e-19; 
c=3e8; 

if abs(x-mean)<=sigma_z 

out = charge/(e*2*pi*sigma_r^2*2*sigma_z*c); 
else out=0; 
end 



% lamdap in m, n0 in cm^-3 

function out = no21amdap(n0) 

c=3e8; 
e=1.6e-19; 
m=9.1e-31; 
e0=8.85e-12; 

out = 2*pi*c/e*sqrt(e0*m/n0*le-6); 



function enhance_plot (fontname, font size, linewid, ma rkersiz) 

% enhance_plot([fontname,fontsize]); 

% 

% Function to enhance MATLAB's lousy text choices on plots. Sets the 
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% current figure's Xlabel, Ylabel, Title, and all Text on plots, plus 

% the axes-labels to the "fontname" and "fontsize" input here where 

% the defaults have been set to 'times' and 16. 

% Also sets all plotted lines to "linewid" and all markers to size 

% "markersiz". The defaults are 2 and 8. 

% 

% INPUTS: fontname: (Optional, DEF='TIMES' ) FontName string to use 

% MATLAB's ugly default is 'Helvetica' 

% fontsize: (Optional, DEF=16) FontSize integer to use 

% MATLAB's tiny default is 10 

% linewid: (Optional, DEF=2) LineWidth integer to use 

% MATLAB's skinny default is 0.5 

% markersiz: (Optional, DEF=8) MarkerSize integer to use 

% MATLAB's squinty default is 6 

% for all inputs, if pass 0, use default 

% 

% Modifications 

% 19-Feb-2002 3. Nelson 

% added linewid and markersiz to help squinting readers 

% 

% 20-Feb-2002 3. Nelson 

% added check for legend. If legend exists, increase the 

% line and marker size, also increase the font to 

% fontsize-2 (2 points smaller than title and labels) 

if ('-exist ( 'fontname' ) | (fontname==0)) 

fontname = 'times' ; 
end 
if ('-exist( 'fontsize' ) | (fontsize==0)) 

fontsize = 16; 
end 
if ('-exist ( 'linewid ' ) | (linewid==0)) 

linewid=3; 
end 
if ('-exist( 'markersiz' ) | (markersiz==0)) 

markersiz = 8; 
end 

Hf=gcf; 
Ha=gca; 

Hx=get(Ha, 'XLabel'); 
Hy=get(Ha, 'YLabel'); 
Ht=get(Ha, 'Title'); 
set(Ha, 'LineWidth ',.75); 
set(Hx, 'fontname' , fontname); 
set(Hx, 'fontsize' , fontsize); 
set(Hy, 'fontname' , fontname); 
set(Hy, 'fontsize' , fontsize); 
set(Ha, 'fontname' , fontname); 
set(Ha, 'fontsize' , fontsize); 
%set (Ha, ' YaxisLocation ' , ' right ' ) 
%set (Ha, 'YaxisLocation ' , ' left ' ) 
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set(Ht, 'fontname' ^fontname); 

set(Ht, 'fontsize' ^fontsize); 

set(Hy, 'VerticalAlignment', 'bottom' ); 

set(Hx, 'VerticalAlignment', 'cap' ); 

set(Ht, 'VerticalAlignment', 'baseline' ); 

Hn = get(Ha, 'Children'); 

n = length(Hn); 

if n > 

typ = get(Hn/Type'); 
for j = l:n 

if strcmp( 'text ' ,typ(j, : )) 

set(Hn(j), 'font name', font name); 
set(Hn(j), 'font size', font size); 
end 
if strcmp( 'line' ,typ(j, : )) 

set(Hn(j), 'LineWidth'^linewid); 
set(Hn(j), 'MarkerSize' ,markersiz); 
end 
end 
end 

legh=legend; 

Hn=get(legh, 'Children'); 
n = length(Hn); 
if n > 

typ = get(Hn/Type'); 
for j = l:n 

if strcmp( 'text ' ,typ(j, : )) 

set(Hn(j), 'font name', font name); 
set(Hn(j), 'fontsize', fontsize -2); 
end 
if strcmp( 'line' ,typ(j, : )) 

set(Hn(j), 'LinelAlidth',linewid); 
set(Hn(j), 'MarkerSize' ,markersiz); 
end 
end 
end 

figure(Hf ); 
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